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A COMPREHENSIVE ANALYTICAL MODEL OF 
ROTORGRAFT AERODYNAMICS AND DYNAMICS 


Part I: Analysis Development 

Wayne Johnson 

Ames Research Center 
and 

Aeromechanics Laboratory 
AVRADCOM Research and Technology Laboratories 

SUMMARY 

The development of a comprehensive analytical model of rotorcraft 
aerodynamics and dynamics is presented. This analysis is designed to calcu- 
late rotor performance, loads, and noise; helicopter vibration and gust 
response; flight dynamics and handling qualities; and system aeroelastic 
stability. The analysis is a combination of structural, inertial, and aero- 
dynamic models that is applicable to a wide range of problems and a wide 
class of vehicles. The analysis is intended for use in the design, testing, 
and evaluation of rotors and rotorcraft, and to be a basis for further 
develoi^ent of rotary wing theories. The analysis is implemented in a 
digital computer prograjn. 


1 . INTRODUCTION 

For the design, testing, and evaluation of rotors and rotorcraft, a 
reliable and efficient analysis of the aircraft aerodynamics and dynamics 
is required. It is necessary to predict and explain the rotor performance, 
loads, and noise; helicopter vibration and gust response; flight dynamics 
and handling qualities; and system aeroelastic stability. Such capability 
is also required as a basis for further development of rotary wing theory. 
This report presents the development of a comprehensive analytical model of 
rotorcraft aerodynamics and dynamics. 
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The analysis developed here is a consistent combination of structural, 
inertial, and aerodynamic models, applicable to a wide range of problems and 
a wide class of vehicles. Typically rotary wing analyses have been developed 
or verifie’ for only a particular type of helicopter or a particular techni- 
cal problem, that reflects the specific interests of the originating organi- 
zation, Tite present model is applicable to articulated, hingeless, giniballed, 
and teetering rotors with an arbitrary number of blades. The rotor degrees 
of freedom Included are blade flap/lag bending, rigid pitch and elastic 
torsion, and optionally gimbal or teeter motion. This analysis is applicable 
to general two-rotor aircraft, including single main-rotor and tandem heli- 
copter configurations and side-by-side or tilting proprotor aircraft con- 
figurations (fig. 1) . The case of a rotor or helicopter in a wind tunnel is 
also covered. The aircraft degrees of freedom included are the six rigid 
body motions, elastic airframe motions, and the rotor/engine speed perturba- 
tions. The trim operating conditions considered include level flight, steady 
climb or descent, and steady turns. The analysis of the rotor Includes non- 
linear inertial and aerodynamic models, applicable to large blade pitch 
an^ ts and high inflow ratio. The rotor aerodynamic model is based on two- 
dimensional steady airfoil characteristics with corrections for three- 
dimensional and unsteady flow effects, including a dynamic stall model. A 
detailed wake model for the rotor nonuniform inflow calculation is developed, 
with a lifting surface theory correction for vortex- induced loads. Available 
prescribed and free-wake- geometry models are used. The aeroelastic stability 
analysis derives linearized equations consistent with the nonlinear rotor 
model. 

The solution of the equations of motion is separated into two parts, 
based on the different time scales involved in rotorcraft dynamics. The 
first part is the solution for the rotor motion and the airframe vibration. 
This motion is periodic, with fundamental frequency fi for the rotor and Ni] 
for the airframe (n is the rotor rotational speed and N is the number of 
blades). The periodic motion is calculated by a harmonic analysis method. 

The second part is the solution for the steady state or slowly varying air- 
frame motion (consisting of the aircraft rigid body and rotor speed perturba- 
tions, and the static elastic deflection of the airframe and drive train) . 
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The assumption that the aircraft motion is quasi-static (compared to the rotor 
speed) allows the periodic rotor solution to be used for transient motions of 
the helicopter as well as for the trim calculations. Most importantly, by 
taking adva :age of the frequency separation of the rotor and aircraft motions, 
an economical calculation procedure is realized. 

The first computation task is the trim analysis, in which the control 
position and aircraft orientation are determined for the specified operating 
condition. The periodic blade motion is calculated, and then the rotor per- 
formance, loads, and noise can be evaluated. The rotor model in the trim 
solution can use uniform inflow, nonuniform Inflow with a rigid wake geometry, 
or nonuniform Inflow with a free wake geometry. The aeroelastic stability, 
flight dynamics, and transient analyses begin from the trim solution. The 
aeroelastic stability analysis sets up a set of linear differential equations 
describing the motion of the rotor and aircraft; the eigenvalues of these 
equations define the system stability. The flight dynamics analysis calcu- 
:..tes the rotor and airframe stability derivatives, and sets up linear differ- 
eati.l equations for the aircraft rigid body motions; the poles, zeros, and 
eigenvectors of these equations define the aircraft flying qualities. The 
transient analysis numerically integrates the rigid body equations of motion 
for a prescribed control or gust input. 

In this analysis all quantities will be dimensionless, based on the air 
density p , the rotor radius R , and the rotor rotational speed fi. 

2. ROTOR MODEL 
2.1 Structural Analysis 

The rotor structural analysis consists of an engineering beam theory 
model for the coupled flap/lag bending and torsion of a rotor blade with large 
pitch and twist. A high aspect ratio (of the structural elements) is assumed, 
so the beam model is applicable. The objective is to relate the bending 
moments at the section, and the torsion moment, to the blade deflection and 
elastic torsion at that section. The analysis follows the work of reference 1. 
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LOCUS OF SECTION 
CENTER OF GRAVITY 
(CG) 


Figiare 2. Geometry of the undeformed blade. 



2.1.1 Geometry.- The basic assumptions are that an elastic axis exists, 
and the undeformed elastic axis is a straight line; and that the blade has a 
high aspect ratio (of the structural elements) so engineering beam theory 
applies. Figure 2 shows the geometry of the undeformed blade. The span 
variable r is measured from the center of rotation along the straight 
elastic axis. The section coordinates x and z are the principal axes of the 
section, with origin at the elastic axis. Then by definition, /(xz)dA = 0. 
Really this integral is over the tension carrying elements, i.e., a modulus 
weighted integral: /xzE dA = 0. This remark holds for all the section 

Integrals in the structural analysis. The tension center (modulus weighted 
centroid) is on the x axis, at a distance x aft of the elastic axis: 

/x dA = x^A and fz dA = 0. Again, these are modulus weighted integrals. 

If E is uniform over the section, then x^ is the area centroid; and if the 

section mass distribution is the same as the E distribution, then the ten- 

s’rr. center coincides with the section center of gravity. 

The angle of the major principal axis (the x axis) with respect to 

'.he hub plane is 6. The existence of the elastic axis means that twist about 

the elastic axis occurs without bending. In general, the elastic torsion 

deflection will be included in 6 . The blade pitch bearing is at the radial 

station r^^. The blade pitch is described by root pitch 6° (rigid pitch 

about the feathering axis, including that due to the elastic distortion of 

the control system), built-in twist 6^ , and elastic torsion about the 

tw 

elastic axis 0^. So 0=0° + 0^^ + 0^, where 0°(ij/) is the root pitch, 

8(rpA) = 6°; built-in twist, = 0; and 0^(r,4/) is the 

elastic torsion, = 0. There is shear stress in the blade due to 0 

only. It is assumed that 0 is small, but 0° and 6 are allowed to be 

® tw 

large angles . 

The unit vectors in the rotating hub plane axis system are "i" , ^ , and 

t , 5 B 

g (fig. 2). The unit vectors for the principal axes of the section (x, r, z) 
are 1 , j, and k; these are for no bending, but Include the elastic torsion in 
the pitch angle 0 . So the principal unit vectors are rotated by 0 from the 
hub plane: 
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2.1.2 Descvi^tion of the bending.- Now the engineering beam theory 

assumption is introduced: plane sections perpendicular to the elastic axis 

remain so after the bending of the blade. Figure 3 shows the geometry of the 

deformed section. The deformation of the blade is described by (a) deflection 

of the elastic axis, x , r , and z ; (b) rotation of the section due to bend- 

o o o 

ing, by <()^ and cf)^; and (c) twist about the elastic axis by 6^, which is 

implicit in 't' and it. The quantities x , r , z , , 6 , and 0 are assumed 

o’ o o Z e 

to be small. 

The unit vectors of the unbent cross section are 1*, The unit 

vectors of the deformed cross section are T s » and "k , where and 

xs ’ xs * xs xs 

k are the principal axes of the section, and r. is tangent to the 
xs xs 

deformed elastic axis. It follows that 

-t- 



Now by definition, j = d r/ds, where r = xT^+(r + r )i + z k and 
s is the arc length along the deformed elastic axis. Hence to first order 

f 

~ ^ ^ - 4 - 

= e’)i^ ^ ^ 


- 7 - 




Figure 3- Geometry of the deformed blade. 



It follows the rotation of the section is 



I / 

X o © 

/ / 


or 


The undeflected position of the blade element is r 
the deflected position is 


xf^ + r^ H- z^, and 



^ ^ -4- ?o*1^ + + X 4:; — ^ -f 


The first term in the deflected position is the radial station; the 

next three terms are the deflection of the elastic axis; the next term is the 

rotation of the section; and the final two terms are the location of the point 

on the cross section. For now the elastic extension r will be neglected. 

o 

The strain analysis is simplified since then to first order, s = r; r^ gives 
a uniform strain over the section, which may be reintroduced later. 


2.1.3 Analysis of strain.- The fundamental metric tensor g of the 


mn 


undistorted blade is defined by: 


dr 






^v, 
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Where ds is the differential length in the material, and x are general 

m 

curvilinear coordinates. Similarly, the metric tensor G of the deformed 

mn 

blade is 




( • ( 1 







Then the strain tensor is defined by the differential length increment 

Z ^vv,w C ^ f 

^vwv« ~ 2^ C 


For engineering beam theory, only the axial components of the strain and 
stress are required. For a full exposition of the analysis of strain, the 
reader is directed to reference 2. 


The metric of the undeformed blade (no bending, and no torsion so 

6* = ®'tw^ obtained from the undistorted position vector = xi^H- rj*+ zTc, 
giving 




ar br ' 




The metric of the deformed blade, including bending and torsion, is similarly 
obtained from the position vector 1? = (x + x )lf + (r + x<t> - z(^ )j^+ (z + z )lf: 

O Z X o 


- 10 - 



t~r 


^ ^ 

&<■ a r 

C \ 4- y ~ ft 4 ' 

4- b' ( t 
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Then the axial component of the strain tensor is 




f r 


z ^ ” *^^**'^ 

^Cl4-'m.^-i — ?,4*x 4-(xI,4-e(i4 

^'*- -i- ( %o -©'Cy + x.'^')^ — X* 3 


The linear strain (for small 


X , z , 6 , ()) , and (f) ) is 
o * o e X* z 




y. — 2 -*- -a^ (x xo 4 - 

^ («x!, -X fti -»- ©L 


The strain due to the blade tension, e^, is a constant such that the 
tension is given by the integral over the blade sectioni 


“T == 




Substituting for and using the results /zdA = 0, /xdA = x^A, and 

/(x^ + z^)dA = Ip = kp^A (where is the modulus weighted radius of gyra- 

tion about the elastic axis) , gives 
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In this expression, the strain due to the blade extension r has been 

o 

included. It follows that the strain may be written as 


+ e*e ^ 


2.1.4 Seation moments.- To find the moments on the section, the second 

engineering beam theory assumption is introduced: that all stresses except 

Orr are negligible. The axial stress is given by « ^^rr* direction 

of 0 is 
rr 



The moment on the deformed cross section (fig. 4) is M = M*^ +M*1^ 

. „ 7^ _ , X xs r-^xs 

+ “z^xs* moment about the elastic axis due to the elemented force o dA 

on the cross section is 

JM = (x t*i + X £ '> •iA 

Integrating over the blade section, there follows the result for the total 
moments due to bending and elastic torsion: 
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Figure 4. Bending and torsion moments on 
the blade section. 
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To has been added the torsion moment GJ6^, due to shear stresses pro- 

duced by elastic torsion. These moments are about the elastic axis. For 
bending it is more convenient to work with moments about the tension center 
at x^: 

/A^ = — ^ ^ 

/ V \:2 ~ \ 

Substituting for snd integrating, the moments are 

<*'?'' — e'c ex.,^ 

~ EXy, ( + ^!v.©e 

Mr 4 kf-r +e!^ E^l»p') ©L •+ 4~^ 
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where 


•:r ,2 = 

-X = \ 

X / 

= \ a-^:» i/\ - V? K 

-x^p = -O (n'^-4-*’^) JA 

^ x'‘+-e*'>^A 

Ipp — I) — Vf ')*’ ■iA 


The integrals are all over the tension carrying elements (i.e., modulus 
weighted). The tension T acts at the tension center at x ; hence the bend- 
ing moments about the elastic axis may be obtained from those about the tension 
center by + xj and 


Me 


( 2 ) 


2. 1,5 VectoT* forrnutation* - Define the section bending moment vector 
and the flap/lag deflection ^ as follows: 


M 


(1^ 

e 


W 


— X *r^ 


( 2 ) 


is not quite the moment on the section, because 
the and components of the moment) . The derivatives of ^ are 


M and M are really 
X z 


xs 




C^L 

c +i -I- 4>x)"^ 
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Then the result for the bending and torsion moments can be written as follows: 

Mr = ^ V.» t-4-©;^ -+-©^V^-T- 


This is the result sought, the relation between the structural moments and 
deflections of the rotor blade. 


Writing the bending stiffness dyadic as 
the purposes of this paragraph neglecting the 


gives 



E-X. 




II 


W 


•sS'®*’ 


El - El "tl* + El kk, and for 
EIxp terms, 


nre;., 


In this form the result appears as a simple extension of the engineering beam 

theory result for uncoupled bending and torsion (the 6' =0 case). The 

tw 

vector form allows a simultaneous treatment of the coupled Inplane and out-of- 
plane bending of the blade, with considerable simplification of the equations 
as a consequence. 

This relation between the moments and deflections is a linearized result. 

Thus the vectors T* and k appearing in El and in are based on the trim 

pitch angle 6 = 0^ + 6 « The net torsion modulus is 

lw 
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where T - Q. pmdp is the centrifugal tension in the blade* For the 

elastic torsion stiffness characteristic of rotor blades, the GJ term 

2 

usually dominates* The k T term is only important near the root for blades 

^ 2 

which are very soft torsionally* The 8^^ Elpp term is only important for 
very soft, highly twisted blades. 

2.2 Inertia Analysis 

This section derives the inertia forces of a helicopter rotor blade. 

The blade motion considered includes coupled flap/lag bending (including the 
rigid modes if the blade is articulated), rigid pitch and elastic torsion, 
gimbal pitch and roll (which are dropped from the model for articulated and 
hingeless rotors) or teeter motion (for two-bladed rotors only) , and the 
rotational speed perturbation- The geometric model of the blade and hub 
includes precone, droop, and sweep; pitch bearing radial offset; feathering 
axis droop and sweep; and torque offset and gimbal undersling. 

2.2.1 Rotor* geometry,- Consider an N-bladed rotor, rotating at speed 
Q (fig. 5). The m-th blade (m = 1 to N) is at the azimuth location ~ ^ 

+ mAiJj, where Aiji = 2 ti /N and = f2t is the dimensionless time variable. 
Because for steady flight the blade motion is periodic, it is only necessary 

4 

to calculate the motion and forces of one of the blades. For this reference 
blade we choose that identified by m = N. The S coordinate system (T^, 
kg) is a nonrotating, inertial reference drame (fig. 5) . The S system 
coordinates are the rotor shaft axes when there is no hub motion. When the 
shaft moves however, due to the motion of the helicopter or the wind tunnel 
support, the S system remains fixed in space. The B system (^» ^ 

coordinate frame rotating with the m-th blade. The acceleration, angular 
velocity, and angular acceleration of the hub, and the forces and moments 
exerted by the rotor on the hub are defined in the nonrotating frame (the S 
system). Figure 6(a) shows the definition of the linear and angular motion of 
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(a) SHAFT MOTION 



(b) HUB REACTIONS 

Figure 6. Notation and sign conventions for the linear 
and angular shaft motion, and the forces and 
moments acting on the rotor huh. 


the rotor hub, and figure 6(b) shows the definition of the rotor forces and 
Bonents acting on the hub. The rotor blade equations of «,tlon will be 
derived in the rotating frame. 

Figure 7 shows the blade hub and root geometry considered. The origin 
of the B .,nd S system is the location of the gimbal (or teeter hinge). For 
articulaten or hingeless rotors, where there is no gimbal, this is simply the 
point where the shaft motion and hub forces are evaluated. The hub of the 
rotor is a distance below the gimbal (gimbal undersling, which Is not 

shown in fl3. 7 ). The torque offset Is positive In the -± direction. 

The arisuth is measured to the feathering aais line (Its projection in 

the hot plane), so the feathering axis Is parallel to the axis, and offset 
from the center of rotation. The precone angle 6p. , gives the orienta- 
tion of the blade elastic axis Inboard of the pitch bearing with respect to 
the hub plane; 6p^^, Is positive upward, and is assumed to be a small angle 
1-he pitch bearing is offset radially from the center of rotation by r n,e 
-.id pitch rotation of the blade about the feathering axis occurs at "r 
d.oop angle SfAj and sweep angle occur at just outboard of 

•h,. oitch bearing; and give the orientation of the elastic axis of 

the blade outboard of the pitch bearing, with respect to the precone. Both 
«FA2 and dpAj are assumed to be small angles; SpAj is positive downward, and 

6FA3 i= Positiva aft. Feathering axis droop and sweep 6 pa, define the 

orientation of the feathering axis with respect to the precone; is 

positive downward. 6FA5 is positive aft. and both are small angles! If 
= <SfA 5 " feathering axis orientation Is just given by the 

precone; if = 6pA2 and 6FA5 “ ^FA3 then the orientation is the same as 

the outboard elastic axis. 

From the root to the pitch bearing (at r = r^^) . the undistorted elastic 

axis is a straight line at the precone angle to the hub plane. The blade out- 

board of the pitch bearing has a straight undistorted elastic axis, with small 
droop and sweep angles. The feathering axis also has small droop and sweep 
with respect to the precone. The entire blade is flexible in bending. The 
portion of the blade outboard of the pitch bearing is flexible in torsion as 
well. The rotation of the blade about the pitch bearing takes place about the 
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Figirre 7. Schematic of the rotor huh and root geometry. 

Only a single, undistorted blade is shown, without 
the gimbal undersling. The gimbal is omitted from 
the model for articulated and hingeless rotors. 



/i 


*FA 



b) top view 
Figure ?. Concluded. 
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local direction of the feathering axis. Incorporation of bending flexibility 
of the blade inboard of the pitch bearing allows consideration of an articu- 
lated rotor with the feathering axis inboard or outboard of the hinges, or a 
cantilever blade with or without flexibility Inboard of the pitch bearing. 

Figure 8 shows the undeformed geometry of the blade. The description 
of the blade for the inertial analysis parallels that for the structural 
analysis (see fig. 2 and section 2.1.1). It is assumed that an elastic axis 
exists, and that the undeformed elastic axis is a straight line; and that the 
blade has a high aspect ratio. Here is the locus of the section center 

of gravity, x^ is the locus of the section aerodynamic center, and x^ is 
the locus of the section tension center. The distances x^^, x^, and x^ are 
positive aft, measured from the elastic axis; in general they are a function 
of r. The corresponding z displacements are neglected. 

The and"k coordinate system is the elastic axis/principal axis 

system of the section. The direction of the elastic axis is and are 

the direction of the local principal axes of the section. The spanwise 
variable is r, measured from the center of rotation. This variable is 
dimensionless, so r = 1 at the blade tip. The section coordinates x and z 
are mass principal axes, with origin at the elastic axis. It is assumed that 
the direction of the mass principal axes and the modulus principal axes is the 
same. The center of gravity is at z = 0 and x = x^. The section mass, 
center of gravity position, and section polar moment of inertia (about the 
elastic axis) are by definition then as follows: 

/ dm = m 
/z dm = 0 
/xz dm = 0 
/x dm = x^ m 
/(x^ + z) dm = Iq 

where the Integrals are over the blade cross section. 
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Figure 8. Geometry of the undeforraed blade. 



The droop and sweep of the blade elastic axis are defined with respect 
to the hub plane axes, so it follows that unless the feathering axis is 
parallel to the outboard elastic axis, these angles vary with the pitch of the 
blade. Let <5 fA2* '^FA 3 * droop and sweep of the blade when the 

pitch angle at 75% radius is zero. Then the following relation can be obtained 
from the root geometry: 


-f- (<^FA2. — 

■+- { Waj — 

-H — ^F*9 


where 6^^ is the blade pitch at 75% radius. 

2.2.2 Rotor motion.- The rotor blade motion is described by the follow- 
ing degrees of freedom: 

(a) Gimbal pitch and roll motion of the rotor disk (omitted for 
articulated and hingeless rotors) , or teeter motion of the blade 
(for two-bladed rotors only) . 

(b) Rotor speed perturbation. 

(c) Rigid pitch motion about the feathering axis and torsion about the 
elastic axis. 

(d) Bending deflection of the elastic axis, including rigid flap and 
lag motion of the blade is articulated. 

Figure 9(a) shows the gimbal motion and rotor speed perturbation in the non- 
rotating frame. The gimbal degrees of freedom are and respectively 
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(a) NONROTATING FRAME 


Figure 9 





(b) ROTATING FRAME 


. Notation and sign conventions for the 
motion and the rotor speed perturbati 



pitch and roll of the rotor disk. The rotor rotational speed perturbation is 
is Figure 9(b) shows the gimbal motion in the rotating frame. The 

degrees of freedom are B and 6 , given by 

(^Ca ~ pCarC 

0^ — -7- ^cc, -a- 

The blade pitch 0 is defined with respect to the hub plane, so only the 

blade inboard of the pitch bearing sees the pitch rotation due to 6 . For 

two-bladed rotors, the teetering degree of freedom may be included. The 

teetering motion is defined in the rotating frame, hence (^1)^ and 

0=0 for this case. 

G 

Figure 3 showed the geometry of the deformed blade. The blade deforma- 
tion is described by twist 6 about the elastic axis, bending deflection 
and z of the elastic axis; and rotations of the section by 4) and 4) due 
to the bending (see section 2.1.2). 

The blade pitch angle 0 is measured from the hub plane to the section 

major principal axis (the x axis). The undeformed pitch angle consists of 

the collective pitch ^qqII built-in twist define ^^qII 

as the pitch at r^. , so 6 (r_.) = 0. The rotation by 0 is not 

FA tw FA coll 

present inboard of the pitch bearing, but there can be pitch if the local 
principal axes with respect to the hub plane, which is included in 0^^ for 
r < r . The pitch of the deformed blade is composed of the root pitch 0^(4^) 

r i\. 

(the blade angle at the pitch bearing, r = rp^, due to control commands, 

control system flexibility, and kinematic coupling); the built-in twist 

and torsion about the elastic axis 6 (r,ij;) (where • only 

e e rA 

0^ produces shear stress in the blade). Thus the blade pitch is 

© » e° 
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e .. + e 

coll con 


Here 6 


coll 


The commanded root pitch angle is defined as 6* 
is the trim value of the collective pitch, which mayiriargrbut is steady^in 
time; and 6^^^ is the perturbation control input (including the cyclic con- 
trol requii :d to trim the rotor) , which is time dependent but is assumed to be 
a small angle. The blade root pitch commanded by the control system is 6^, 
while e° is the actual root pitch. The difference (0° - is the rigid 
pitch motion due to control system flexibility or kinematic coupling in the 
control system. Hence, the blade pitch may be written as 


= c ^ 




-rvT^ "+ ( 


The pitch angle may now be separated into trim and perturbation terms. 


e = 


6 + e, where the trim term is 

m 








•V-W 


and the perturbation is 




= (©'’-©*•') -f 


ccm 


The trim pitch 6^ is a large, steady angle; the perturbation pitch 
is small angle since all the components i 
root (r = is then 




small. 

The 

pitch 

at the blade 





o 


Ai* o 



-f- 

e 
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For the rigid pitch motion the notation is used: 


c 

fo ^ ^ 


Cer® — + ©OOV^ 


(This notation is consistent with that for the modal expansion of the elastic 
torsion 0^, as described below.) Note that p^ is the total rigid pitch 
motion of the blade, including the control angle O 

con 

2.2.3 Coordinate frames The rotating hub plane coordinate frame is 
obtained from the nonrotating hub plane frame by rotating about the z - axis: 

■4 = 

^6 = "^4 

"it = 


The blade coordinate frame is obtained by rotating by the angle 

“ ^FA 2 ^t>out the x-axis, by the angle the z-axis, 

and by the angle 0 about the y-axis : 






c-mrA d La — 

^3 




i - 


^^9 -f- 


-i. 




The cross section principal axes for the deformed blade are then 
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-h 

- a - ^4 


The vector is tangent to the deformed elastic axis. 

2.2*4 Blade acoclcvation,- The distance from the rotor hub to the 
center of gravity of the blade section is: 


^ L - Xp, _ _ c-Vj _ 

L- •+ '■f ft 3 

L <■ + (Sfc — 

-*-®° L 

- (».-e-x.<V)^^^ (r-r^,V 

— + (ifAj ~WAs^"^M^**fA'>3 



U L-VtA ~ (r~C^^') -r\ 

L - «fa - Cr-n»')ifA^ H- 1-( -I- iFA,') ■] 

L'- + «fa(5c 
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where 


for elastic torsion (k ^ 1) and 

y, = — 

■+- c«.-e-y,i^lf + 


for rigid pitch. Then the velocity of the blade section, relative to the 
rotating frame is 

\f - l-^^\ — -t* (f’n +^Ae<i^ -♦-'^ <■?& 

— ^ C t 

>« — ^ 

- £. f‘ 

Neglecting the squares of velocities, the acceleration relative to the rotat- 
ing frame is 

or = {£'^\ - -Ul ( «''^S -+ 3^* ©o') 

-f- ( - — -So^ 


For the blade Coriolis acceleration the radial velocity component *5^ • is 
required, including the effect of the change in the radial position of the 
section due to bending: 

-+■ ( + ifik, -i-s 

then ^ 

— ^6 -(. c _ (c- 

-*- tft . ( ^ 

^ vxi'>3 

The acceleration of the blade is required with respect to an inertial 
frame, specifically the S system. The B coordinate frame rotates at a con- 
stant angular velocity = nkg with respect to the S frame. The shaft 
motion is composed of linear and angular displacement of the origin of the 
S frame. The acceleration, angular velocity, and angular acceleration of the 
S system have the following components in the nonrotating, inertial frame: 
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— i 

— 

^ -J 


-3 

= 

■* . 


T^o = 


+ -^s 4- S/. 


t 

* t. 


It is assumed that a , Gi , and U) are all small quantities. 

o o o 

The acceleration (a^^) and velocity (y^) of the blade relative to the 
B frame have been derived above. Now the acceleration of a blade point in 
inertial space will be derived, in terms of the motion of the shaft, the rota- 
tion of the rotor, and the blade motion in the B frame. From the result for 
the acceleration in the rotating coordinate frame (the S frame, rotating at 


rate ), there follows: 
o 




-a 

r ) -j- (->„x r 


to 


where a and v are the acceleration and velocity relative to the 

r,s r,s ^ . 

S frame. The B system rotates at angular velocity ^ with respect 

the S frame. Hence with 9. constant and no angular or linear acceleration of 
the B frame with respect to the S frame, there follows: 

o^c^s = <^r 2. sc- X r 4* JTC X ( X r 3 


— V, 


+ ;^xc 


where 
Thus : 


a and v are the acceleration and velocity relative to the B frame, 
r r 

» = cxo-l-cxr + 2 -5c •+• -Sc X ( X c J 

4 2cixr 4- 2-io^ X 4- Uo -hva^xr 
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To first order in the velocity and angular velocity. 


this becomes finally: 





cxr -4- 2. X vr 


■P ^ 

^ X r 


X 






or in dyadic operator form, with n 




-I- ZS? ( -f>]4 


■c IX terms In a are respectively the acceleration of the origin, the 
Coriolis acceleration due to the angular velocity of the origin, the angular 
acceleration of the origin, the relative acceleration in the rotating frame, 
the relative Coriolis acceleration, and the centrifugal acceleration. 

For the blade bending and torsion equations, the following components of 
the acceleration will be required: 
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^ (:JxyK)fK 

Vt-o 

-+- SZ^ li ^ H- Wik, ~ 

-H 2^'^-^ -^r 

-h 252 ^ 

-f "2. 32 r Og ( c,«-aH> -4- o<^ 

-»- r- (' Ti ( S<^s-4» — -4T!^ 


For the hub moment the angular acceleration Is required: 

"r vi"^ s»^ -h -H c ( -^rcT^ — y«^ 

-4- ^ -S^ — £1. ( ^ 

H_fcO 

H- S£* <""ii tr •♦" *"^i ^ ^ -♦■ 
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H- 252.(cyt 

^ — (e- 

-+-(^-e— Xo 1 ?'> ^ 

-I- 2 -S^cT^j *^.Vr 

4-'mi r*irv — r (>J2^ v^NV ^ 

-+- Z: 5 ?_r^ U ( »<x - 4 - 

-V r^ (u ( ^ 5 wvv — ^ 

and for the hub force we can use 

« = -<-%ii» •^‘■ps.T^t -+- U,-€ y 

- ( <-% -^"^. (^u-t— *.-^V 

-s^.* L-e9TJ,.(x.-e+*.-i 

+‘fi* (-V* -+^''— “'■’^^^3 
-l-'^ ( Vs~.'V — 

+ ^ -H Cji, Sv~v\> 

-I- 2 SZ-r- (:-^g 
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The approximation r = has been used in all cases to evaluate the hub 

motion terms. 

2.2.5 Aerodynamic forces.- The aerodynamic forces acting on the blade 

section at the elastic axis are F , F , and F (see fig. 11 in section 2. A), 

z X r 

These are the components of the aerodynamic lift and drag forces in the hub 
plane axis system (the B frame) . F is in the hub plane, positive in the drag 
direction; F^ is normal to the hub plane, positive upward; and F^ is the 
radial force, positive outward. There are also radial components of F^ and 
F^ due to the tilt of the section by blade bending; here F^ is just the 
radial drag force. Thus the aerodynamic force acting on the section at the 
deformed elastic axis is: 

fr — f i f fc ■+• ifA, — + ) 

- Fx ( - ^4 -4- tfAj V ^ 


The section aerodynamic moment about the elastic axis is positive nose 

upward (so These section aerodynamic loads are Integrated 

over the blade span to obtain the total forces and moments. 

2.2.6 Force and moment equilibrium^- The equations of motion for 
elastic bending, torsion, and rigid pitch of the blade are obtained from 
equilibrium of inertial, aerodynamic, and elastic moments on the portion of 
the blade outboard of r: 

O ^ 

— AAg -4- Ma — 


where 
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Where is the structural moment on the inboard face of the deformed cross 

section (so -M_, is the external force on the outboard face) ; M is the 

A 

total aerodynamic moment on the blade surface outboard of r; and M is the 

total inert al moment of the blade outboard of r. The structural moment M 

E 

IS obtained from the engineering beam theory for bending and torsion (section 
2.1.5), from the control system flexibility for rigid pitch, or from the hub 
spring for gimbal or teeter motion. Alternatively, Mg may be viewed as the 
force or moment on the hub due to the rotor (so -Mg is the force on the 
rotor). Mj. is the inertial moment of the blade outboard of r, about the 
point r^(r), obtained by integrating the acceleration times the blade density 
(dm dp) over the volume of the blade: 

For bending of the blade, engineering beam theory gives 

'^E = -+ ) Me 

Therefore the operator Is applied to Mj and Tl also . For 

bending, moments about the tension center (x = x^) are required. Then the 

desired partial differential equation for bending is obtained from 

The ordinary differential equation for the k-th bending mode of one blade is 

obtained by operating with (...)dr, where is the flap/lag bending 

mode shape (see section 2.2.15). For elastic torsion, engineering beam theory 

gives “Jxs * operator is applied to m!|^ and “m^. For 

torsion, momn. Li, about the section elastic axis (x = 0) at r are required; 

also, elastic torsion involves only the blade outboard of r . The desired 

FA 

partial differential equation for torsion is then obtained from 3M /Sr. The 
ordinary differential equation for the k-th torsion mode is obtained by 
operating with /• 5j^(...)dr, where is the elastic torsion mode shape 

(see section 2.2.15). The equation of motion for the rigid pitch degree of 
freedom p^ is obtained from equilibrium of moments about the feathering axis, 

’ There M is the moment about the feathering axis (x = 0) 

- 38 - 


I 



at 


^ = ^FA’ 


is the direction of the feathering axis, including 
perturbations due to blade bending: 




r ^ 


The elastic restraint from the control system flexibility gives the restoring 
moment about the feathering axis, completing the desired equation of motion. 


The total rotor force and moment on the hub (at the gimbal point) are 
obtained from a sum over the N blades of and the force and 

moment due to the m-th blade: 













*^(m) (m) 

Since -F and -M are the forces on the blade, from force and moment 
equilibrium of the entire blade it follows that 


F 


(wn') 
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The hub force and moment are required In the nonrotating hub plane frame (the 
S system), with components defined as follows: 


■f = VA-e^ 

M = -Qis 


(see fig. 6). 


The equations of motion for the gimbal degrees of freedom and 6 

are obtained from equilibrium of moments about the gimbal, M =1* • “m and 

M is the total moment (from all N blades) about the 
gimbal point, in the nonrotating frame. The equation of motion for the teeter 
r’e?:ree of freedom is obtained from equilibrium of the moments about the 

-ee- r hinge from both blades, in the rotating frame. The equation of motion 
for the rotor speed perturbation is obtained from equilibrium of the 

shaft torque moments. Q = f ^ The drive train couples the torque 

perturbations of both the rotors, hence this degree of freedom is best con- 
sidered with the other motions of the helicopter body. 

2.2.7 Bending equation.- The equation of motion for blade bending is 
obtained from 








where M is 




the moment about the tension center (x = x ) at r, and 

C 

Cti!»c4 +“£ - Ctt M 
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Considering first the blade outboard of the inertia moinent is 


C ' r “1 

'c L — cc J 


X <x ^ ^ 


a m 4 j^ 


so 


25 ^ 


— -+- rv..-e -+- '^c'tV X ^ 


2 ' “1 
^ AA -?x^wv — 

ac*- ° ^ 


L ^ ‘ — >“c^ 


V A VV> 


1 " 


-? . M 




Finally: 


h AAt 


^ r 


= (-e^ -»-lc't^ - [cxo-t 


_,r» 


-h fc^'c.— ^ ‘^*®‘ 

^ L c*»‘^--'»*'^-'*tt'!r \c i'“ 1 

_ jfx.-i+i.-i'i' t?.-e-i^.-i->xi 
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The last term in this result will be neglected since it is order (c/R)^ 
smaller than the first term. Including the case r < r^^, which only intro- 
duces an effect of the droop and sweep, the result is 

^ (z,') ^ 

-4- C ('Kc— 

-t- [ (^.-e -xjiy si 

- s(c-o 


wh^e 6 (r) is the delta function (an impulse at r = 0) . Operating witli 
o '^k and integrating the second and third terms by parts gives: 




Sc 

"+■ Ss.'lx' C^Xc— Vx"S^ if 

In- • t — ■K»'i -Ht'iY 

So 

3o 't*' (Xc— VW. 

— s', -X 

- (K'^* s' 
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Finally, the torsion terms are introduced; the inertia and centrifugal forces 
directly due to the blade bending motion are extracted; and the fourth term 
above is again integrated by parts. Thus the inertia force is: 






sLr =r 




— yvi ^ 

— ^ ^ *br 

■+- L \'«. “j 4^ 

— (irpAj^ “♦• Waj”^^ Sfp^ Vv*4x- 

^ ® ^^c'-^ai^'^"3Vs.| 

SX*" ® ^<nl S**v^lr 


The structural moment (from section 2.1.5) is 


sl> 


[^ Ex^^-tf + y* 3" ^ 
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The aerodynamic moment about the tension center (x = x^) at r, due to the 
blade loading at the elastic axis at station p is 

SO 

<)' 4c - 3 > • ^ 


2.2.8 Elastia torsion equation.- The equation of motion for elastic 
torsion is obtained from 






4r- 


where M is the moment about the elastic axis at r, and 


The inertial moment is 






V V r (x.+xvt -4- 1 







Operating with 
integrations in 



( *^ -4- (%o-^ ) '‘t 4^ 

— Csi-d -h- > X. ^ 4^ 

^ I 

-— -4 - ”i^ * ^ 

— (i®‘C — l'§'^ vv*4j?^ 


ft ^1 (•••)dr and changing the order of the p 
FA K 

the second term gives 


and r 


^ 4jt r=^ \ ^v. <4>r 

’^FA ^ 

Hfa ^ 

^fA 
^*'eA 

— (lo't ~'5ILo‘fc')V ^\ 


where 


Finally, introducing the torsion terms by expanding the unit vectors, the 
inertia force is: 

O' 

^ C ® <i-r 

'Va 

— ^ ‘ (xx'^ 4-r* 


Where Ig = / (x + z )dm is the section pitch moment of Inertia, about the 
elastic axis. In the centrifugal acceleration we have neglected a number of 
terms due to the blade torsion motion which are the same order as the pro- 
peller moment, but which are normally much smaller than the structural moment. 
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2 1 

With the centrifugal tension T - n S pindp , the structural moment 


(from section 2.1.5) is: 

‘fA ^ 




fA 

~ ^ ^ ^W> ®e, 1 

- C (EXx,t - EX 

. . 


Finally, the aerodynamic moment about the elastic axis at r is 


M, = 


SO 


hlA^, 

2 ) 1 - 


— (>o't + 






'^%<. — "%%s ^ 


‘ (xo-?-+- ^o'^V 


— M. 
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Hence 




where 




2.2.9 Rigid pitch equation.- The equation of motion for rigid pitch is 
obtained from + Mp^^ = where 

“ ^PA • M C “Spa **■ *♦• M 


tnd M xs the moment about the feathering axis at r = r . The inertia 

FA 

moment is 


AA. 




So 




FAn 


= ” >^X4- — ( 4v^4j-- 

FA VPA ) 


C -X,1^ - X JL - 

^ 4c-* 
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where 


^fA ‘Fa 

Introducing the torsion terms by expanding the unit vectors, the inertia force is 

, (' XgS (t^*s-4;^*B><Lr — •■^'^•'tr 

JTf* ® If* 

- Spp^l Iff. - Was'2» ■+ 

* f^-e-x.-i-xi.-^ -u.-e-x,i')l(j*'> 

■*■ ( ifK-y'ti ~ * 

(’^W^ — >.o1^ AxT^"^3 

The aerodynamic moment about the feathering axis at r^^ is 
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The aerodynamic and inertial moments about the feathering axis are 
reacted by moments due to the deformation of the control system. The restor- 
ing moment acting on the blade, about the feathering axis, is -M . It is 
given by the product of the elastic deformation in the control systL, and the 

control system stiffness K : 

6 

where the rigid pitch p^ consists of the kinematic coupling and the blade 
commanded pitch angle: 


The first two terms are the lateral and longitudinal cyclic pitch control 
inputs; the next terms are feedback from the governor, and kinematic coupling 
due to the rotor mast bending. The term -Kp^q^ is the kinematic pitch/ 
bending coupling due to the control system and blade root geometry, where q 
is the 1-th bending degree of freedom (introduced below). Similarly, Kp^ is 
the pitch/flap coupling for the gimbal or teeter motion. For the rigid flap 
motion of the blade, this coupling is usually expressed in terms of a delta- 
three (6^) angle, such that K^, = tan 6 ^. Finally, the term is the 

pitch change due to the rotor azimuth perturbation with a fixed swashplate. 

For rigid control system (K^ very large) the rigid pitch equation of motion 

reduces to p = p . 

o 

Including control system damping in the restoring moment gives 
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where is the viscous damping coefficient. For consistency with the 

elastic torsion equations, the control system stiffness can be written in 
terms of the nonrotating natural frequency of the blade rigid pitch motion, 



and the damping coefficient in terms of a structural damping coefficient g^ : 

C-© * 

Then the structural pitch moment is 

NVpAg = c ‘Jo (fo — (p«* ^ 

* ^ A 

2.2.10 Root force. -The net force of the m-th blade acting on the hub is 
= F. ~ F-. The inertial force is. 



and the aerodynamic force is 

^ ^ C -4- ^ 

The components of the total hub force In the nonrotating frame are 

i = IP--' = H-e, 


2.2.11 Boot moment.- The net moment of the m-th blade acting on the 


rotor hub is 


M 


(wv> 
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The inertial moment is 


} o 


and the aerodynamic moment is 


7a^ = si ^ 


The components of the total hub moment in the nonrotating frame 


are 




M — £T /^^ 

V»A =. > 






Note that the (torsion) component of the root moment in the rotating 
frame is neglected compared to the (flap) and (lag) components. 

The flapwise root moment in the rotating frame gives the pitch and roll 
moments on the gimbal: 



E 

v-vs 1 


^ o 

- £ 

s 1 

My = 

Ni 

<£: 

■«?4 = 

^ AA 

& 1 

The inertial moment is 





=x s' ( 

C ) v-^ 
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and the aerodynamic moment is 


2.2.12 Girribdl equation.- The equations of motion for the gimbal degrees 
of freedom are obtained from the pitch and roll components of the total rotor 
hub force. Allowing for a gimbal spring and damper in the nonrotating frame 
reacting the rotor moments, the equations of motion are 





« 

-4- C<a, 

-t- kc. 

0 

11 


The gimbal hub spring and damper constants can be written 

-O 

^ ^ ^ 

R 2 

Where I = / r m dr and I, is a characteristic inertia of the blade, and 
o o b * 

V is the rotating natural frequency of the gimbal flap motion. To allow 

Lt 

for different longitudinal and lateral hub spring rates, v and v can be 

LxL \jo 

used for the B ^ and ^ equations. 

Ot bb 

2.2.13 Teeter equation.- The equation of motion for the teeter degree 
of freedom of a two-bladed rotor is obtained from equilibrium of flap moments 
about the teeter hinge. Allowing for a teeter spring and damper in the 
rotating frame, the equation of motion is 
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-= o 


and K^ are the damper and spring constants about the teeter hinge. In 
terms of the natural frequency and damping coefficient, we may write 

= 2 

where I « /\^mdr. 
o o 

The teetering moment is the root flapwise moment from the two 

b] ades : 

2M-r = ^ , C-i'T 

Wi jrl 

Where again ^ ^ Sl ^ ^ 

T4 

2.2.14 Modal equations.- Consider the equilibrium of the elastic, 
inertial, and centrifugal bending moments. From the results of section 2.2.7 
these terms give the following homogeneous equation for bending of the blade: 

— SZ!*’ b V — S2. 

—4— v*« c 

-5^- 


I 



This equation may be solved by the method of separation of variables. Writing 






it becomes 




This is the modal equation for coupled flap/lag bending of the rotating blade. 
It is an ordinary differential equation for the mode shape rt(r) ; this mode 
may be interpreted as the free vibration of the rotating beam at natural 
frequency v. 


This modal equation, with the appropriate boundary conditions for a 

cantilever or hinged blade, is a proper Sturm-Liouville eigenvalue problem. 

It follows that there exists a series of eigensolutions "rf (r) of this equa- 

2 ^ 

tion, with corresponding eigenvalues . The eigensolutions or modes are 
orthogonal with weighting function m; so if i k, 




These modes form a complete series, so it is possible to expand the rotor 
blade bending as a series in the modes: 

The bending modes are normalized to unit amplitude (dimensionless) at the tip: 

1^(1) I = 1- 

Consider the homogeneous equation for the elastic torsion motion of the 
nonrotating blade, i.e., the balance of structural and inertial torsion 
moments. The results of section 2.2*8 give 
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The equation for the torsion motion of a rotating blade, including centrifugal 
forces and some additional structural torsion moments could be used instead. 
For the torsional stiffness typical of rotor blades these terms have little 
effect howe er, and the nonrotating torsion modes are an accurate representa- 
tion of the blade motion. Solving this equation by separation of variables, 
we write 0^ = ^(r)e , so 

CGsT “+* ^ ^ Qi 

This equation is a proper Sturm-Liouville eigenvalue problem, from which 

it follows that there exists a series of eigensolutions C, (r), and corres- 

2 ^ 

ponding eigenvalues “ 1...“). The modes are orthogonal with weighting 

function I , so if i k 

O 

21© AX- Ci 

The modes form a complete set, so the elastic torsion of the blade may be 
expanded as a series in the modes : 

©t “ piCt") 

c ^ * 

These modes are the free vibration shape of the nonrotating blade in torsion, 
at natural frequency The torsion modes are normalized to unity at the 

tip, = 1. 

2.2.15 Modal expansion.- The bending and torsion motion of the blade is 
expanded as series in the normal modes. By this means the partial differ- 
ential equations for the motion (in r and t) are converted to ordinary 
differential equations (in time only) for the degrees of freedom. 

For the bending we write; 

1= I ' ' 
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where are the rotating, coupled flap/lag bending modes defined above. 

These modes are orthogonal and satisfy the modal equation given above. The 
variables are the degrees of freedom for the bending motion of the blade. 

For the blade elastic torsion we write 

where are the nonrotating elastic torsion modes. These inodes are 

orthogonal, and satisfy the modal equation given above. The variables 

p (i > 1) are the degrees of freedom for the elastic torsion motion of the 

^ o o c 

blade. The degree of freedom for rigid pitch motion is = 6 = (6 - 6 ) 

+ e . For rigid rotation about the feathering axis, the mode shape is 
con 

simply 5 =1. Thus the total blade pitch perturbation is expanded as the 

o 

series : 

0 - r 

L = O 

The total blade pitch 0 (mean and perturbation) is then: 

© = + © + - 4 - t f! ('•t') 1/ 


The partial differential equation for bending of the blade is obtained 
from The ordinary differential equation for the k-th bending 

mode (the equation) is then obtained by operating with ’ (•••)‘i^^ 

(which has already been done in section 2.2.7). The modal equation is used 
to introduce the bending mode natural frequency into the equation, replacing 
the structural and centrifugal stiffness terms, and the orthogonality of the 
bending modes decouples the inertial and spring terms as fo7_lows; 
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where 



C ’ 7 

No I*. 


The partial differential equation for torsion of the blade is obtained 
from 3M^/dr. The ordinary differential eqviation for the k-th torsion mode 
(the pj^ equation) is then obtained by operating with . . .)dr (which 

has already been done in section 2.2.8). The modal equation is used to 
replace the structural stiffness term with the torsion mode natural frequency, 
and the orthogonality of the modes decouples the inertial and spring terms as 
follows: 
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where 




^5> Aj- 


2.2.16 Lag damper^- Articulated rotors usually have a lag damper, which 
has an important influence on the blade loads. Therefore a lag damping term 
is added to the blade bending equation of motion as follows: 
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where and Is the lag damping coefficient (Ij^ is a char- 

acteristic Inertia of the blade, used to normalize the inertial constants as 
described in the next section). The quantity . TTj^'(e) is the slope of 
the k-th bending mode in the lagwise direction, just outboard of the lag hinge. 
The manner in which the lag damping enters the equation of motion is obtained 
by a Galerkin or Rayleigh-Ritz analysis. The lag damper results in a bending 
moment at the lag hinge. Thus it is necessary to evaluate moments at the 
blade root by integrating along the span, which has in fact been our practice. 

Note that structural damping has also been included in the bending 
equation, modelled as equivalent viscous damping. The structural damping 

6g (nqual to twice the equivalent damping ratio) in general is 
different for each degree of freedom. Structural damping is included in the 
torsion equations in a similar manner. 


Consider also a nonlinear lag damper, for which the lag moment opposing 
the motion is proportional to I ^ at low lag velocity (hydraulic damping) 
and constant at for lag velocity above (friction damping): 


- 




where 

5 


i5, 



Hence the term 


is added to the right-hand side of the bending equation. Here linear damping 

Is included on the left-hand side still, but only to Improve the convergence 

of the solution; so the C t term must be subtracted from M 

^ lag. 
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2,2,17 GvccoitationaZ forces^^ The acceleration due to gravity is 
g = gic^ = S^5p^£» where g is the gravitational constant, is the vertical 

vector, and R is the coordinate transformation matrix between the rotor 

br 

shaft axes (S frame) and the aircraft body axes (F frame, see section 4,1.2). 
In terms of the aircraft trim pitch and roll angles, the vertical vector is 





Tt 






-+* 




(see section 4.1), The gravitational forces acting on the rotor blades may 

be accounted for by substituting *a for a^ , the hub linear acceleration, 

o o 

— % 

Thus the components of g in the S frame are subtracted from the components 
of the hub acceleration in the nonrotating shaft axes: 

-*■ -V 


2.2.18 Equations of motion.- The rotor blade equations of motion are 
now obtained by substituting for the expansion of the bending and torsion 
motion as series in the modes of free vibration. Names are given to all the 
inertial constants. Also, the equations of motion, hub forces and moments, 
and inertia constants are normalized at this point in the analysis, using the 

characteristic blade inertia I, , and the blade Lock number y = pacR^/L is 

^ R ^ 2 

introduced. (A good choice for this characteristic inertia is I, = / r mdr.) 

b o 

The inertia constants are divided by with this normalization denoted by 

a superscript The blade equations of motion are divided by . The 

hub forces and moments are divided by , so they appear in rotor coeffici- 

ent form. The equations of motion for blade coupled flap/lag bending and for 
blade rigid pitch/elastic torsion are thus: 
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The inertia constants are defined in section 2. 2.19* 

In rotor coefficient form, the rotor hub force and moment are 
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or 
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So the blade root force and moment are resolved in the nonrotating frame, and 

1 ^ 

-her: filtered by the hub operator tt £ . The harmonics of the forces in 

^ in=l 

the nonrotating frame can be related directly to the harmonics of the rotating 
forces; the solution of the support equations of motion requires however the 
hub forces and moment in the time domain. 


The components of the blade root force and moment in the rotating frame 
are as follows: 


^ fvx = ( EiA 




ka** / V. *» X 
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The inertia constants are defined in section 2*2. 19, 


Note that the total hub forces due to the rotor linear acceleration are 
simply 



where is the normalized mass of a blade. Since the rotor mass is included 

in the aircraft mass, these hub linear acceleration terms should be omitted 
when and are evaluated for the aircraft equations of motion. 

These terms should be retained however when evaluating the actual blade root 
forces. 

Similarly, since the rotor weight is included in the aircraft weight, 
the corresponding gravitational force terms are omitted. 

Dividing by (1/2)NI^, the gimbal equations of motion are 

V -+ ccT ^ = o 
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where 


O-'Q 


Z 
— \ 

2. 


o-a 



to 

z 

v-v •i. \ 
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Dividing by 

— '5 


NI^, the teetering equation of motion is 



where 



- 1 _ 


"2, 


C 'I - 


^ o 


The equation of motion for the rotor speed perturbation i is obtained 

s 

from equilibrium of the rotor torque. The speed perturbations of the two 
rotors are coupled by the helicopter transmission, so the equation of motion 
for is best derived with the body equations. 

Finally, the aerodynamic forces required for the blade equations of 
motion and the rotor hub reactions are as follows: 
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2.2.19 Imrtial constants.- The normalized Inertial constants required 
for the blade equations of motions and the hub reactions given above are 
defined as follows: 
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We have used the relation 


o 


for elastic torsion (k ^ 1) , and 

>(„ = — u.-t — n;:i 

■*■ 'ct* (r-r^^'-i 

— -•- Cii»^-SfAn'5Lfc ••-f4fAJ-S^'Ae'^^6')(»■-^f(^ 

for rigid pitch (k = 0); or for k > 0 
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where 
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Also then. 






A*. 







for the aerodynamic coefficients. 

The blade inertial and structural properties (m, , x , El, I , GJ, etc.) 

X C 0 

will be defined at a series of radial stations, r^, with linear variation 
between. 


The blade bending and torsion mode shapes will be evaluated at M + 1 
equidistant radial stations: r = 0, Ar, ...MAr where Ar ■= 1/M. The 

inertial coefficients are then calculated by numerical integration (using the 
trapezoidal rule) over these radial stations. 


A concentrated mass at the blade tip (r = 1) will be allowed, with a 
corresponding center of gravity offset. This tip mass contributes an' equiva- 
lent distributed mass, as follows: 


= ;:4'co ( ') 



(.I'i 
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where Ar is the segment length for numeral integration. Alternatively, 
the tip mass can be included in the distributed mass directly (to avoid 
difficulties with the (x^ - x^)m terms evaluated at r = 1) . 

The total mass of the blade can be specified, so 


Alternatively, a point mass can be added at r = 0 to account for the weight 
of the hub. 


2.2.20 Aerodynamic spring and damping,- To improve the convergence of 
the solution for the blade motion, spring and damping forces should be included 
on the left-hand- side of the equations of motion. The required perturbation 
aerodynamic forces are: 






f ft 4- f ^ 

4- ^ -f- ^ Mf, ». 


These terms will be added to both sides of the equations of motion, so they 
need not be exact values of the damping and spring forces, but only close 
enough to achieve good convergence (see section 5.1). The damping terms are 
needed to avoid unrealistic resonant amplification of the harmonics near the 
natural frequency, and the spring terms help obtain the correct phase of the 
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response quickly. Following the aeroelastic analysis (section 6.1.4), the 
following expressions are used for the aerodynamic coefficients. 



- 80 - 




where 
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2.3 Blade Bending and Torsion Modes 

2.3.1 Coupled bending modes of a rotating blade,- Equilibrium of the 
elastic, inertial, and centrifugal bending moments on the blade gives the 
differential equation for the coupled flap/lag bending of the rotating blade 
(see section 2.2.14). For free vibration— the homogeneous equation with 
harmonic motion at the natural frequency v — we obtain the modal equation 
for bending of the blade: 
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x^k is the bending deflection (mode shape) , 


Here TfCr) = z - 
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is the bending stiffness dyadic SI - Qk^ is the rotor rotational speed, and 

V is the natural frequency of the mode. The boundary conditions are as 
follows : 

(a) at the tip (r = R) : Elq" = (Elq")’ = 0 

(b) and at the root (r = e) : n = n' = 0 for a cantilever blade; 

n = 0 and Elq" = ' for an articulated blade. 

The root boundary condition is applied at the offset r = e to allow for 
hinge offset of an articulated rotor, or a very stiff hub of a hingeless rotor. 
Different offsets can be used for the out-of-plane and Inplane motion 

(e^ and e^) . With the hinge springs at an angle 0^ from the hub plane, the 
hinge spring dyadic is 

-f- ^ -4- 

-+ (^L 

where Kp is the flap spring and is the lag spring constant. 

It is useful to be able to use for the pitch angle of the structural 
principal axes the effective angle 6^ “ (R6 . The parameter (R is zero for 
no structural coupling of the inplane and out-of-plane blade motion; and 

-fi2- 
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;5P = 1 for complete coupling. For the hinge spring pitch angle 0^, an input 
value can be used; or can be used; or more generally 6g + ®so* 

This differential equation is an eigenvalue problem for the mode shapes 
n and the natural frequencies v. The equation and boundary conditions 
constitute a proper Sturm-Liouville problem. It follows that a series of 

-A, 

eigensolutions or modes n^(r) exists with corresponding natural frequencies 
v^; and that the modes are orthogonal with weight m. Hence if i k. 



The frequencies satisfy the energy balance relation: 

V — 




St f 


Xr 


The modal equation will be solved by a modified Galerkin method follow- 
ing reference 3. This approach works better for large radial variations in 
the bending stiffness than does the Rayleigh-Ritz method in standard form. 
Write the differential equation as 

M** — ( Sr ^ - o 


1 


r- 


Ei: 



M 


o 


with boundary conditions M = M’ = 0 at r = 1, and n = 0 and M = 
r = e. The deflection and moment are expanded as finite series in the func- 


tions f^ and g^ 
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It is required that each of the functions and ^ satisfy the boundary 
conditions; then the sum automatically does. Since a finite series is 
required for numerical calculations, this will be an approximate solution. 
For best numerical accuracy the functions f^ and must be chosen so that 
the lower frequency modes can be well represented by the truncated series. 
Substituting these series into the differential equations and operating with 


gives 

o 

^■ci - c Ai i<- — o 


Integrating by parts and applying the boundary conditions gives 

^ = T»-A' 1'^ - t + C t' M 
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so the- first equation becomes 

I -i 


/ 

Tc.; 


^ SjlC\ix 
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Hence the problem reduces to a set of algebraic equations for c = Lc^] and 
■d = [d.] 

c d + l>c — =-=> 

c£T c — A> 

or 

Cc^K'd^ -e b = c> 


For simplicity the functions used for the moment expansion are 
Then the coefficients of the matrices are 


K, = V 

r' _^U -i M 
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(Note that using = El f^" would give C = . A so d*- and this 

solution would reduce to the standard Galerkin form.) 

^ The eigenvalues of the matrix b'^CA-V + D) are the natural frequencies 
V of the coupled bending vibration of the blade; and the corresponding 
eigenvectors c^ give the mode shape •?. As a final step, the modes are 
normalized to unity at the tip: (-S(l)| . 1. This modified Galerkin approach 

equivalently replaces the Rayleigh energy expression for the natural frequency 

(given above) by 

d. ex';;« 


The blade nonrotating modes and frequencies can be obtained using 


't, 

_v 






Ss (c^ 


A convenient set of functions for f^ are the bending mode shapes of 
a nonrotating, uniform beam. Such functions will satisfy the required boundary 
conditions, and furthermore are orthogonal (necessary for good numerical con- 
dltioning of the Galerkin solution) . Let be the series of eigensolutions 
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4 4 4 

of the differential equation d w/dx = a w with appropriate boundary condi- 
tions, Using these functions for both out-of-plane and inplane deflections 
gives 


V, 





=: 


j=. 




where x = (r - e)/(l - e). 

For a uniform hinged blade, the nonrotating mode shapes are: 

2, ^ 


w ^ 


where a is the solution of 

Loo^a O.. = 



The first mode (a = 0) is w = x. For a uniform cantilever blade, the non- 
rotating mode shapes are 

^ ^ C ^ C. 

2- ^ Srs^^ S i/k A 

^ o«a Ok -4- 

2 ^ &AaaQa. Oi 
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where a Is the solution of 


Ok 



The values of a for the lowest modes are given in the table below. 


Mode Hinged 


Cantilever 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


0 

3.926602313 

7.068582747 

10.21017612 

13.35176878 

16.49336143 

19.63495409 

22.77654674 

25.91813940 

29.05973205 


1.875104069 

4.694091134 

7.854757439 

10.99554074 

14.13716839 

17.27875953 

20.42035225 

23.56194490 

26.70353756 

29,84513021 


The centrifugal force is required for the bending mode calculation. 

With the section mass defined at radial stations (i = 1 to M) the centrifu- 
gal force is 




M 


£1 



where r^_j^ 1 r < r^. Then for linear variation of the section mass between 
the stations r^, the integrals can be evaluated as follows: 
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where 
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2.3.2 Articulated Blade Modes.- For an articulated blade the modal 
differential equation need not be solved if the higher bending modes are not 
required. Rigid flap and lag motion about the hinges gives the two lowest 
frequency modes : 



Note that separate hinge offsets may be used for flap and lag motion. The 
natural frequencies are obtained directly from the energy relation, as 
follows ; 
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2.3.3 Torsion modes of a nonrotating blade.- Equilibrium of the elastic 
and inertial torsion moments (see section 2.2.14) gives the model equation 


with the boundary conditions ^* = 0 at the tip (r •= R) and ^ = 0 at the 

root (r = rp^). The modes are orthogonal with weight I_, so if i ft k 

6 

si 


The frequencies satisfy the relation 
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These are the nonrotating torsion modes, so the solution is independent of 
the rotor speed or collective pitch. 


The equation is solved by the modified Galerkin method, as described in 
detail above for the bending modes. Write the differential equation as 

i _ 


“0 


U> 


M 


S' - 


Expand the torsion deflection and torsion moment as series: 

"T — ■£ . 'ji ( 


where the functions f^ and satisfy the boundary conditions on J and T. 
Substitute these series in the equations, operate with 


integrate by parts and use the boundary conditions 




to obtain 



Hence the problem reduces to a set of algebraic equations for ^ « 



c'^c — A A = '=> 


or 


A"' c"^ 



For simplicity, the functions used for the torsion moment at g = fl. 
the coefficients of the matrices are 

'■«> = 'SrV s' 

~ 

(Note using g^ = GJf^' would give the standard Galerkin result.) 
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The eigenvalues of the matrix B (C A C ) give the natural frequencies 
of the torsion vibration, and the corresponding eigenvectors for give the 
modes. Finally, the torsion modes are normalized to unity at the tip, 

f(i) = 1. 

A convenient set of functions to use for is the solution for the 

torsion modes of a uniform beam; 

== \ 3 

These functions satisfy the boundary conditions, and will often be close to 
the true mode shapes. 

2.3.4 Kinematic pitoh/bending coupling.- The kinematic pitch/bending 
coupling and the pitch/gimbal coupling have a significant role in 

the rotor dynamic behavior. The definition of Kp^ is the rigid pitch motion 
due to a unit deflection of the i-th bending mode: Kp^ = -d6/dq^. For an 

articulated rotor, the first "bending" modes are rigid lag and flap motion 
about the hinges. The pitch/flap coupling is often defined in terms of the 
delta-three angle i K = tan It is possible to simply input these 

kinematic coupling parameters to the dynamics analysis, if values are available 
from either measurements or some other analysis. It is also desirable to be 
able to calculate the coupling from a model of the blade root geometry. 

Figure 10 is a schematic of the blade root and control system geometry 

considered, showing the position of the feather bearing, pitch horn, and 

pitch link for no bending deflection of the blade. The radial locations of 

the feather bearing and pitch link are r and r respectively; the length 

ro r n 

of the pitch horn is orientation of the pitch horn and pitch link 

are given by the angles (}) + 6-_ and • Control input produces a vertical 

Jr H / j t ti 

motion of the bottom of the pitch link, and hence a feathering motion of the 
blade about the pitch axis. Bending motion of the blade, with either struc- 
tural flexibility or an actual hinge inboard of the pitch bearing, produces 
an inplane or out-of-plane deflection of the pitch bearing. With the bottom 
of the pitch link fixed in space, a pitch change of the blade results. 
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CONTROL 

INPUT 


Figiire 10. Schematic of blade root and control system 
geometry for calculating the kinematic 
pitch/bending coupling. 
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The vertical and inplane displacements of the pitch horn (the end at 
due to bending of the blade in the i-th mode are: 


~ (^i (Q&) 


The kinematic pitch/bending coupling is derived from the geometric constraint 
that the lengths of the pitch horn and pitch link are fixed. The result is: 

VCf. = 


Similarly, for a gimballed (or teetering) rotor the pitch/flap coupling is: 



-4- 4^U ^ 


_ A%Vc.JU»' t 






2.3.5 Blade pitch definition,- Outboard of the trim pitch angle 

is given by the collective and twist angles, while inboard of r^^ it is 
given by just the twist angle: 
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(see section 2.2.2). It is convenient to use the collective pitch value at 
75% radius, 0^^. Then 


© = 



-H 






which requires 6^ (r = .75) = 0 (but no change to 0 for r < r ) For 

LW p^/ . 

a rotor without a pitch bearing, it is more appropriate to maintain continuity 
of 0 by adding a linear term inboard of r^^; 


0 = ^ 



^ d 

1 ‘Va®’' 

H- 



For the structural and inertial analysis the pitch angle is multiplied by the 
structural coupling parameter (R. 

The twist distribution is required at the radial stations for 

which the inertial and structural properties are defined; and at the radial 
stations at which the aerodynamic forces are calculated. The aerodynamic 
twist definition can include the zero lift axis pitch 0^^^ (see section 2.4.1). 
Frequently, a linear twist distribution is used, for which 

-Vw.» ^ 
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2.4 Aerodynamic Analysis 


In this section the aerodynamic forces and moments on the rotor blade 
are derived. The general case of a rotor in high or low inflow, axial or 
nonaxial flight is considered, including the effects of reverse flow and 
large angles. Lifting line theory (i.e., strip theory or blade element theory) 
is used to calculate the section loading from the airfoil two-dimensional 
aerodynamic characteristics, with corrections for yawed and three-dimensional 
flow effects are required. The unsteady aerodynamic lift and moment are 
obtained from thin airfoil theory, and a dynamic stall model accounts for the 
unsteady aerodynamic phenomena at large angles of attack. 

2.4.1 Section aerodynamic forces.- A hub plane reference frame is used 

for the aerodynamic forces. All forces and velocities are resolved in the hub 

plane (i.e., the B coordinate system). The hub plane reference frame is fixed 

with respect to the shaft, hence it is tilted and displaced by the shaft 

motion. Figure 11 illustrates the forces and velocities of the blade section 

aerodynamics. The blade pitch angle is 0, measured from the reference plane. 

The velocity of the air as seen by the moving blade has components u^, u , 

2 2 1/2 ^ 

and u^, resolved with respect to the reference frame; U = (u + u ) is 

i p ^2^ 

the resultant air velocity in the plane of the section; and ^ = tan 
is the induced angle. The section angle of attack is 

o/ ^ © -4- ^ 

where e is the pitch of the aerodynamic zero-lift axis of the section 

relative to the structural/inertial principal axis at pitch angle 6 (6^^ may 
vary along the span, and should not therefore be included in the definition 
of the section aerodynamic coefficients as a function of a; 6^^ can however 
be included in the aerodynamic twist distribution, if 0^^ is defined 
separately for the Inertial /structural pitch and for the aerodynamic pitch) . 
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LOOKING SECTION VELOCITIES 

INBOARD 


Figiire 11. Rotor blade section aerodynamics. 
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The velocity is in the hub plane, positive in the blade drag direction; 

is in the hub plane, positive radially outward along the blade; and u^ 
is normal to the hub plane, positive down through the rotor disk. The aero- 
dynamic forces and moment on the section, at the elastic axis, are defined 
as follows: L and D are the aerodynamic lift and drag forces on the section, 

respectively normal and parallel to the resultant velocity U; and are 
the components of the total aerodynamic force on the section resolved with 
respect to the hub plane, normal to and in the plane of the rotor; F^ is the 
radial drag force on the blade, positive outward (the same direction as 
positive Uj^) ; and is the section aerodynamic moment about the elastic 

axis, positive nose up. The radial forces due to the tilt of and F^ are 
considered separately, hence F^ consists only of the radial drag forces. 

The section lift and drag are 

e Cjj -f. 

b — Z ^ 


where U is the resultant velocity at the section, p is the air density, 
and c is the chord of the blade. (The air density can be dropped since all 
quantities are actually dimensionless, based on p, and R.) The section 
lift and drag coefficients, c^ and c^, are functions of the section angle of 
attack and Mach number: 

sJ- — s -+ = < 3 - 4 - 

M U 


where the tip Mach number (the rotor tip speed fiR divided by the 

speed of sound). is the unsteady aerodynamic lift force. The radial 

drag force is 
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This radial drag force is based on the assiiraption that the viscous drag force 
on the section has the same sweep angle as the local section velocity. The 
moment about the elastic axis is 


M.. = — * N'l.i 

•= \ ccji -V -t- 


where is the distance the. aerodynamic center is behind the elastic axis, 

is the section moment about the aerodynamic center (positive nose up) , 
and is the unsteady aerodynamic moment. 

The components of the section aerodynamic forces relative to the hub 
plane axes are then 

L ^ ^ L- U-r ^ 

= L ~ C /iX 


Substituting for L and D, and dividing by a, the two-dimensional lift-curve 

slope, and by c , the mean section chord (which enter the Lock number Y 

m ^ 

also), we obtain: 
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The net rotor forces required are obtained by integration of these section 
forces over the span of the blade: 

= V. ^ - li 

lx- - ^ 
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where 


*V “ f ^ ^ + W*. - Wai^ — (— +- 

C f,"!* 


To numerically integrate the aerodynamic loads over the blade span, 
define K radial segments by the boundaries 




where r = 1. For the k-th segment, the airloads are calculated at the 


center: 


^ = i (r^ -H 




Then the spanwise integration is approximated by a summation over all segments; 




In summary, the rotor blade aerodynamic forces are evaluated as follows. 
First the section velocity components and pitch angle are evaluated, and then 
the angle of attack and Mach number. Next the section aerodynamic coefficients 


-lOP- 



are obtained (see section 2. A. 4), and from them the section force components 
and moment. Finally, the section forces are integrated over the rotor radius 
to obtain the required generalized forces. 

2.4.2 Blade velocity.- The air velocity seen by the blade section is 
due to the rotor rotation, the helicopter forward speed, the rotor and shaft 
motion, and the wake induced velocity. The rotor is rotating at speed 
The velocity of the air as seen by the rotor disk has the following dimension- 
less components in the shaft axis system: y , positive aft; y , positive from 

X y 

the right; and y^, positive down through the disk: 

^ 

Often the lateral velocity component y^ = is assumed to be 

zero in the rotor aerodynamic analysis, and indeed it is small for most flight 
conditions. An exception is the case of sideward flight. An alternative to 
including y is to rotate the shaft axes until ^ • “v? = 0, but that would 

y ® 

imply a redefinition of the rotor zero azimuth position for every flight state. 
Such a redefinition of \jj is not desirable since it changes the values of 
parameters such as the control system phasing, and even changes the definition 
of the harmonics of the rotor motion. Hence it is preferable to directly 
include the effects of the lateral velocity in the analysis. 

The rotor wake-induced velocity is X. = v./S^R, normal to the rotor disk 

11 

and positive downward. A simple model may be used, such as a uniform or 
linear variation over the disk, or calculated nonuniform induced velocities 
may be used. For the latter case, all three components of the wake induced 
velocity (in shaft axes) will be considered: 

The rotor advance ratio y and Inflow ratio X as conventionally 
defined are here 
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Vi, 


These are the dimensionless inplane and normal components of the total velocity 
seen by the rotor disk. The hub plane angle of attack and yaw angle are then 

= TkC' 




N/ 


V = 




-N 


At 




N>; 

-snc-ft. 


4— 


Here V is the helicopter velocity, with angle of attack relative to 

HP 

the hub plane is positive for forward tilt of the rotor disk) . The 

advance ratio p is zero for hover and axial flow, and y > 0 for helicopter 
forward flight. 

The aerodynamic gust velocity has components u^ • v„, and w„ in the 
shaft axis system, normalized by dividing by the tip speed fiR. The longi- 
tudinal gust Vg is positive from the front, the lateral gust v^ is posi- 
tive from the right, and the vertical gust w^ is positive upward 

(^gust " " ^G^s '^G^ relative to the rotor). This gust velocity is 

evaluated at azimuth angle and radial station r on the rotor disk 
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(at r = r(cos + sin relative to the hub). The quasisteady 

shaft motion and the gust velocity at the rotor hub will be included in the 
advance ratio component and (see section 4.1.2). 

The blade and shaft motion have been defined in the inertial analysis 
(section 2.2). The resulting velocity components in the rotor shaft axes are 
thus : 

lA-r = r 4- 4“ 

— X ^ 

-4- (/* i cv'y V" -4.V^C».'P 

— )— C/*' p. ^ ^ 

I'i.o-i— v,»'^ ^ — r^. 

— (\£4/ij')Q (5-.*^— xltV ^ 

-4* >U'+ 

^/^vS.v^H' Co(^ 
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+ ( ?a. -'««. -A.»S -/.^o(,'^ H- r C K-s- - 

and the pitch angle is 

oo 

IS o 


In body axes, the trim velocity vector is fixed with the reference frame, and 
would therefore tilt with it. With inertial axes however, a tilt of the rotor 
by the shaft motion results in a small change in the directions of the com- 
ponents of y as seen in the reference frame. All the ya terms in the 
expressions above for u^, Up, and Up result from such tilt of the inertial 
■xes relative to the trim velocity vector. The aircraft body yaw, pitch, and 
roll will be defined as body axis motion however. Hence the body Euler angles 

are not to be included in the evaluation of a , a , and a for the blade 

X y’ z 

velocities. 

2.4.3 Induced velocity.- For the case of uniform inflow, the rotor wake- 
induced velocity is obtained from the momentum theory result 


= T7 


vV.cS. ♦ 


2 2 2 

where X = y + X and y = y + y . Empirical correction factors k, and 
z i X y h 

•Cj are included for the effects of nonuniform inflow, tip losses, swirl, 
blockage, etc., in hover and forward flight. An iterative solution of this 
equation for X^ is necessary; 
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C- + 


• -^-/* •'H "I 

-1 




3/z. 


with 



CT-r 

2. J /zv^ -*-yA*>/vc^‘*- 


to Start the solution; 3 or 4 iterations are usually sufficient. For the 
vortex ring and turbulent wake states this momentum theory result is not 
applicable. Thus if 

-f. ’S' k^Y <■ 


I 'll I 


where 

= 'I ^ 

The wake-induced velocity is reduced when the rotor disk is in the 
proximity of the ground plane. The effect of the ground will be accounted 
for using the following approximate expression from reference 4 for the ratio 
of the induced velocities in and out of ground effect: 


the following expression is used instead: 






1^’ 
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where z is the height of the rotor hub above ground level, normalized by 

the rotor radius; and e is the angle between the ground and the rotor wake 

(e = 0 for hover and s approaches 90® in forward flight), which accounts 

for the effect of forward speed. Note that ground effect is essentially 

negligible for altitudes greater than the rotor diameter (z > 2) or at forward 

1/2 

speeds ji > 2(C^/2) . This expression compares well with test results, 

down to an altitude of about one-half rotor radius (see reference 4). The 
rotor wake-induced velocity in ground effect is thus 






- O- 


0<»E 


^ 2 t height of the helicopter center of gravity above ground 

level; and (x^, y^, be the components of the rotor hub position relative 
to the center of gravity, in a body axis system (the F frame, see section 
4.1,5). Then the altitude of the rotor hub above ground level is 

Z = '^AC.L. — 

The vertical (T^) is defined relative to the body axes by the trim pitch and 
roll Euler angles (6^^ and see section 4.1.1). The angle between the 

rotor wake and the vertical is 

= C A - V -f, ^ 




I 



where R „ is the transformation matrix between the shaft and body axis 
S F 

coordinate frames. 

As a first approximation to the rotor nonuniform induced velocity dis- 
tribution, a linear variation over the disk is considered: 




where X, is the mean value of the induced velocity, calculated as described 
1 

above. Typically k is positive, roughly 1 at high speed; and k is 

X y 

smaller in magnitude and negative. Both k and k must be zero in hover. 

X y 

Based on references 5 to 7 we will use 


R 


X 




S 





-+ ui 


with typically f =1.5 and f =1.0. There will also be an inflow variation 
X y 

due to any net aerodynamic moment on the rotor disk. The differential form 
of momentum theory gives 

— ZC/wK^r'os«»<V’ -h 2 -Csax ^ 

including an empirical factor f^. 

With twin-rotor aircraft it is also necessary to account for the rotor- 
rotor aerodynamic interference in the wake-induced inflow velocities. The 
induced velocity at each rotor will be expressed as a linear combination of 
the isolated rotor induced velocity. Let X^j and ^^ 2 ^^ trim induced 

velocity of the two isolated rotors, calculated as above. Then the trim 
inflow ratios are 
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\, = 











Vi, 

CseR.'^^ 


Here and are the rotor-rotor aerodynamic interference factors. 

Separate values are used for the interference factors in hover and forward 
flight, with a linear variation from p «= 0.05 to 0.10. 


In summary, the isolated rotor mean Induced velocity is calculated from 
the advance ratio and thrust, 


where = 1 out of ground effect. Including the rotor-rotor interference 

and the linearly varying induced velocity components, the inflow ratios are 
then 
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for rotor //I and rotor //2. 

2.4.4 Section aerodynamic characteristics.- The section aerodynamic 
characteristics required are the static lift, drag, and moment coefficients as 
a function of angle of attack and Mach number: c (a, M) , c (a, M) , and 

Xf Q 

c^(a, M) • Most often rotor loads analyses use two-dimensional airfoil test 
data in tabular form. The aerod 3 mamic description of the blade also requires 
Bzl> zero lift angle of the section relative to the structural/inertial 
principal axis at pitch angle 0; and x , the distance the aerodynamic center 

A 

(in normal flow) is behind the elastic axis. The strict definition of 6 

ZL 

is actually the pitch of the axis corresponding to a = 0 in the airfoil data 

used. Similarly the strict definition of x is simply the location of the 

A 

axis about which the moment data c are given. It is convenient to use the 

m 

zero lift axis and the aerodynamic center, but the most important considera- 
tion is that the definitions of 0^^ and x^ be consistent with the zero angle 
of attack and moment axis definitions in the airfoil data used. 

The angle of attack a is defined in the range -180 to 180, with the 
same sense as 0. The lift, drag, and moment as a function of angle of attack 
are defined as in two-dimensional airfoil tests, where a is varied by pitch- 
ing the airfoil; and the lift is always positive vertically; the drag is posi- 
tive in the direction of the free stream; and the moment is positive nose up. 
For the rotor blade in reverse flow then (u^ < 0) , a positive pitch G or 
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positive (down) normal velocity Up gives an angle of attack near -180; which 

gives positive and c^; which are down lift L and forward acting drag D* 

The section moment is given about a fixed axis of the section* In reverse 

flow the aerodynamic center shifts to near the three-quarter chord (from near 

the quarter chord in normal flow) so it is expected that the c data will 

1 ^ 

c 1 

show a nose up moment contribution of AM = -rL or Ac = -; 7 C„ in reverse flow 

L m z 31 

(see fig. 12). 

The steady, two-dimensional airfoil data (c^, c,, and c as a function of 

Z a m 

a, M, and r) will be used in tables of the following form. The data will be 
defined at a finite set of angle of attack points. To facilitate interpola- 
tion, these points will consist of several groups, with the same angle of 
attack increment within each group. Then the set of angle of attack points 
can be specified by the a at the boundaries between the groups, and the 
indices of these points: N^, to and to (for N^“l groups). 

The organization will be similar for the variation with Mach number. For the 
radial variation, the blade will be divided into segments with the same 
section, defined. by the r at the boundaries: N^, and r^^ to 

segments. Hence the data set for the lift coefficient has the form 

N^, (k = 1 to Ng^) , o<. (k » 1 to \) 

Dk (k= 1 to N„), (k=l to N„) 

Nj., (k= 1 to Np+1) 

c^(i) for i = ^ ■+* Ja 

(((Ja- 1 ^ ^ 

and similarly for the drag and moment coefficient data sets. 


The data will be linearly interpolated over angle of attack and Mach 

number. The boundary point definitions determine the values of a and M 

for all points in the data set. Consider the angle of attack variation. The 

boundary point definition of for i « 1 to implies that the angles 

of attack for points between the boundaries a , and are 

ni i+1 


cy 


n; 





« —Mi. 
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a, deg 


Figure 12. Sketch of section aerodynamic characteristics. 
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for j “ 0 to n^) . 

Hence given 

a, we search for 




It follows immediately that 
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e>< ^ 


a 


where 
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i _ -- 
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n; 


([a] means the greatest integer In a; l.e. , Integer arithmetic). With a 
function c defined at and at linear interpolation then give 




If a < a , set j = 1 and a = 0; if a > Onvi i set j = -1 and 

X a 3i “ 

= 1. Similarly, for a given Mach number M search for k such that 
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then calculate 
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If M < , set j =1 and m. ^ = 0; if M > 


m 


j+1 


*"N„ 


set 




-1 and 


m 


3+1 
such that 


= 1. The appropriate radial station is determined by searching for j 


XT < r- 

i'' 




r. 


-M 


The aerodynamic coefficient is evaluated at the four corners, and then the 
interpolated value is 


d 
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2.4.5 Tip flow oorreations.- Three dimensional flow effects at the blade 
tips significantly alter the wing loading. Principally it is necessary to 
correct the blade element theory section loading calculation for the lift 
reduction and compressibility relief near the tip. The standard tip loss 
correction assumes that the blade has drag but no lift outboard of radial 
station r = B. Hence for a radial segment extending from r^ to r , the 
lift coefficient is multiplied by the factor 


W\OoV 







The moment and drag coefficients are not altered. 


B = .97 can be used, or 

fe - i 


N 


For the tip loss factor 


Alternatively, the tip losses can be accounted for by multiplying the blade 
element theory lift by the Prandtl function: 




_l 




An effective tip loss factor can be evaluated from this function: 

S — i — aju^ 


hence 


V — nY 
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The three dimensional flow at the blade tip increases the critical Mach 
number of the tip sections, compared to the two dimensional flow character- 
istics. This compressible tip relief may be accounted for by reducing the 
effective section Mach number by the factor 






'T\~ 


The factor f must be specified at each blade station, for the lift, drag, 

M 

and moment. 

Swept and tapered tip planforms are defined in the present analysis by 
the blade chord, aerodynamic center, pitch angle and zero lift angle, and 
center of gravity distributions (c, x^, ® and sweep of the 

blade elastic axis at the tip is neglected however. The tip planform should 
also be considered in choosing the tip loss factor and compressible tip relief 
factors for the rotor blade. 

2. A. 6 Yaued flow aorreation.- Yawed flow over the blade section may be 
accounted for using the equivalence assumption for swept wings: that the 

yawed section drag coefficient is given by two-dimensional airfoil character- 
istics, and the normal section lift coefficient is not influenced by yawed 
flow below stall. Since the wing viewed in a frame moving spanwise at a 
velocity V sin A (where V is the wing velocity, yawed at angle A) is 
equivalent to an unyawed wing with free stream velocity V cos A, except for 
changes in the boundary layer, there should then be no effect of spanwise 
flow on the loading below stall. Accounting for the effective dynamic pres- 
sure and angle of attack of the yawed section relative to the normal section 
leads to 





Cc4 2^^ CXW -A- 


(cJ 
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for the section aerodynamic coefficients in terms of two-dimensional airfoil 
characteristics. These results are largely verified by the experimental data 
for yawed wings. The section yaw angle is given by 





In reverse flow ( |a\ > 90°) the angle of attack correction is 


= 




for the drag, and 


— L -I- \‘3o* CK 

for the lift and moment. 

2.4.7 Dynamic stall model.- Dynamic stall is characterized by a delay 
in the occurrence of separated flow due to the blade motion, and high transient 
loads induced by a vortex shed from the leading edge when stall does occur. 
These features are modelled by the following procedure adapted from reference 8. 
McCi-nskey (reference 9), and Beddoes (reference 10) have found that the dynamic 
delay correlates fairly well in terms of the normalized time constant 
T = At V/c. Their results for lift and moment stall are 

^ 1 , 

McCroskey 4-3 o.CT o-Z 

Beddoes 5* H ± O lo 2.45 ^ o.S- 
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or approximately = 4.8 and “ 2.7 (a constant is also required 

for the drag stall delay). Hence the section lift will be evaluated at the 
delayed angle of attack 


where Ai|j = fiAt = t c/ |ii j (radians). A maximum value of the angle Increment 
Li L Li T 

(a - a^) should be specified In order to avoid difficulties at small values 
of u^. Thus 




— 


- ( tl- I 5 ^ 








vvcx 


is used. The lift coefficient below stall should not be affected by the 
dynamic stall model, rather the stall delay should extend the linear range 
above the static stall angle of attack. Hence the corrected lift coefficient 
takes the form 


Including the yawed flow correction this becomes 


= 










Here Ac^^ Is the lift Increment due to the loading edge vortex used at 
dynamic stall, which is discussed below. Similarly a delayed angle of attack 
Is calculated for the drag and moment from appropriate time constants and 

T , and the corrected section aerodynamic coefficients are 
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including the yawed flow correction. In reverse flow (90 - | a | < 0) the 
lift coefficient correction should be 






£ ->to + 1^0* S«jpo(^ 










4- 


and for the drag and moinent coefficients 

C ^ ^*2>* ')cjaa^I-A. -4- l%c? "#" 


When the blade section angle of attack reaches the dynamic stall angle 
“ds’ ® leading edge vortex is shed. As this vortex passes aft over the air- 
foil upper surface it induces large transient loads. The experimental data 
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of reference 11 show that the peak incremental aerodynamic coefficients depend 

* 

on the pitch rate at the instant of stall, a c/V, approximately as follows: 


Ac^ 


— 











/V 

> .OH 
<.oa 

• oz < o^c/V'^ * 

^c/V> ’OS’ 






<3^c./V < -OS 
oicA >-o5 


with ~ ^^^ds ” present model of the dynamic stall 

loads it is assumed that the incremental coefficients due to the shed vortex 

(Ac., Ac , and Ac.) rise linearly to the above peak values in the small 
£ m a 

azimuth increment (typically 10® to 15®), and then fall linearly to 

zero in the time again. Hence the model involves impulsive lift and 

nose down moment increases when dynamic stall occurs, which produce the blade 

motion and loads characteristic of rotor stall. After these transient loads 

decay the blade section is assumed to be in deep stall, and dynamic stall is 

not allowed to occur again until the flow has reattached. Plow reattachment 

takes place when the angle of attack drops below the angle Generally 

a dynamic stall angle about three degrees above the static stall angle gives 

good results. Different values of a , , Ail/ , , and a can be used for the 

ds ds re 

lift, drag, and moment characteristic if necessary to adequately model the 
dynamic stall of an airfoil. The calculation of the vortex induced lift in 
dynamic stall is outlined in figure 13. The drag and moment are calculated 
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Figure 13. Outline of calculation of dynamic stall 


vortex- induced lift coefficient. 


I 
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in a similar fashion, except that the drag is not multiplied by sign of 
(a (90 - j a |)); and 1/2 Ac^, is added to the moment in reverse flow. 

As an alternative dynamic stall model, consider that developed in 
references 12 to 14. They introduce an effective angle of attack of the form 



— 


WV</v% 





where i is a function of Mach number and the airfoil section obtained from 

Xj 

oscillating airfoil tests. This angle used in place of in 

the expression for c. given above, with Ac =0. Similar corrected angles 

of attack are calculated for the moment and drag coefficients, using appropri- 

ate factors x,, and x^- For an NACA 0012 airfoil, reference 14 gives 
M D 

^ C A/\ 


(in radians) . 


A no-stall model can be implemented by using for the smaller of the 

actual angle of attack a and a maximum angle of attack a in the linear 

^ max 

range (say 10^) : 




I vwx-fjd. (\oL\ ^ \’%cT — 


\^\ 


The incremental coefficients (Ac., Ac , and Ac.) should be set to zero as 

Jc m d 

well. 

In summary, the following procedure is used to calculate the section 

aerodynamic coefficients. First the Mach number correction for tip flow is 

applied: Th^ section coefficients c^, c^, and c^ are calculated 

from a, a, A, and first the yawed/delayed effective angle of attack is 

calculated; then c.._, c,^_ and c ^.^for a __ and M -- are obtained from 

£2D d2D m2D eff eff 

two-dimensional steady airfoil tables; c^ 2 j) a = 0 or 180° is also required; 
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finally the section coefficients are evaluated • Next the dynamic stall vortex 
loads evaluated from a and a. Finally the 

tip loss correction is applied to the lift coefficient. 

The aeroelastic analysis (see section 6.1.4) requires the derivatives of 
the section coefficients with respect to angle of attack and Mach number. 

These derivatives are evaluated by applying the above procedure with small 
increments in a and M (not the purpose of evaluating these 

derivatives, a, A, and the dynamic stall vortex loads are held constant. 

2.4.8 Unsteady lift and moment^-- The thin airfoil theory result for the 
unsteady aerodynamic lift and moment about the pitch axis for the rotairy wing 
is 


where x^^ is the distance between the aerodynamic center and the elastic 
axis: 






'•■K. 



normal flow 
reverse flow 


(here x^^ must in fact be the position of the aerodynamic center) ; and in 
the double sign the upper one is for normal flow and the lower one for 
reverse flow, ± - sign (V). Here w = u^ sin 6 - Up cos 0 is the upwash 
velocity normal to the blade surface (with no order c terms); B = dw/dx is 
the gradient of the upwash along the chord, as due to a pitch rate; and 
V = Up cos 0 + Up sin 0. 
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In this result the order c lift and order c moment terms have been 
neglected. The virtual mass terms (aerodynamic forces due to the section 
pitch and heave acceleration) can also be neglected. The sign changes in 
reverse flow have been accounted for in this result. Radial flow effects are 
included in the slender body pressure terms (from the radial derivative w') 
and in the contributions to the upwash w. The time derivative w includes 
terms due to the time varying free stream. Corrections for real flow effects 
on the lift“Curve slope and aerodynamic center have been included (thin 
airfoil theory gives a = 27 t and the aerodynamic center at the quarter chord) . 

For stalled flow, these unsteady aerodynamic forces can be set to zero 

(L^^ = ~ 0) . The unsteady forces at high angle of attack are accounted 

for in the dynamic stall model for c^, c,, and c . 

£ d m 

Finally, the velocities required for these unsteady aerodynamic forces 
are as follows: 

V = -4- & 

B= © -•- (ici. -CeoT? 

• . 

and from w = u^ sin 6 Up cos 0, Up * z + ” ^"^o ^ * 

there follows 

tA-|* ^ 

V ( & ^ ^ ^ ^“T ^ 

-4- U|J_ 
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with 




= u cos - y sin 
X y 


and using 


^ fa 4-T^«(-t--^ — Vo“?-V 

^ - Vo"^ 
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2.4.9 Circulation.- The blade bound circulation is required for the 
wake induced velocity calculation: 

r = % u«=-cj, -V 


Thin airfoil theory gives for the unsteady circulation (below stall) 

ik fe Cl 


(see section 2,4.8) • 


2 . 5 Environment 


The aerodynamic environment of the helicopter is defined by the speed of 

sound c , and the air density ratio to sea level standard p/p . (The blade 
^ o 

Lock number is calculated using p^.) One approach is to input values of c^ 
and p/p^. Alternatively these parameters can be calculated from the altitude 
h for a standard day; or from the pressure altitude and temperature x. 


For a constant temperature lapse rate, the density ratio and speed of 
sound are obtained from the following expressions. 




"Ti/T-o 

c T-, -t-'cbZ-n, 


standard day 
given temperature 


A- ^ 
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For the case of temperature and altitude specified, the density altitude is 





Alternatively, the air density and temperature can be specified directly. 
Then the equivalent altitude can be obtained from 







1 

/ A 3 

Mo To ) 



The required constants are given in the table below. 


Constants English units SI units 


dimension h 
dimension x 
g/5R 

(g/CR-1)"^ 

Tq/5 



g 


(g/5R) 


ft m 

Oj o 

5.256115 

0.234956 


145442 ft 

518.67 °R 

459.67 °R 
1116.45 ft/sec 

32.17405 ft/sec^ 
.002378 slug/ft^ 


44330.8 m 
288.15®K 
273.15'^K 
340.294 m/sec 
9.80665 m/sec^ 
1.225 kg/m^ 


0.190255 
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2.6 Normalization Parameters 


It has been the practice here to deal with dimensionless quantities 

based on the air density, rotor speed, and rotor radius (p, Qy and R) . 

In addition, the equations of motion, the inertial coefficients, and the 

aerodynamic forces have been normalized using the following parameters: 

I , a characteristic moment of inertia of the blade; c , the blade mean chord 
b Tn 

and a, the blade two-dimensional lift-curve slope. The values of these 
parameters have no influence on the numerical problem and its dimensional 
solution; they only affect the values of normalized, dimensionless quantities 
It is convenient to use the blade Lock number y the rotor solidity a 

as primary parameters. Then and c^ are obtained from 



For this purpose, the lift curve slope is set to a value of a - 5.7. 
The Lock number will be defined for standard sea-level conditions (y^) ; 
then y = y (p/p ) . 
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3. ROTOR WAKE ANALYSIS 


3.1 Nonuniform Wake- Induced Velocity 

3.1.1 Rotor vortex wake.- Conservation of vorticity on a three- 
dimensional wing requires that the bound circulation is trailed into the 
wake from the blade tip and root. Radial variation of the bound circulation 
produces trailed vorticity in the wake, parallel to the local free stream 
direction at the instant it leaves the blade. Azimuthal variation of the 
bound circulation will produce shed vorticity, oriented radially in the wake. 
The strength of the trailed and shed vorticity is determined by the radial 
and azimuthal derivatives of the bound circulation at the time the wake 
element left the blade. The lift and circulation are concentrated at the tip 
of the rotating wing, due to the larger dynamic pressure there. Consequently 
the trailed vorticity strength is high at the outer edge of the rotor wake, 
and the vortex sheet quickly rolls up into a concentrated tip vortex. The 
Lormation of this tip vortex is influenced by the blade tip geometry. With 
.oiuare tips, much of the roll up has occurred by the time the vortex leaves 
^he trailing edge. The rolled up tip vortex quickly attains a strength nearly 
equal to the maximum bound circulation of the blade. The tip vortex has a 
small core radius, depending on the blade geometry and loading. The vorticity 
in the tip vortex is distributed over a small but finite region, called the 
vortex core, due to the viscosity of the fluid. The vortex core radius is 
defined at the maximum tangential velocity. The vortex core is an important 
factor in the wake induced velocity, since it limits the maximum velocity 
induced near a tip vortex. Only a limited amount of data on the vortex core 
radius is available, particularly for rotary wings. There is an inboard 
vortex sheet of trailed vorticity in the wake, with opposite sign as the tip 
vortex. Since the gradient of the bound circulation is low on the inboard 
portion of the blade, the root vortex is generally much weaker and more 
diffuse than the tip vortex. 

The trailed and shed vorticity of the rotor wake is created in the flow 
field as the blades rotate, and then convected with the local velocity of the 
fluid. This local velocity consists of the free stream velocity, and the 
wake self-induced velocity. The wake is transported downward, normal to the 
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disk plane, by a combination of the mean wake induced velocity and the free 
stream velocity. The wake is transported aft of the rotor disk by the inplane 
component of the free stream velocity. The self-induced velocity of the wake 
also produces substantial distortion of the vortex filaments as they are 
convected with the local flow. Thus the wake geometry basically consists of 
distorted interlocking helices, one behind each blade, skewed aft in forward 
flight. 

The strong concentrated tip vortices trailed in helices from each blade 
are the dominant feature of the rotor wake. Due to its rotation, a rotor 
blade encounters the tip vortex from the preceding blade in both hover and 
forward flight. These tip vortices produce a highly nonuniform flow field 
through which the blades must pass. In hover the tip vortex is convected 
downward only slightly until after it encounters the next blade. The vortex 
produces a large variation in the tip loading on the following blade there- 
fore, with a substantial influence on the rotor hover performance. In for- 
ward flight the rotor wake is convected downstream, so the tip vortices are 
swept past the entire rotor disk instead of remaining in the tip region. The 
close vortex/blade encounters occur primarily on the sides of the disk, where 
the blades sweep over the vortices. The resulting large azimuthal variation 
in the induced velocity produces a large higher harmonic content of the blade 
loading. Nonuniform inflow is thus an important factor in the vibration, 
loads, and noise of the rotor in forward flight. In a tandem helicopter, 
the rear rotor also encounters the wake of the front rotor. 

For close vortex/blade encounters, the induced loading varies rapidly 
along the blade span. Lifting line theory does not give an accurate predic- 
tion of such loading. Thus lifting surface theory is required to accurately 
estimate the vortex-induced loads on a rotary wing. The most economical 
approach is to use lifting line theory with a correction factor for close 
vortex/blade encounters, based on a lifting surface solution for an infinite 
aspect-ratio, nonrotating wing encountering a straight, infinite, constant- 
strength free vortex. In the present analysis this correction will be 
incorporated as a factor reducing the induced velocity as required to obtain 
the correct loading by lifting line theory. Note however that with this 
approach the actual blade angle of attack at vortex/blade interactions will 
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be larger than calculated for the lifting line theory loading solution. 

Direct application of lifting surface theory to the rotary wing is usually 
impractical with current computation techniques and machines. An examination 
of measured rotor airloads indicates that the vortex Induced loading is 
generally high when the blade first encounters a vortex, but decreases as the 
blade sweeps over the vortex. There is evidently some phenomenon limiting 
the loads (see reference 15 to 17) . Local flow separation due to the high 
vortex— Induced radial pressure gradients on the blade appears at present to 
be the most likely explanation for the reduction in loading after the initial 
encounter. Bursting of the vortex core induced by the blade is also a possi- 
bility. Another possibility is that the vortex interacts with the trailed 
wake it induces behind the blade, with the effect of diffusing the circulation 
in the vortex. Note that the latter two phenomena, involving a change in the 
vortex Itself, will also influence the loading if the vortex encounters yet 
another blade of the rotor. Following reference 16, the phenomenon limiting 
vortex induced loads after the initial encounter will be modelled by increas- 
ing the core radius of a segment after it encounters the blade, with upstream 
propagation along the vortex to produce the loads reduction. An increase in 
core size is a convenient means to reduce the influence of the vortex; the 
exact physical explanation for this phenomenon is at present speculative. 

A possible model for the tip vortex viscous core is solid body rotation, 
which implies that all the vorticity is concentrated within the core radius 
’^c at the point of maximum tangential velocity) . Measured vortex 

velocity distributions show that the maximum tangential velocity is much less 
than r/2ir r^, implying that a substantial fraction of the vorticity is out- 
side the core radius. Reference 16 suggests using a circulation distribution 


based on measured velocity distributions of vortices from nonrotating wings; 
the corresponding vorticity distribution is 
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where r is the distance from the vortex line. In this case half the 
vorticity is outside the core radius. Along a line at right angles to the 
vortex and a distance h above it (as in a blade/vortex intersection) , this 
vorticity distribution produces a downwash with peak value 



a distance \/ h + r either side of the intersection; compared to w 
“ r/47Th at a distance h from the intersection for a vortex with no core. 
(Note that as far as the downwash velocity is concerned, this core effect is 
equivalent to moving the vortex away from the blade, to an effective distance 

n T 

he ■“ such a simple interpretation will be useful in the lifting 
surface correction.) The peak tangential velocity with this vorticity 
distribution is F/Att r^* half the value obtained with all the vorticity 
concentrated within the core radius. 

The rotor wake induced velocity is calculated by integrating the 
Biot-Savart law over the vortex elements in the rotor wake. The wake 
strength is determined by the radial and azimuthal variation of the bound 
circulation. For the wake geometry a simple assumed model, experimental 
measurements, or a calculated geometry can be used. With the helical geometry 
of the rotary wing wake, it is not possible to analytically evaluate the 
induced velocity, even if the self-induced distortion of the wake is neglected. 
A direct numerical integration of the Biot-Savart law is not satisfactory 
either, because the large variations in the induced velocity at close vortex/ 
blade encounters requires a small Integration step size for accurate results. 

It is most accurate and most efficient to calculate the rotor nonuniform 
inflow with the wake modelled using a set of discrete vortex elements. For 
each vortex element in the wake the induced velocity at a point in the flow 
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fxeld is evaluated by an analytical expression, and the total induced velocity 
is obtained by trimming contributions from all elements. The tip vortex is 
well represented by a connected series of straight-line vortex segments. The 
inboard trailed and shed vorticity can be modelled using rectangular vortex 
sheets, or a lattice of discrete straight-line vortex segments (with a large 
effective core to limit the Induced velocity close to individual line seg- 
ments) . A large core vortex line element might in some cases be a better 
model than a sheet for the inboard trailing vorticity, if the inboard wake 
has partially rolled up to form a root vortex. The inboard wake is less 
important to the nonuniform Inflow calculation than the tip vortices, so a 
more approximate model may be used. The approximations involved in modelling 
the rotor wake using a set of discrete vortex elements include replacing the 
curvilinear geometry by a series of straight-line or planar segments; a 
simplified distribution of vorticity over the individual wake elements 
(linear variation, or even constant strength); and perhaps physical approxi- 
mations such as the use of line elements to represent the inboard vortex 
sheet. The development of a practical model involves a balance between the 
accuracy and efficiency resulting from such approximations. 

3.1.2 \Jake model. - The blade bound circulation will be calculated at 
discrete points on the rotor disk radially and azimuthally. Assuming a 
linear variation of the bound circulation between these known points results 
in a wake model consisting of vortex sheet panels (see fig. 14) . Assume that 

the blade bound circulation F(r, ^p) is given at the radial stations 

r^ (i = 1 to M) and at the azimuthal stations ij). = jAij; (j = 1 to J, 

Ai); = 2ir/j) . Let (j) be the age of vortex elements in the wake (((i^^ = kAi^, 
k = 0 to ») . The strength of the trailed and shed vorticity of a wake element 
is determined by the bound circulation of the blade at the time the vorticity 
was created. Consider a wake panel of age <p = to arising from the 

blade between radial stations r^ and r^^^ (fig. 15). The strength of the 
vorticity in this panel is determined by the bound circulation at the time 
it was created, which is known at the four comers. The bound circulation 
corresponding to the panel leading edge is that at time ’I' “ where is 
the current blade position (dimensionless time) and = 1tAi|i is the age 
of the panel at the leading edge. The bound circulation corresponding to the 
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Figure 14. Wake model with bound circulation calculatec? 
at discrete points on rotor disk. 
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panel trailing edge is that at '1^ - ^k+1’ earlier than the leading edge. 

The difference between the bound circulation at r. and r^ . ^ defines the 

1 i+1 

trailed vorticity strength 6, which is constant radially along the panel 
assuming a linear variation of the bound circulation from r^ to When 

the bound circulation varies azimuthally however the trailed vorticity 
strength 6 is different at the panel leading and trailing edges; a linear 
variation of 6 in the direction of the trailed vorticity will be used. 
Similarly, the difference between the bound circulation at ^ and 

'P defines the shed vorticity strength y, which is constant azimuth- 

ally along the panel (for a linear azimuthal variation of the bound circula- 
tion) but varies linearly from the left to the right panel edges. 

A vortex sheet panel in the wake may be economically approximated by 
shed and trailed line vortices located in the middle of the panel, with a 
large core to avoid the induced velocity singularity near a vortex line. A 
vortex lattice model of the rotor wake is produced by collapsing all the 
wake panels to such finite strength line segments. Since the line segments 
are in the center of the sheets, the points at which the induced velocity 
and bound circulation are evaluated lie at the midpoints of the vortex 
lattice grid, both radially and azimuthally. Positioning the collocation 
points midway between the trailed vortex elements (radially) is a standard 
practice of wing theory utilizing the vortex lattice wake model, in order to 
avoid the singularities at the lines; positioning the collocation points mid- 
way azimuthally is required to correctly obtain the unsteady aerodynamic 
effects of the shed wake (see ref. 18) . Simply collapsing the shed and 
trailed vorticity in the wake panels to lines, the strength of the line seg- 
ments will vary along their length as described above. The shed and trailed 
line segments will cross in the middle of the panel. As a further approxima- 
tion, a stepped (piecewise constant) variation of strength can be used 
instead of the linear variation, with the jump in strength occuring at the 
center of the segment where it crosses the other vortex line. Such a vortex 
lattice wake model with constant strength line segments corresponds to a 
stepped distribution of the blade bound circulation, azimuthally and radially 
(with the jumps occuring midway between the points where the circulation is 
evaluated.) 
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The rotor vortex wake quickly rolls up at the outer edge to form a 

concentrated tip vortex. Because of the dominant role of the tip vortex in 

the wake flow field, it is important to model these rolled up tip vortices in 

the induced velocity calculation. The lesser role of the inhoard wake vortic— 

ity also allows a more approximate model to be used for it. Let T (tIj) be 

max 

the radial maximum of the blade bound circulation. It is assumed that in the 
far wake, where the rollup process is complete, that all of the bound circula- 
tion ^max concentrated in the tip vortex. The tip vortex will be 
modelled by a vortex line segment with a small but finite core radius. When 
^max azimuthally, the tip vortex strength varies along its length. 

Furthermore, the inboard portion of the wake will be modelled by a single 
sheet panel, with trailed and shed vorticity as described above. This far 
wake model may be viewed as corresponding to the circulation distribution 
sketched in figure 16. The linear variation from F = 0 at the root to 
^ ~ ^max defines the single inboard sheet, and the sharp drop 

^max defines the tip vortex line. (This circula- 

tion distribution should not be associated with the actual bound circulation 
at the rotor blade. Rather it is an approximation for the vorticity distribu- 
tion in the far wake, which is determined by the rollup process. Since an 
analysis of the rollup is not attempted here, the actual vorticity distribu- 
tion over the inboard sheet is not known. An approximation involving constant 
strength determined by the known maximum bound circulation is appropriate 
therefore.) This far wake model is computationally efficient, since it 
depends only on the maximum bound circulation F 

max 

The rollup process may not be complete by the time the tip vortex 

encounters the following blade. The induced loads will be significantly 

lower if the tip vortex has strength less than the maximum bound circulation. 

Therefore the tip vortex rollup must be included in the wake model. Figure 16 

sketches the radial circulation distribution assumed, which produces the 

model for the rolling up wake. The circulation goes from zero at the root 

^max radial station r^^; to f^^^^ F^^ at the tip. Thus there is a 

line tip vortex of strength f F , and two inboard wake panels . The 

rollup process will take place over the wake from (j) = 0 to (|» = (|> . The 

RU 

position of the maximum circulation and the rollup fraction will vary linearly 
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^MAX 



ROLLING UP 
WAKE 



NEAR WAKE 


Figure l6. Equivalent circulation distribution for models 
of far wake, rolling up wake, and near wake. 
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from and at (|i - 0 to r = 1 and f = 1 at 4> = ^ analysis 

of the rollup process is not part of the present work, so the parameters 
'*’rU’ ^RU’ ^RU prescribed Inputs to the calculation procedure. 

Note that lae velocity induced by the rolling up wake will also depend only 
on the single bound circulation value I 

max 

Just behind the reference blade, where the induced velocity is being 
calculated, it is the detailed radial and azimuthal variation of the wake 
vorticity which is important, not the rollup process (except for the influence 
of the rollup on the tip loads) . Hence for the near wake of the reference 
blade the full vortex panel representation is retained. The corresponding 
radial distribution of the circulation is also sketched in figure 16 for the 
near wake; in this case it is the actual blade bound circulation distribution. 
The tip vortex rollup is often partially complete at the blade trailing edge, 
so a line vortex at the tip is included, with strength equal to a fraction 
^NW calculated bound circulation at the most outboard radial station, 

-he complete model of the rotor wake is shown in figure 17. 

The very first panels of the near wake require special consideration. 

In order to correctly calculate the unsteady aerodynamic effects, the shed 
wake is stopped a quarter chord behind the bound vortex (ref. 19). The 
singularity near the side edges of the trailed vortex sheets presents a 
difficulty in calculating the induced velocity at a point due to the immedi- 
ately adjacent panels. Thus if the induced velocity is to be calculated near 
a junction between two panels, they should be replaced by one panel with the 
collocation point well away from the edges of the single panel. This diffi- 
culty can be also avoided by using line vortex elements for the trailed 
vorticity in the near wake, or by moving the panel side edge away from the 
collocation point. Finally, the front edges of the individual panels should 
all be aligned with the bound vortex. 

When calculating the Induced velocity at points off the rotor blade, as 
at another rotor or for the airframe aerodynamics, the near wake model need 
not be used. Often calculating the induced velocity away from the rotor will 
require a consideration of more wake spirals than are needed for points on the 
rotor disk. 
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TIP VORTEX 



Figure 17 . Sketch of wake model for nonuniform induced 
velocity calculation. 
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For computing the Induced velocity, the tip vortices will be represented 
by a connected series of straight vortex line segments, with a small viscous 
core radius. Normally a linear variation of the strength along each segment 
will be used; a piecewise constant variation corresponding to a stepped bound 
circulation distribution is also possible. The inboard wake panels will be 
represented by planar, rectangular vortex sheets, with shed and trailed 
vorticity varying linearly along its length. For computational economy, the 
vortex sheets can be replaced by line segments in the middle of the sheets, 
with a large core radius (which in this case does not have physical signifi- 
cance, rather a viscous core is a convenient means for eliminating the 
singularity near a line used to represent a sheet element; unless the inboard 
trailed vorticity does rollup to form a diffuse root vortex) . If the induced 
velocity is to be calculated near the side edge of a vortex sheet element, it 

can be replaced by a line element in order to avoid the edge velocity singu- 
larity. 

3.1.3 Geometry.^ A nonrotating tip path plane coordinate frame with 
origin at the rotor hub will be used for the Induced velocity calculation. 

The solution process will iterate between the Induced velocity calculation, 
and the harmonic blade motion and helicopter trim solution (using uniform 
Inflow to start the cycle) . Thus the hub plane orientation (S system) will 
be updated based on a new induced velocity estimate. In contrast, the tip 
path plane orientation is well defined by the helicopter or rotor trim, hence 
is less sensitive to changes in the induced velocity estimate. Also, the 
rotor wake geometry is simplest when defined relative to the tip path plane. 
The tip path plane tilt relative to the hub plane is given by the first 


harmonics of the tip deflection 


fl‘^(l) : 


f ^ 
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The rotation matrix 



will transform position and velocity vectors from the shaft axes (S system) 
to the tip path plane axes (T system) . 

The Induced velocity is required at the radial stations r^ along the 
rotor blade. The radial stations at which the induced velocity and bound 
circulation are evaluated will each be a subset of the blade loading radial 
stations, but the two sets need not be identical. From section 2.2.4, the 
postion vector of the rotor blade is 

-V V -i- 9w.h? '^5'^ X-- 
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Transforming to the tip path plane (neglecting second order terms) gives 
then 

^ ^ + roasN'^ 

— c(^c -♦'^s ^ 


The induced velocity is to be evaluated at the points r7(r ) along the 
rotor blade. It is useful in the computation of ^ to have the option to 
suppress the inplane deflection, to suppress all harmonics except the mean, 
or to linearly interpolate the geometry between the root and the tip. 

The wake induced velocity is also required at points in the flow field 
. ?f the rotor disk: 

a. at the wing/body, horizontal tail, and vertical tail for the rotor/ 
airframe aerodynamic interference; 

b. at the other rotor hub for rotor /rotor aerodynamic interference; 

c. at an arbitrary point in the flow field; 

d. and at the reference blade of the other rotor, for detailed rotor/ 
rotor aerodynamic interference. 

For the first two, only the mean value of the Induced velocity will be used 
in the present analysis. The induced velocity distribution over the disk of 
the other rotor can be used in the present analysis only if the two rotors 
have the same rotational speed (see section 5.1.11) ; so for the single main 
rotor and tail rotor configuration the rotor/rotor interference can be 
accounted for only in terms of the induced velocity at the rotor hub. The 
position vector of the wing/body is 
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where and are the position of the rotor hub and the wing/body center 

of action, in the body axes (F system). The vectors for the horizontal tail, 
vertical tail, and other rotor hub are similar. For the induced velocity of 
rotor //2, these vectors must be multiplied by For the induced 

velocity at the disk of rotor //2 due to the wake of rotor #1, the position 
vector is 

^2/1 “ L" 


where ^'1^21 azimuth angle of the reference blade of rotor #2 when 

= 0 for rotor #1 (see section 5.1.5) and 
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The position vector for the induced velocity at the disk of rotor #1 due to 
the wake of rotor //2 is 



(Tt-r, -I- i=^, 





Cc-„. 


^©"tVoor* 




If the rotor rotates clockwise it is necessary to change the sign of the 

\ ® 

component of r (between the R^^ and R^^ rotation matrices) . 

The geometry of the tip vortex behind the reference blade will be defined 

by the vector r <J>), where \\) is the present azimuth angle of the blade 

w 

and (J> is the age of the vortex element. The wake geometry is required at 
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ths dxsciT 6 t 0 azimuth positions ® £rAi{/ and walc6 agas (p^ ® V Ai|i whsrB Z 

ranges from 1 to J (one revolution of the blade, with = 2ir/j) and k 

ranges from zero to the specified number of wake spirals for the induced 

velocity calculation. The tip vortex geometry behind the other blades of 

the rotor can be obtained for r^ at the appropriate azimuth angle. The 

tip vortex elements are created at the blade tip (r^ at radial station 
.. b 

r = 1) , convected with the free stream velocity and distorted by the 
self— induced velocity in the wake. The rotation of the wing together with 
convection by the free stream velocity produces the basic helical geometry 
of the rotor wake. As at the rotor disk, the induced velocity throughout 
the wake is highly nonuniform. The actual position of the wake elements, 
determined by the Integral of the local convection velocity, is thus highly 
distorted from the basic helical form. The resulting wake geometry is 

^here D(4^, (|)) is the distortion due to the wake self-induced velocity (note 
0) = 0) and the free stream convection velocity is 


B f 1 


^ 1 -^y. 1 


\ j 



relative to the tip path plane. 

Similarly the geometry of the inboard wake sheet will be defined at the 

root and tip edges, trailing from the blade position T at radial stations 

b 

^ “ ’^ROOT r = 1 respectively. The distortion D will be different for 
the tip vortex and the inboard sheet. Because of the dominant role of the 
lip vortices, the most important information in the wake geometry is the tip 
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vortex position, and a less accurate definition of the inboard sheet geometry 
is often acceptable. 


The induced velocity calculation may require the wake geometry beyond the 

point where the stored distortion ends. For this portion of the wake rigid 

geometry will be used. Consider the distortion (})) when the age (J> is 

greater than the age of the last element in the known distortion, 6- 

last 

= The wake geometry will be extrapolated from using 

only vertical convection due to the mean induced velocity: 


Note that the azimuth angle of the blade at the time the wake element was 
created, ~ has been held constant. 

For the self-induced distortion of the rotor wake geometry, the following 
models are considered: 

a. rigid or prescribed wake, with contraction and two-stage convection; 

b. and a calculated free wake geometry (section 3.2). 

In the rigid wake geometry it is assumed that all elements in the wake are 
convected downward by the mean induced velocity at the rotor disk, giving 



relative to the tip path plane, where 


Note that this distortion is independent of the azimuth angle The con- 
vection velocity is the mean induced velocity at the rotor disk, includ- 

ing ground effect and rotor /rotor interference in general. This model can be 
generalized to a two-stage convection: 


- 14 ?- 



-k.+ 


N 

4>£= 

with 

- ^Ai 


determined by the constants and To improve convergence, 

X = 1/2 (X^ + should be used in place of X^. Including contraction 

of the wake gives for the distortion 



tr \ 

br ^ \ 

J 


where the radial displacement (also independent of is 

= - Cl-cT^*'^ >Cl-K^')r-. 


^ 1 the tip vortex and the outside edge of the inboard sheet , and 

® ’^rOOT inside edge.) Hence the rigid wake geometry is determined 

by the parameters f^, f^, K^, and K^, which may be different for the tip 
vortex and inboard sheet. Alternatively, the constants and can be 
specified, instead of f^ and f^. 

Landgrebe (ref. 20) developed a prescribed wake geometry model for a 
hovering rotor from experimental model rotor flow visualization data. The 
model consists of contraction and two-stage convection as defined above, with 
the constants as follows: 
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a. tip vortex 


k ^ ^ \ -4- • %/Crr ' 

b. sheet tip edge 

k, = 1-S5 7 ^ 

= \.=i 

c. sheet root edge 

K. = 

1/ — — ^«ooZ^ 


d. radial contraction 

Ks = 

K 4 = 


“H 2.1 C.-T 

-1^ 


where 6 is here the blade linear twist rate in degrees. Kocurek and 
tw 

Tangier (ref. 21) revised the tip vortex geometry based on experimental data 
for low aspect-ratio two-bladed rotors, obtaining 

k., =. B + c 

vl .. \/vv\ ^ 

^ [c^ -N C— J 

)C^ = 4.0 
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with 

-=- .ooo“7 2-^ d-tw 


c = 

2*3 


— 

\ *o 

— e. 

w =• 

o-e 

— - Ol“7 2. ©-rw 


In reference 22 it was found that the prediction of measured rotor hover 
performance was improved when the wake geometry was prescribed based on the 
blade maximum bound circulation rather than as above. Hence 








can be used in place of C^; in general must be averaged over ^ as 

-wcli. These prescribed wake models were developed for a hovering rotor. To 
apply them in more general flight conditions we can use 







with the mean induced velocity including ground effect and rotor/rotor 

interference (the blade loading C^/o is retained.) 

The wake geometry arrays will be organized as follows. The rigid or 
prescribed wake geometry is defined by and at <J) - kAi/;, k = 1 to 
(independent of tp) . So the structure of the array is 

D(k), for k = 1 to 

for and D^, for the tip vortex and the two inboard sheet edges. The 
convection rate is also required, for extrapolation of the geometry 

beyond (p = The free wake geometry is defined by 4>) . The first 
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subscript in the array will be the age (}) = for k = 1 to the 

second subscript will be the blade azimuth angle ip = ZAip, f. = 1 to J. So 
the structure of the array is 


D(n), n = (J? - l)Kp^.Q + k for ((k = 1 to , X = 1 to j) 


The free wake geometry will be used for the tip vortex only (see section 3.2). 


In the near wake and the rolling up wake, the position of a panel corner 

at an arbitrary radial station p is required (r < p < 1) . Linear 

ROOT 

interpolation between the root and tip edges of the inboard sheet gives 


- 




COoX 












The geometry of the near wake panels should include the increment 


Ar 




I- 


to account for the blade bending (the variation with wake age is neglected). 

The first panels of the near wake are aligned with the bound vortex. Let 
r^ and ^ be the position vectors of the right and left front corners of a 
near wake panel, obtained from the wake geometry at (J) = 0. Let and r^ 
be the position vectors of the right and left rear corners, obtained from the 
wake geometry at (p - Then a rectangular panel aligned with the bound 

vortex is obtained if and r^ are replaced by r^ + Ar and r^ + Ar 
respectively where 


Ar 




( n-r3~ > 


]• (' 
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For the shed vortlcity in these first panels, the front comers are also 
moved aft by a quarter chord, which is accomplished by adding 





to and with Ar given above. The blade chord at the induced velocity 
radial station will be used; and the shed wake panel will be omitted entirely 
if c/4 > |Arl . 

That portion of the first tip vortex segment extending from the bound 
vortex to the trailing edge (a length 3c/4) should be perpendicular to the 
bound vortex. Let r^ and r^ be the position vectors of the vortex segment 
end points, at the blade tip and at the first downstream point respectively. 
Let r^ be the position vector of the first blade point inboard of the tip. 

Then this line vortex segment will be replaced by two segments extending from 
to 

-a 3c 


and from there to 



where 


Ar 




and c is the tip chord. Constant strength will be assumed for the portion 
from the bound vortex to the trailing edge. 


Ij 


r(r^, \|ij) . 


Th^ solu- 


3.1. A lyiduc^d VGZocity ocLZoziZcii'tori^— The blade bound circulation is 
calculated at discrete points on the rotor disk: V 

tion is periodic, so the azimuthal points cover one revolution of the blade 
^or j = 1 to J (A\/; « 2 tt/J). The radial stations r. (i = 1 to M) 
will be a subset of the aerodynamic loading radial stations. Except for the 
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near wake, the vorticity strength in the present wake model actually depends 

only on the maximum circulation F., defined as the value of T,. with maxi-- 

J 13 

mum magnitude over all radial stations r^ at a given azimuth (the 

computation will allow the use of the maximum over the radial stations out- 
board of station r^ ) . 

Gmax 

Summing the contributions from all vortex elements in the wake gives 
the induced velocity as the product of the blade bound circulation and 
influence coefficients: 







H- 







The second term is due to the near wake (extending from ({) = 0 to ({) = 
behind the reference blade at azimuth angle = AAijj) • A set of influence 
coefficients is obtained for each point in the flow field at which the 
induced velocity is calculated: at points distributed radially and azimuthally 

over the rotor disk, and at points off the rotor disk (see section 3,1*3). 

The influence coefficient arrays will be organized as follows. Consider 
C.(r); the first subscript is the index due the azimuth angle of the bound 
circulation ( j = 1 to J) . The second subscript is the index over all the 
field points r* at a given azimuth angle (k = 1 to MR) . The third subscript 
is the index over the azimuth angle of the field points (J. = 1 to J) . So the 
structure of the array is 

for (((j = 1 to j), k = 1 to MR), 1 to j) 


The field points at a given azimuth angle consist of the induced velocity 

points along the rotor blade span; perhaps the induced velocity points along 
the blade of the other rotor, or at the hub of the other rotor; and perhaps 
the points at the wing/body, horizontal tail, vertical tail, or an arbitrary 
field point. The organization of the array for the near wake influence coef- 
ficients C.. (t is similar, except there is an additional subscript which 
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is the index over the circulation radial stations (i = 1 to M) , and the index 
over the azimuth angle of the bound circulation covers only the near wake 
O - >■ - <^0 l) : 

* + + i 

for ((((1 . 1 to M), j - Ji,- to ), k - 1 to MR). 

1 to j) 

Also, for the near wake the field points at a given azimuth station consist 
of only the induced velocity points along the rotor blade span (no points off 
the rotor disk) . 

The calculation of the influence coefficients proceeds as follows. The 
outermost loop involves the dimensionless time ij;, which is also the azimuth 
angle of the reference blade. The solution is periodic so the induced velocity 
is evaluated for ^ = 0 to 2tt (at the discrete points \1^ = £ = 1 to J, 

A'J> - 27t/J) , For a given the position vectors at which the induced 
velocity is required can be evaluated: at the radial stations along the 

reference blade; at the wing/body, horizontal tail, vertical tail, other 
rotor hub, or an arbitrary point; and at the radial stations along the 
reference blade of the other rotor. 

Next there is a loop over all the blades of the rotor (m = 0 to N-1; 
m = 0 is the reference blade) . The azimuth angle of the m-th blade is 

~ ^ 2tt/N = (5, + mJ/N)Ai|^. Finally there is a loop over the wake age 

4> = kAijj (k = 0 to the maximum extent of the far wake, which may be different 
when calculating the velocity at points on or off the rotor disk) . 

The blade specification plus the wake age determines the vortex panel 
being considered, extending from <J) to <() + Ai{/ behind the m-th blade. Given 
1/;^ and <p, the position vectors of this wake panel can be evaluated: the 

end points of the tip vortex line segment, and the four comers of the 
inboard sheet (at the side edges, as described in section 3.1.3). The wake 
strength at the panel leading edge is determined by the bound circulation 
at - <1>, and the strength at the trailing edge by the bound circulation 
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at - 0 “ Alp. These azimuth angles define to which influence coefficients 
the induced velocity of this panel contributes. 

The wake age determines whether the panel considered is part of the near 
wake, the rolling up wake, or the far wake models (as described in section 
3.1.2). The near wake model is only used behind the reference blade (m = 0) . 
The near way is not used in calculating the velocity at points off the rotor 
disk. 

The far wake model consists of a tip vortex line segment and a single 

inboard wake panel. The line segment has strength The sheet is due to 

a circulation distribution linear from zero at the inside edge to F at 

max 

the outside edge. The induced velocity expressions for a line segment and a 
rectangular vortex sheet then give the contributions to the influence coef- 
ficients . 

The rolling up wake model consists of a tip vortex line segment of 

strength fF , where 
max 


and an inboard wake sheet divided into two panels at radial station 




ffcu - Ci- 


Linear interpolation between the side edges gives the wake geometry at p. 
The circulation corresponding to these panels goes from zero at the inside 
edge, to F at p, to fF at the outside edge. The induced velocity 
expressions for these vortex elements then give the contributions to the 
influence coefficients. 


The near wake model consists of a tip vortex line segment of strength 


^NW % 


F„ . is the bound circulation at the most outboard radial 
Mj 

station) ; and separate inboard wake panels between the bound circulation 


radial stations. Linear interpolation between the side edges gives the wake 
geometry at r^ (i = 1 to M). An increment accounting for the blade bending 
must also be added to the position vectors in the near wake. The circulation 
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corresponding to these panels goes from zero at the inside edge, to at 

r^, to at the outside edge. The induced velocity expressions for 

these vortex elements then give the contributions to the influence coefficients. 

The inboard wake panels in the near wake directly behind the blade 
(4) = 0 to A4^) require special consideration (see section 3.1.2). First, the 
position vectors at the rear corners of each element are adjusted so the 
front edge of the rectangular vortex sheet is exactly aligned with the bound 
vortex. When evaluating the induced velocity near the junction of two trail- 
ing vorticity sheets, they can be replaced by a single sheet. The leading 
edges of the shed vorticity sheets must be moved a quarter chord behind the 
bound vortex. However, if line segments rather than rectangular vortex sheets 
are used for the trailing or shed wake, the above modifications are not 
required . 

Finally, the contribution of the bound vortex of each blade is calculated. 

The bound vortex is a straight line segment extending from the root to the 

tip of the blade at azimuth angle with strength varying from zero at the 

root to r at the tip. The contribution of the bound vortex of the 
max ^ 

reference blade to the induced velocity at the reference blade is not 
included. 


By this procedure the influence coefficients are calculated for a given 

wake geometry. Then from a circulation estimate at some stage in the blade 

motion and helicopter trim solution, the induced velocity *v can be evaluated. 

The vector v is at the induced velocity radial station, in the tip path 

plane coordinate frame. The aerodynamic analysis requires the induced velocity 

at the loads radial stations, in the hub plane coordinate frame 

(X=Xi - XT - Xk). Transforming to the S system gives 
X s y s z s 
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using the current values of the tip path plane tilt angles (6 and B ) . The 

c s 

induced velocity is then calculated at the loads radial stations by linear 
interpolation. 

The induced velocity for the rotor /airframe aerod 3 maraic interference (at 
the wing/body, horizontal tail, vertical tail, or an arbitrary point) is 
required in the body axis system (the F frame) : 



For the induced velocity from rotor //2, a factor is required as 

well. The induced velocity for the rotor/rotor interference (at the hub or 
over the disk) is required in the shaft axis system of the other rotor. The 
interference velocity at rotor //2 due to the loading of rotor #1 is then 



and at rotor #1 due to rotor //2 




These coordinate rotations normally should not be included in the influence 
coefficients because the updated values of the tip path plane tilt angles are 
to be used in the matrix R^^ (although the rotation of the induced velocity 
by the small angles and B^ should not be important) . The factor ^^R^/^^R 2 
can be incorporated in the influence coefficients however. For a clockwise 
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rotating rotor, the sign of the component of the induced velocity must 

be changed (between the and rotation matrices) . 

The interference induced velocity due to rotor #1 is calculated at 
rotor #2 for the azimuth angles = S,Ai|; + (^ = 1 to J) . When 0 

it will be necessary to linearly interpolate the velocity to the azimuth 
angles ip = kAt|i (k = 1 to J) . Similarly, the induced velocity due to rotor it 2 
at rotor ifl must be interpolated from to ip - kA\p. The 

velocity moved off the rotor disk by rotor #2 is also calculated for 
ip = £A\]; - and must be interpolated to \p = kAi|/. 

3.1.5 Ground effect.- Ground effect can be included in the nonuniform 
induced velocity calculation by introducing an image element for every 
vortex element in the rotor wake. The image element position is obtained 
by reflecting the actual wake element position about the ground plane, and 
changing the sign of the vorticity. Let be the distance the rotor 

hub is above the ground (see section 2.4.3). The position of the image element 
is required in the tip path plane axes relative to the hub. First the posi- 
tion vector of the actual element is rotated to earth axes; then the origin 
is shifted to the ground, the sign of the component is changed, and the 

origin is shifted back to the hub; finally the vector is rotated back to the 
tip path plane axes: 

= L2.«A6u‘'^ +■ 3 

where 

The Induced velocity of the vortex element is calculated, and subtracted from 
the induced velocity contribution of the actual element. The actual element 
is below the ground plane if 
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(This can occur if the wake geometry does not allow for ground effect..) In 
this case the induced velocity contributions of both the actual element and 
its image are set to zero. 

3.1.6 Hover or vertical flight (axisymmetrie geometry) The nonuniform 
inflow calculation can be simplified in hover due to the axisymmetry of the 
wake geometry. For the hover case the influence coefficients will be the 
same for the induced velocity at all azimuth angles, except for an azimuth 
shift and axis rotation: 


e. = 








I 


Even in hover the rotor may have a net pitch or roll moment if the center of 
gravity is offset from the shaft (with offset hinges or a hingeless rotor) . 
Hence in general the hover case will not involve induced velocity and bound 
circulation independent of azimuth angle. These considerations apply to the 
general vertical flight case as well. 

An accurate calculation of the induced velocity of a rotor in axial 
flight usually requires consideration of the wake very far from the rotor 
disk. The detailed wake model described above is required only close to the 
disk however. Very far from the disk a more approximate and more efficient 
model will be used, obtained by spreading the vorticity vertically over the 
distance h between successive sheets, as sketched in figure 18. The axial 
convection velocity in the far wake is taken from the prescribed wake model: 

giving for the spiral axial spacing 

iis = 'Z»tr V 
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Figure 18. Far waJce model for hover or vertical flight. 


( 
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The tip vortex elements are spread vertically to form a vortex sheet with 
axial and spiral components. There is a corresponding axial root vortex from 
the inboard trailing vorticlty. The shed vorticity is spread vertically to 
form a vortex sheet. This wake model extends L turns (an axial distance Lh) 
beyond the last spiral of the detailed wake model from each blade. 

The influence coefficients of this very far wake model are calculated as 
follows. For a panel in the last spiral in the wake, the geometry is 
specified by the location of the ends of the tip vortex line segment, and the 
corners of the inboard panel. The wake strength is determined by the bound 
circulation corresponding to the panel leading and trailing edges: 

The geometry of the sheet vorticity on the wake boundary is obtained from 

- V 

the position of the tip vortex segment, r^^ and (vectors to the segment 

ends) . The vectors to the sheet vector are then 





The induced velocity of the tip axial vorticity is given by the trailed sheet 
vorticity solution with 
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(see section 3.1.8). The induced velocity of the tip spiral vorticity is 
given by the shed sheet vorticity solution with 


r; = 1-r^ 

= o 
*• o 


Let r^ and r^ be the position vectors of the Inboard wake panel (the Inside 
and outside edges respectively, at + Ai|;/2) . Then the vectors to the 
inboard shed vortex sheet are: 

H 

Cj = Q - -1- q .11 

^ ft - 


and the induced velocity is given by the shed sheet vorticity with 


r. = ur; 
q = uq 

v; = Lf;.. 


(see section 3.1.8). Finally, the vectors to the ends of the root axial 
vortex are 
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7^ — LXV. 


and the induced velocity is given by the line vortex solution with 




(see section 3.1.7). 

3.1.7 Finite length vortex tine element*- Consider a straight vortex 
line segment of length s, as shown in figure 19. The vortex segment has 
linearly varying circulation, between and at the end points (T^ is 
the circulation at <f> and is the circulation at (j) + Aip) . The induced 
velocity is required at the point P, defined by the position vectors and 
r^ from the ends of the segment. The vectors *r^ and r^ can be in any con-- 
venient coordinate frame; the components of the induced velocity will be 
obtained in the same coordinate frame. The Biot-Savart law gives the induced 
velocity due to this line segment: 



where r is the vector from the element da on the segment, to the point P; 
and r = I r |. The coordinate a is measured along the vortex segment, from 


to S 2 : 


S. — 




— C 






= S, S 


where s is the length of the segment: 

T \ ^ ^ i X IL X — * 
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Write r = r - ae, where r is the minimum distance from the vortex line 
m m 

(including its extension beyond the end points of the segment) to the point 


and e 

is the unit vector in 

the direction of the vortex: 






s 

A 

e 






— i 

! vectors r and e are perpendicular, and 


s- = c;- 


The vortex strength varies linearly along the segment: 

s 


r ^ 


It follows that 


/iw — 




s 

r„ 










[" ^ 

L Fr 


- '1 
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>»aL 




n Kf, 






-4^-rr ^ <"v 
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for the induced velocity of the vortex line segment with linearly varying 
circulation. 


Consider also the induced velocity of a line vortex segment with a 
stepped circulation distribution. The distance from the midpoint of the 
segment to the point P is 




and the midpoint is located at 


S 


3 


( it 



The line segment has constant strength from s^^ to s^, and 
strength F^ from to 82* Applying the above result (with 
both values of the line segment gives 



^x*^3 r 











constant 

F = 0 ) to 

6 



The influence of the vortex core is accounted for by multiplying the 
induced velocity of the line segment by the factor 

yvk v/v\ ^ ^ 1. ^ 


for concentrated vorticity (solid body rotation) 



or by 
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for distributed vorticity (from ref. 16). The core radius r is the loca- 

c 

tion of the maximuni tangential velocity* 

In references 23 and 24 a lifting surface theory solution was developed 
for the vortex induced loads on an infinite aspect-ratio, nonrotating wing 
encountering a straight, infinite vortex at an angle A with the wing 
(fig. 20). The vortex lies in a plane parallel to the wing, a distance h 
below it, and is convected past the wing by the free stream. The distortion 
of the vortex line by interaction with the wing is not considered. In linear 
lifting surface theory, the blade and wake are represented by a planar distri- 
bution of vorticity. This model problem was solved for the case of a sinus- 
oidal induced velocity distribution, with wave fronts parallel to the vortex 
line. An approximate solution was obtained by fitting analytical expressions 
to the numerical solution for sinusoidal loading. The vortex induced velocity 
distribution can be obtained by a suitable combination of sinusoidal waves of 
various wavelengths, and the same super-position gives the vortex induced 
loading from the sinusoidal loading solution. The approximate solution is 
not valid for extremely small wavelengths, but the range of validity is 
sufficient to handle the cases arising in rotary wing applications. For the 
velocity induced along the wing span by a vortex of strength F: 



^ je;*- 


(where r and h are here divided by the wing semichord) the approximate 
lifting surface solution (from ref. 23) for the section lift is: 

" 1 Cr< + Ca.-4-CoV 
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Figure 20. Lifting surface theory solution for vortex- induced 
loads . 
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For the incompressible case, the coefficients in this expression are functions 
of the vortex angle A: 


= 5.12 + 1.88 (A/90^) 
b^ = -cosA 
b£ = J (A/90° - 1) 

I 

= .544 (-cosA) 4 .07 sin2A 

= -.434 - 1.09 (1-sinA) 4- .607 (1-slnA)^'^^ 

a^ = .0084 4- .0069 (-cosA)^’® 

c’ = 5.9 
o 

= 1.683 4 .27 (1-sinA)*^® - .154 (1-sinA)^'^ 

= 1.417 4 .366 (1-sinA)'®'^ - .392 (1-sinA)^'^ 

= .91 4 .93 (1-sinA)^*® - 1.025 (l-sinA)^*''*^ 


(A is in the range 90® to 180®; the solution for A from 0 to 90® is 
obtained by symmetry considerations) . The corresponding lifting line theory 
solution for the vortex induced loading is 

, _ MP i 

J 1 +- Oi.-vCo')*’ 


— cx 


i 
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~ ~«662, Cj^ “ 1«296( and •• tt/ 2. This lifting suirface solution 
will be used in the present analysis in the following manner. For each line 
segment it will be determined whether it is close enough to the blade for 

surface effects to be important is more economical to apply such 

a test than to always use the correction). If so, the induced velocity con- 
tribution of the line segment as calculated above, will be multiplied by the 
ratio 

The parameters required to apply this lifting surface correction for 

vortex Induced loads are h, rsin A, and A, From the minimum distance ~r ' 

m 

between the point P and the vortex segment (in the tip path plane coordinate 
frame) , and Including the influence of the viscous core on the Induced 
velocity, the vortex/blade separation is 

Av = f- 7 

.-A ^ 

Let i* « - r^)/s be the unit vector along the vortex segment; and 

j* = l!J,cos<;; + J^sinij; the unit vector along the blade. Then the intersection 
angle is 

JL = 

-A. yT I 

c^-a -A— = — \ ^ V 


Finally, the distance from the vortex line is 


r i'Wi -A. 




- 4 - ^ 
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The lifting surface correction will be used if the distance from the segment 
midpoint is less than a specified distance d : 

X/S 

Cj = I < 


Typically d^^ should be around 10c (see ref* 24) • 

The use of a larger viscous core radius after a blade/vortex interaction 
will be allowed, to model such affects as vortex induced stall or core burst- 
ing, which limit the induced loads (see section 3.1.1). Let (ij^) be 

inter 

the age of the tip vortex segment which first intersects the following blade, 

with the generating blade at azimuth angle Then a larger core radius r, 

b 

will be used if the line segment age is greater than 4), (\p) . The transition 

b 

at (j), occurs ini 
b 

inter* then propagates up the tip vortex at a rate 


tially a fixed increment after the intersection at 




{ 


CH'^ -4- 

+u CH- ~ aV — A.H* 


“ ( ‘•’-■►-LX- 


(only j = 1 to J need be tested if V, < 1, and cf, = 0 for all if 

b b 

V- > 1) . The initial, small tip vortex core radius r is used if .i < 4^i_ 
b ’ . c k+1 ^b 

the large core radius r. is used if (t), > ; and if 4), < the 

b k b k b k+1 

core radius is obtained by linear interpolation between r^ and r^^. The 
initial, small core radius r^ is always used for the near wake. 

The wake age 6. at the first blade/vortex intersection can be 

® inter 

determined by examining the projection of the tip vortex wake geometry 

(4^, 4^) on the disk plane. Consider a line segment extending from at 
w i 

4>, to r^ at 4^, , and the m-th blade at azimuth angle ^ + m27r/N. The 

vortex segment line is defined by 






C r. 
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with o “ 0 at and a = 1 at t 2 ; while the blade line is defined by 




with p = 0 at the hub and p = 1 at the tip. The intersection of these 

two lines is obtained by equating the and components of and 

V T T vortex 


blade ' 


r, . T 


I 7? — i ^ ^ X 

( . L.r - r,* Lt } — 


( fi, •■^1 _ 7^ .-JT > ^ 


which gives 


r. i-T j-r 


H»v- 


Lt ^ ‘■^T 






• <-T- ?r — r,. gy 

■ — CTJ — 7^”^ CAoVv" 


(There is no intersection if the denominator is zero) . 
tion of the vortex and rotor blade is 


The vertical separa- 


it. - i 


A blade/vortex intersection is defined to occur if 0 < a <1 and 

inter 

0*2 < P inter 1^1 ^ well. The intersections will be 

identified by examining each line segment in the tip vortex behind the blade 
at i)/, beginning with the wake age = Aij;. The segment will be tested 
against the blade index m such that the magnitude of 






- 17 ?- 



is iriinimized. If an intersection occurs, then 
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If an intersection does not occur, the next line segment is tested, 

3, 1,8 Rectangular vortex sheet. ^ Consider a planar rectangular vortex 
sheet element, as shown in figure 21. The induced velocity is required at an 
arbitrary point P in the flow field, defined by vectors from the four 
comers. The strength of the sheet is defined by the circulation values at 
the four corners and at and at and at the out- 

side edge, and at the inside edge). This vortex element is approxi- 
mated by a planar, rectangular sheet with sides s and t (fig. 21). The 
point P is defined by a vector r^ from the center of the sheet. The 

^ ^ A 

orientation of the sheet is defined by orthogonal unit vectors e and e 

► sc 

parallel to the sides of the sheet, and the normal unit vector e = e x e : 
^ n s t 


ir ^ 






A. 

es s 


Ju 

z 


( ‘■3 




1 


(t 


<^Z 


V 

<■1 ^ 


The vorticity strength is 6 in the 


A 

e. 


direction and y in the 


A 

e 


t ' s 

direction, varying linearly along the length of the vortex filaments. The 


minimum distance from P to the sheet or its extension is 


r : 
m 


^ ..3 A A 


A 


The vector r is perpendicular to the plane of the sheet and intersects it 
m 

at a point M. A coordinate system (a, t) will be used on the sheet plane, 
with origin at M so the center of the sheet is at 
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^ — So = 

T * = 


A. 


— • e.. 


The edges of the rectangle are then given by a = s ± s/2 and x 

o 

The distance from a point on the sheet to P is 


t ± t/2. 
o 


The linearly varying vorticity of the sheet is 




or in terms of the circulation at the comers: 

t - n-r> ^ ~ — *L 




'5 


_ — s 


r. - T 3 ^ 




Note that conservation of vorticity gives ^y/da = -96/8t or 

The Biot-Savart law gives the induced velocity of this vortex sheet: 


A V — — 
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Finally, the induced velocity in terms of the circulation at the four comers 
is : 

^TT 5't 

- tfs.- i'> i. + (S.- 3-3 _ £. 


£- <^5.+ I'^Jc X, - Cs.+ 4,'>t^Xj+r^ttl3 


where the first two terms are due to the trailed vorticity, and the last two 
terms are due to the shed vorticity* This may be written as 




( +^5^J + cr, ] 
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where 


V-ti = 



-ArrSt. 

^ "to ^ S ^“3 


''n 

\ 

irt 

r ^ ^ -t A 

L i Z ^ "^3 3 


it 

1> 

I 

S“tr 

r ^ ^ -r- , A 


-3 

\ 

A-wS-t 




There is a logarithmic singularity in the velocity induced at the vortex 
sheet side edges (the ^ singularity 

will be avoided by replacing the trailed or shed sheet vorticity by a line 
segment with a large core radius. Hence if 

^ ^ lt.«\ ^ -f cJvs 

then the point P is near a side edge of the trailed vorticity; if 

C~^0 ^ ^ civi ^ 

then it is near a side edge of the shed vorticity; and if 

\rl < civs 

for any comer, the point P is near the side edges of both the shed and 
trailed vorticity. 

An economical approximation is to replace the vortex sheet by a line 
segment, with either a linear or a stepped circulation distribution, and a 
large core size to eliminate the large induced velocity near the line. The 





strength and position of the line segments are determined from the circula- 
tion and position of the four comers of the sheet: for the trailed vorticity 




r -Cl 


= 

r, - 
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i 

: shed vorticity 




- v: 



r; 
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Note that with line segments it is not necessary to approximate the sheet 

geometry by a planar rectangle. The core radius can be specified arbitrarily; 

or r = s/2 can be used for the trailed vorticity and r = t/2 for the 
c c 

shed vorticity. 


3.2 Free Wake Geometry 

A method was developed in reference 16 for calculating the free wake 
geometry of a single rotor in steady state flight out of ground effect, which 
will be adapted for use in the present analysis. The wake model for the 
free wake calculation consists of line segments for the tip vortices; and 
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rectangular sheets or line segments for the inboard shed and trailed vake 
(similar to the far wake model used here for the induced velocity calculation; 
see section 3.1.2). The near wake or rolling up wake as described above are 
not considered. 

Only the distorted geometry of the tip vortices is calculated in the 

analysis of reference 16. The rigid or prescribed wake geometry is thus 

still used for the inboard vorticity. The distortion of the tip vortex 

geometry from the basic helix is defined in reference 16 as a vector D (^, 6), 

s 

giving the displacement of the wake element with current age 6 which was 
created when the blade was at azimuth angle i|^. A tip path pletne coordinate 
frame is used, with the x axis to the right (the advancing side), the y axis 
aft, and the z axis down. The procedure for calculating the wake geometry 
consists basically of integrating the Induced velocity at each wake element. 

The outer loop in the calculation is an iteration on the wake age 6. The 
induced velocity q('l') are calculated at all wake elements for a given age 
6, and all azimuth angles ij;. Then the increment in the distortion as the 
wake age Increases by Aif/ is : 

An efficient calculation of the wake geometry requires many variations in 
this basic procedure. Reference 16 adapted the near-wake and far-wake scheme 
for reducing the computation. The first time the induced velocity is evaluated 
at a point in the wake, the contributions from all wake elements must be 
found. For subsequent evaluations of the induced velocity at that point, 
only the induced velocity, due to the nearby wake elements are recalculated. 

The other major consideration for minimizing the computation is the matter 
of updating the induced velocity calculation. At a given point in the wake 
geometry calculation, there is a boundary in the wake between the distorted 
geometry and the initial, rigid geometry. The distortion has been calculated 
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between the rotor disk and the boundary; downstream of the boundary the wake 

is undistorted. As time increases by Av, the entire wake is convected 

downstream, and the rotor blades move forward by Aip, adding new trailed and 

shed vorticity at the beginning of the wake. If there were no distortion of 

the wake during the time Aijj, then the induced velocity at a given wake 

element would not change except for the contributions from the newly created 

wake vorticity Just behind the blade. Thus the normal calculation procedure 

consists of calculating the induced velocity at the boundary, by just adding 

at each step the contribution from the new wake directly behind the blade. 

Of course, the wake does distort as it is convected and as the estimate of 

the distortion improves, thus it is necessary to update the calculation of 

the induced velocity in the wake. In boundary updating, the induced velocity 

is calculated at the boundary still, by summing the contributions from all 

elements in the wake. In general updating, the induced velocity is recalcu- 

lated at all points in the wake upstream of the boundary. Boundary updating 

is typically done every 90° on the front and rear portions of the helices, 

and every 45° along the sides where the distortion is greater. General 

updating is typically done every 180°. General updating can not be done 

often if the amount of computation is to be kept low, but it does improve 

the accuracy and convergence. Numerous techniques for secondary improvements 

in the efficiency and accuracy were also included. The distorted wake 

geometry is required for m revolutions, where m decreases with forward 

speed. A single iteration of the free wake analysis consists of calculating 

the distortion D (i|^, 6) for = 0 to 27T and 6 * 0 to 27rm. Usually two 
s 

iterations are sufficient to obtain the converged solution for the wake 
geometry . 

The present analysis requires the wake geometry in the form of D(ij;, ((j) 
wliere is the current azimuth angle of the blade and cj> is the wake age 
+ 6 and 6 respectively in the notation of ref. 16). The present analysis 
uses a tip path plane coordinate system with x aft, y to the right, and z 
upward . Hence 
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is the distortion as used in section 3.1.3. The rotor velocity components 
V and V are required relative to the tip path plane. The hub motion 
and gust velocity at the hub can be Included in these velocity components. 
Ground effect and rotor /rotor interference can be accounted for by using 
an effective normal velocity: 

O - -t- Vi-T 

where f„„ is the ground effect factor defined in section 2.4.3, and 

ufc int 

is the interference induced velocity. The total uniform Induced velocity 

/'"h xf “♦"V 

is required to define the wake geometry at the start of the calculation. 
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4. AIRCRAFT MODEL 


4.1 Aircraft Configuration Definition 

A general two rotor aircraft is considered, with rigid body and elastic 
motion of the airframe. Aerodynamic forces on the wing/body, horizontal tail, 
and vertical tail are modelled, including aerodynamic control surfaces. The 
drive train connecting the rotors is modelled, with engine dynamics and a 
rotor speed governor. The case of a rotor or helicopter on a flexible support 
in a wind tunnel is included in the model. 


4.1.1 Aivovaft orientation: flight path and trim Euler angles,- The air- 

craft flight path is specified by the velocity magnitude V, and the orienta- 
tion of the velocity vector with respect to earth axes (fig. 22). The veloc- 
ity vector has a yaw angle (positive to the right) , and a pitch angle 

9_ (positive upward). Thus the climb and side velocities of the aircraft are 

r r 

V ^ = Vsin and V . , = Vsin cos0__. The aircraft attitude with 

climb FP side FP 

respect to earth axes is specified by the trim Euler angles, first pitch 6 

FT 

(positive nose up), then roll (f)__ (positive to the right). Airplane conven- 

r i 

tion is followed here for the coordinate systems — x positive forward, y 
positive to the right, and z positive downward (see reference 25). Between 
the earth axes (E system) and the velocity axes (V system) there are the rota- 
tions and Between the earth axes and the body axes (F system) there 

are the rotations G and 4^ . Thus the rotation matrix between the V system 
and the F system is: 


CSft 




- SeFTC«p(> 




FV 


^*F-r FT ^FP 

— C©pT 

— ^Fp 


~ ^»FT ^+Ft ^FF 




^•Fx^^F-r^afP ^fF 

Cftpp 

^©FT^FT^»FT ^ fP 
~ 

-+*C®pt^^Ft ^FP 
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The trim calculation determines the Euler angles 6™ and <{) (and perhaps the 

i i ^ r i 

flight path climb angle 0„ also). 

r r 

The velocity of the aircraft is V = so the components in the body 

axes are: 






The acceleration due to gravity is 


== gk or in body axes 






— ir^©FT- 


4,1.2 Rotor position and orientation,- The rotor hub position is speci- 
fied in the body axes relative to the aircraft center of gravity position, 

"r , = (xit, + 1 • The rotor orientation is defined by the rotation 

hub F F hub 

matrix between the shaft axes (S system) and the aircraft body axes (F system). 


- 185 - 





The position and orientation of the rotors relative to the body axes are fixed 
geometric parameters. The aircraft velocity is + = v3* 
so the shaft axes components of the velocity seen by the rotor^are 


I 











4f •'fv 



I 


The hub plane angle of attack and yaw angle may then be obtained from 
and the advancing tip Mach number from 

V < ' + yu| ‘ 


The sign of the lateral velocity must be changed for a clockwise rotating 

rotor; and for rotor #2 the velocity components must be multiplied by 
The quasistatic hub motion and the gust velocity at the rotor hub will be 
included in the advance ratio components: 





, J 


1 / 
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For a helicopter main rotor, the orientation with respect to the body 
axes will be specified by a shaft angle of attach 


(positive for rearward 




Le ^ (positive to the 

right) . 

— c® 


S© 






— 


The orientation of a tail rotor will be specified by a cant angle (j) 

R 

(positive upward), and a shaft angle of attack 6 (positive rearward). The 

K 

tail rotor thrust is to the right for a counterclockwise rotating main rotor, 

and to the left for clockwise rotation. Thus the definition of the tail rotor 

shaft axes depends on the main rotor rotation direction* Let have the 

mr 

value +1 for a counterclockwise rotating main rotor, and n = -1 for a 

mr 

clockwise rotation. Then the rotation matrix for the tail rotor is: 




r-c^ 


S^> S>e» 






L-Sa 

c© 



: and rotor 

of a tilting proprotor 

aircraft can be tilted by 


an angle 


where 


a = 0 for axial flow (airplane configuration), and 
r P 

= 90® for edgewise flow (helicopter configuration) . The rotor orientation 

is also described by a cant angle <J> (positive inward in helicopter mode, 

zero in airplane mode), and a pitch angle 0 (positive nose upward) which is 

K 

the angle of attack of the shaft with respect to the body axes when ot =0. 

P 

Thus the rotation matrix is : 

— 5 ^ 5 ^ 
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The rotor hub location r^^^^ for the tilting proprotor aircraft is defined 
by the pivot location r and the mast height h, so 

■m \7r\T ° * 



-3 




C ^ So< 

-C+S4 <1-0 


4.1,3 Wind tunnel case,^ For the case of a rotor or rotorcraft in a 
wind tunnel, the forces and moments on the body are reacted by the model sup- 
port system, so trim of the body forces is no longer a concern. The flight 
path and trim Euler angles can be set to zero (9__ = A =0 \i; =0') so 

the wind, earth, and body axes coincide (R^ = I). The wind axes and body 
axes are therefore the tunnel axes system, with the x-axis directed upstream, 
the y-axis to the right, and the z-axis vertically downward. The rotor 
orientation is specified by the matrix ^gp above. To accommodate the case 
of a wind tunnel with a turntable, the matrices can be post-multiplied 

by the matrix 






C9C4I 

Cj® Sh» 



Cj® 

0 





where is the turntable yaw angle, positive to the right, and 6^ is the 

test module pitch angle, positive rearward. 

4.1.4 Gust velocity The aerodynamic gust velocity will be defined 
relative to the velocity axes, with longitudinal component positive rear— 

ward, lateral component v^ positive from the right, and vertical component 
w^ positive upward = -u^l^ - v^jJ - w^k^) . The components in the body 

axes are then 


- 188 - 




The components in the rotor shaft axes are 



For a clockwise rotating rotor, the sign of is changed; for rotor //2 the 

gust velocities must be multiplied by Hence define 



including also the factor S^R^^/f^R^ and the sign change for a clockwise 
rotating rotor as appropriate. 

4.1.5 Aircraft description.- The aircraft geometrical description con- 
sists of the location of rotor //I, rotor //2, and wing/body, the horizontal 
tail, and the vertical tail relative to the center of gravity. The orienta'- 
tion and position of the aircraft components will be defined in a body axis 
system (the F frame) with origin at an arbitrary reference point, as in 
figure 23. Given the dimensional positions relative to the reference point, 
for example 


center of gravity: FSCG, WLCG, BLCG 

rotor //I: FSRl, WLRl, BLRl 
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Figure 23 • Bef'ini'tion of aircraft, geometry. 
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Then the coordinates of the location of rotor itl relative to the center of 
gravity (in the F frame) are 

X = (FSCG - FSR1)/R 
y = (BLRl - BLCG)/R 
z = (WLCG - WLRD/R 

and similarly for the rotor #2, wing/body, horizontal tail, and vertical tail. 

The mode shapes of the airframe elastic motion are described by the six 

components of linear and angular hub motion, in the F frame: £, and at r^ , 

k k hub 

(for each rotor) . Assuming that the generalized coordinate has dimensions 

of m or ft, it follows that the generalized mass has dimensions of kg 

or slug; that the hub linear motion 5^ is dimensionless; and that the hub 
angular motion has dimensions rad/m or rad/ft. These elastic vibration 

modes can be arbitrarily scaled; if and y^ are multiplied by a factor S, 

then should be multiplied by and the solution for qj^ will be 

divided by S. 

For the case of a wind tunnel with a turntable, the geometry will be 
defined for zero yaw angle, relative to a reference point at the center of 
the rotation. Then 



where R^^ is defined in section A. 1.3. 

4.1.6 Pilot controls.- The control variables Included in the rotor- 
craft model are collective and cyclic pitch of the two rotors, and the air- 
craft controls, which consist of engine throttle 6^, wing flaperon angle 6^, 

wing aileron angle 5 , elevator angle 6 , and rudder angle 6 . The control 

a e r 

vector is thus 
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The pilot's controls consist of collective stick 6^ (positive upward), 

lateral cyclic stick 6^ (positive to the right), longitudinal cyclic stick 

6 (positive forward), pedal 6 (positive yaw right), and the throttle 6 : 
s p t 

N/p — ^ 3 

For the purpose of trimming the helicopter, a linear relation between the 
pilot's control inputs and the rotor and aircraft control variables is assumed 

— T- i 

V = Vp -V- 


where v is the control input with all sticks centered (v = 0) , and 
o P 

is a transformation matrix defined by the control system geometry. 


T 

CFE 


The control transformation matrices for the single main-rotor and tail- 
rotor, the tandem main rotor, and the side-by-side main-rotor configurations 
are given below. The K’s are gain factors in the control system, and the 
Aij/’s are swashplate azimuth lead angles. The main rotor, front rotor, or 
right rotor is assumed to be rotor //I; and the tail rotor, rear rotor, or left 
rotor is rotor #2. The parameter Q here takes the value +1 for counter- 
clockwise rotation of the rotor, and = -1 for clockwise rotation. 
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4,2 Aircraft Analysis 

The aircraft motion consists of the six rigid body degrees of freedom 
and the elastic free vibration modes, A body axis coordinate frame with origin 
at the aircraft center of gravity (the F system) is used for the description 
of the motion. Airplane practice is followed for these axes — x is forward, 
y is to the right, and z is downward (ref. 25). The coordinate frame used 
is not a principal axis system however; moreover, the airplane practice of 
aligning the x-axis with the trim velocity is not followed, since for rotor- 
craft it is necessary to consider the hovering case (V = 0) . 

The parameters of rotor are used in making quantities dimensionless 
and in normalizing the aircraft equations of motion. It is assumed that 
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rotor //I is the main-rotor of a single main rotor and tail rotor helicopter; 
the front rotor of a tandem rotor helicopter; or the right rotor of a side-by- 
side rotor helicopter. With the hub forces in rotor coefficient form it is 
convenient to normalize the equations by dividing by the characteristic inertia 

A. 2.1 Degrees of freedom.^ The linear and angular rigid body motion of 
the aircraft is defined in the body axes (F system). The linear degrees of 
freedom are x (positive forward), y (positive to the right), and z 

r r r 

(positive downward) . These variables are dimensionless based on the rotor 

radius R; thus the velocity perturbations are normalized using the rotor tip 

speed QR rather than the forward speed V as is airplane practice. The 

angular degrees of freedom are the Euler angles (yaw to the right), 0„ 

r r 

(pitch nose up), and (p ^ (roll right). Then the linear and angular velocity 

r 

perturbations are 

Up ~ F ^ ^ y ^ 


where 






\ ^ 
O 




For the elastic motion of the aircraft in flight, the displacement u 
and rotation at an arbitrary point *r^ are expanded in a series of the 
orthogonal free vibration modes: 



are 


The first six degrees of freedom are the rigid body motions: ^32* ••^55 
respectively 6 

for K > 7 are the elastic modes of the aircraft. Orthogonality implies that 
the elastic vibration modes produce no net displacement of the aircraft center 
of gravity, or rotation of the principal axes. 

For the rigid body motions the mode shapes are simply 

Ct, ■••Ul =• [ x3 

C t ••• 1 = L Kt. o ;] 

since 

ja = o?p yr s -yp — 


p, ij)p, Xp, Yp, and Zp. The generalized coordinates 


4 . 2.2 Hub motion,^ The rotor equations of motion require the six com- 
ponents of the hub linear and angular motion in the shaft axis system: 





• • • 
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The matrix R«_ transforms the motion from body to shaft axes. The moment 
SF 

arm of the rotor hub about the aircraft center of gravity is in body axes, 
The total velocity of a point is the sum of the trim and perturbation 


velocities, u = V + Iq 


sk 


^ , in body axes. The rotor equations require the 
k 


velocity components at the hub in an inertial frame however (the S system) , 

and the Euler angle rotations between the body and inertial axes introduce 

perturbations of the trim velocity V, So the perturbation velocity becomes 

u = a„ X V + Zq , c, , where in the S system 
F SK K. 


Op V = 








\ /* ^*3 “*■ >** 3 / 


These contributions to the hub velocities (x^, y^, z^) cancel the terms in 
the blade velocity due to the Euler angle rotation of the inertial axes rela- 
tive to the trim velocity (the pa terms in u^, Up, and u^) . Thus the 

evaluation of the hub rotation (a , a , a ) for the aerodynamic analysis should 

X y z 

not include the body Euler angle contributions, as discussed in section 2.4 .2. 

Finally, the rotor hub acceleration is u = o)p x V + where the 

first term is the inertial acceleration due to the rotation of the trim 
velocity vector by the body axes angular velocity. This additional contribu- 
tion of the Euler angle velocity to the hub linear acceleration, in the shaft 
axes system, is 



which can be written i o = ^ with 


(L — 


o 
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For rotor //2 the linear hub displacement must be multiplied 

by account for the differences in normalization, c being based 

on rotor #1 parameters while a is based on rotor it2 parameters in this case. 

For a clockwise rotating rotor it is necessary to change the signs of y, , a , 

n X 

and These conversions will be included in the definition of c and c, by 

a, multiplying rows 1, 2, and 3 by ^^^^2 

b. changing the signs of rows 2, 4, and 6 

as appropriate. In addition, the derivatives of the hub motion of rotor //2 
must be corrected for the different time base, by multiplying the velocities 




. 2,^2 


by and the acceleration by 

5 ^. 


o< “ 


St: 




- = 


(See section 5.1.5 for the time scale correction in the case of harmonic body 
motion.) 

4.2.3 Pitch/mast-bending coupling.- Flexibility between the rotor swash- 
plate and hub will produce a blade pitch change due to elastic motion of the 
airframe. This coupling between the rotor pitch and mast bending is accounted 
for by introducing kinematic feedback of the following form: 

4.2.4 Equations of motion.- Following the usual steps of airplane flight 

dynamics analysis (see ref. 25), the linearized rigid body equations of motion 

are obtained by equating the angular and linear acceleration to the net moments 

and forces on the aircraft: lolt = IM and M(ut + w* ^ V) = I^. In terms of the 

r F F 

body axis degrees of freedom, including the gravitational forces, the equations 
are : 
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■ o 

r ^FT ■» 


Here M is the aircraft mass, including the rotors, and I is the moment of 
inertia matrix: 

I -I 

X xy 

I I 

xy y 

I -I 

xz yz 

= 0 if lateral symmetry is assumed) . These equations are dimension- 
less, and have been normalized by dividing by the characteristic inertia 

(i'iNI,),. Thus M* = /R^) and I* = l/(^NI, ) . Note that the gravitational 

b D b 2 

constant g is also dimensionless, based on the acceleration Q R; and 

M*g = ylC^/za^ where W = Mg is the gross weight. 




- 199 - 


f r 



For the elastic degrees of freedom, since orthogonal free vibration modes 
are used the equations of motion are simply 




V ( It 


where M, is the generalized mass including the rotors (in normalized form 
* ^ 2 

= M^/ (^NI^/R )), 0 )^ is the natural frequency of the mode, and g^ is the 
structural damping coefficient. 

* 

The generalized forces are due to the hub reactions of the two 

rotors, and the aerodynamic forces on the aircraft. Since the rotor mass is 
included in the aircraft inertia, the hub linear acceleration terras should 
not be included in the evaluation of the hub forces for these equations of 
motion. The aircraft aerodynamic forces are considered in section 4.2.6. 

Similarly, the rotor gravitational forces are not included in the rotor 
hub forces, since the rotor weight is included in the aircraft gross weight. 

4.2.5 Hub forces.- The aircraft generalized force due to the rotor hub 
reactions is 








M 




Normalizing by dividing by gives then 




t-vlii Js*?K 


^ 2Cw 

r-A 

V 2Ct 

4^' 
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or Q = c F, with c defined above for the hub motion. 

For rotor //2 it is also necessary to account for the differences in 
normalization, Q and c being based on rotor //I parameters while F is based 
on rotor //2 parameters. Thus the force coefficients of rotor #2 must be 
multiplied by 

c N Xu /rY, 

and the moment coefficients by 

<r “X V 'll) I 


For a clockwise rotating rotor it is also necessary to change the signs of 

C, , C^, and C, . Using c corrected for rotor //2 normalization and the rotor 
Y Q Mx 

rotation direction as required for the hub motion, it is then only necessary 
to multiply the matrix by 

T 

to obtain c for rotor //2. 

A. 2.6 Aircraft aerodynamic forces,- The aircraft aerodynamic forces 
considered are those acting on the wing/body (WB) , horizontal tail (HT) , and 
vertical tail (VT) . The generalized forces for the aircraft rigid body 
degrees of freedom, due to the aerodynamic forces and moments, are as follows. 
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Here L, D, and Y are respectively the aerodynamic lift, drag, and side 

forces; M , M , and M are the roll, pitch, and yaw moments on the wing/body; 
X y z 

and q is the dynamic pressure* The horizontal tail cant angles is 4^ _ 

ni 

(positive to left), and the vertical tail cant angle is (positive to 

right) . The moment arms of the aerodynamic centers of action about the air- 
craft center of gravity are in the body axes, r = xT^ + ytp + zT^ . The factor 
2Y/aa results from normalizing the equations by dividing by 
parameter A is the rotor disk area. 


The aircraft aerodynamic analysis thus requires the wing/body lift, 
drag, and pitch moment (L/q, D/q, and M^/q) as a function of angle of attack 
a and of the flaperon deflection 6^; and wing/body side force, roll moment, 
and yaw moment (Y/q, M^/q» ^ function of sideslip angle B and 

aileron deflection 6^; the horizontal tail lift and drag (L^^/q, 
a function of angle of attack and elevator deflection 6^; and the verti- 

cal tail lift and drag (L^^/q, Dy^/q) as a function of angle of attack a y^ 
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and rudder deflection 6 , These force and moment characteristics have 

r 

dimensions of length-squared and length-cubed respectively. 

The aircraft aerodynamic forces depend on the air velocity seen by the 
components and on the aircraft control positions. The air velocity consists 
of the trim aircraft velocity, the perturbation linear and angular rigid body 
contributions, the gust velocity and the rotor-induced aerodynamic interference 
velocity. In body axes (the F system) the total velocity is thus 



which must be evaluated at the wing/body, at the horizontal tail, and at the 
vertical tail. The angular velocity of the aircraft is 



The rate of change of angle of attack is also required (a = • The air- 

craft controls consist of flaperon, elevator, aileron, and rudder (6^, 6^, 6^, 

V- 

The aerodynamic interference velocity due to each rotor is required. 

With a nonuniform induced velocity calculation, X is the mean value of the 
wake velocity calculated at the position of the fixed aerodynamic surface 
(see section 3.1). The complete time history of the velocity, required to 
evaluate the mean, can be useful information itself. 

As a simple model for the aerodynamic interference, the rotor-induced 
velocities at the wing/body, horizontal tail, and vertical tail can be 
obtained as a linear combination of the mean induced velocity at the two 
rotors : 
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assuming that the induced velocity is normal to the disk plane (-k^ direction) . 
The K factors account for the maximum fraction of the aerodynamic surface 
which is affected by the wake, and the fraction of the fully developed wake 
velocity which is achieved. Typical values would be K = 1.5 to 1.8 (or zero 
for no interference) . The C multiplicative factors account for the decrease 
in the wake induced velocity away from the wake surface, using the following 
expression: 

Cl “ vvvococ C • ^ • 4- ^ 


where Z is the perpendicular distance from the aerodynamic surface to the 
nearest wake boundary (£ < 0 if the surface is inside the rotor wake cylinder) . 
Consider the geometry sketched in figure 24. The aerodynamic surface is lo- 
cated at (r - r^) relative to the rotor hub. The unit vector along the wake 
center-line is 


e ^ 




sf 




^•3 


+ V 


and we write 










(times R, /R^ for rotor //2) . Now the unit vector in the - ^ plane, 

perpendicular to e is 







Figure Zk , Model for rotor/airfrajne aerodynamic interference. 
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So the distance from P to the wake center-line is 


^ ^ 


= It X 


A point on the edge of the rotor disk is a unit vector perpendicular to k ; 


the point also in the e/? - r^ plane (perpendicular to "e - r„))is then 




X (" e K ("c I 


The distance of the wake edge to the wake center-line is thus 


= t- ^ - - e-ts 


I ^ 

^ t K(' r-<>^ 1 


So the distance from P to the wake boundary is 

^ — X(. -An ■= lexfr-<vt')l 1 - 


I t-T^ I 


= C I - 


It-fj 

six’- +«-*•’ 
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and j 

. - > 

C VMNA/>t ^ ^ ^ 


From the velocity components at the wing/body, the angle of attack and 
sideslip angle are: 




vsia 


= -G»~r' ^ 

u 


and the dynamic pressure is 




i 


VO 




The aircraft aerodynamic interference at the tail will be accounted for by an 
angle of attack change c and a sideslip angle a (positive when decreasing 
a and 3 at the tail) so the net velocity components are 


( 


U -h 
V — 
w — 


W ^ V \ 

u j 

lA / 


Then the horizontal tail angle of attack and dynamic pressure are 

gu*> 

U + ^ w -4- <T V 


o^H-r = 


ft V ” 


( ClA W 

Z- 


C' 


+(v-<ru^^-»- (w-~ € 

’ vsj 
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and similarly for the vertical tail 

T — > Cv — 

» V -I- — » 


u -+- € vv —4- r- V 




where ({5^^ and are the tail surface cant angles. The time-varying non- 

uniform inflow will increase the mean dynamic pressure in the wake: 

- 4 " a—"^ 


where 0 is the mean-square wake velocity perturbation, at the wing/body or 
tail as appropriate: 

cr- ^ ^ 


For best results, experimental data should be used for the aircraft aero- 
dynamic characteristics, including the airframe interference effects. As a 
simple model, the following expressions can be used for the wing/body: 


ts _ 


^V*03C ^ g «. • 
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and N 

are 

zr 

often negligible) ; 

UVT 
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L-VT 

— 




and for the horizontal and vertical 


i'v-r'^ -4- ^ 


To account for stall, a = sign a * min (lal, a ) is used in the wing/body 
lift and pitch moment, and in the tail forces. Here ^VT 

the zero lift angles relative to the reference body axis system; and 6^ is 
a wing flap angle. 

The wing/tail interference is evaluated from 
The area f - can be estimated from the wing area, span, and chord (S , S , c ) 

^ , WWW 

and the horizontal tail length by 

r , ^-725 . 2 .^ 

= 2-2 C5v,x/0 


(from ref. 26). A lag in the wake velocity at the tail is also included: 

9> 


Ae 




^ Vt 


2)c/ V 




(ref, 25). The wing-induced velocity could be obtained from the first order 
differential equation 
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•r ^ H- € 






The usual approach for airplane flight dynamics analysis is to write this 
equation as 




3 £ ^ 

\ -v*r i 


36 






Oof 


Using T = H„„/V for the time lag gives the above result. From reference 25, 
Hi 

the sideslip interference angle is approximately 



V 


where is the vertical tail position (positive upward). 

In summary, the aircraft aerodynamic forces are calculated as follows. 

The aerodynamic environment is defined by the helicopter trim velocity, 
perturbation linear and angular velocity, gust velocity, rotor-induced inter- 
ference velocity, and the aircraft controls. The total velocity components 
are calculated at the wing/body, horizontal tail, and vertical tall; from 
these the angle of attack and dynamic pressure are calculated. Then the aero- 
dynamic forces and moments on the aircraft are calculated. Finally, the 
generalized aerodynamic forces are evaluated. 

4.2.7 Aircraft aerodynamics — high frequency.- The aerodynamic model 
described in the preceding sections deals with the steady forces on the air- 
craft, and the stability derivatives for the rigid body motions involved in 
flight dynamics. Such a model would not be appropriate however for the high 
frequencies of rotor-induced vibration, for either rigid body or elastic 
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motions of the airframe. An accurate analytical evaluation of the generalized 
forces at high frequencies would require a sophisticated model of the wing/ 
body and tail aerodynamics, including the effects of the rotor wake-induced 
flow field, for the normal modes of elastic vibration of the airframe. Such 
an analysis is not attempted in the present investigation. 

The only generalized aerodynamic forces considered for the airframe 
elastic modes are the direct damping and control forces. In dimensional form, 
the equations of motion are then 












F^kqj^ and constants (with dimensions of length-squared and 


where 

length-squared per radian respectively) that depend only on the airframe 
characteristics; F, 




VV 


is the damping force divided by ^pV: 

/V 


and ^qj^^ control force derivative divided by 






'b6> 


The dimensionless form, normalized by dividing by then 


- (st) 
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4.3 Transmission and Engine Analysis 

The rotor rotational speed degree of freedom can be an important factor 
in the helicopter dynamics, and the rotor torque perturbations can produce 
significant drive train loads. A model is required which accounts for the 
coupling of the two rotors through the flexible drive-train, and for the 
engine damping and inertia. The drive train dynamics will be described by the 
rotor speed, the interconnect shaft torsion, and the engine shaft torsion 
degrees of freedom. The equations of motion are derived from the balance of 
rotor and engine torques (in the nonrotating frame) . A model for a governor 
with throttle or collective feedback of the rotor speed error is also con- 
sidered. 

4,3.1 Eng'iriB modeZ * -T\\e. engine model includes the inertia, damping, and 
control torques: 


The engine speed is and is the torque on the engine. The engine 

rotary inertia is is the engine speed damping coefficient, i.e., the 

torque per unit speed change at constant throttle setting: 




^ \ 


The variable 6^ is the engine throttle control position- is the torque 

applied due to a throttle change at constant speed: 



I 






Thus and can be obtained from data on the engine power as a function 
of throttle position and engine speed. 
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The engine damping may be related to the engine trim operating condition 


by 


k. 



where k is a constant depending on the engine type. This approximation is 
applicable to a wide variety of engines. The constant takes the value k = 1 
for a turboshaft engine (ref, 27) or for a series DC electric motor (ref. 28). 
It takes a value k = 1/(1 - r\) for an Induction electric motor or an armature 
controlled shunt DC electric motor (ref. 28; n is the motor efficiency). For 
a field controlled shunt DC motor, the only damping is mechanical and the 
damping of the load, so < = 0 (ref. 28). For a synchronous electric motor 
there is a spring on the rotational speed due to the motor, so the above model 
is not applicable (ref, 28) • Generally, the inertia of the engine or motor is 
more of a factor in the dynamics than the damping. 

The normalized engine damping and throttle coefficients are: 


'E 






- 


where r„ is the ratio of the engine speed to the rotor speed. When the 

E 

throttle control is only involved in the governor, the parameters required 

for Q is not really 8P /36 , but rather just over all loop gain of the 
t E t 

governor — the torque perturbation due to a rotor speed error. 

The transmission losses may be viewed as a viscous damping source, with 
a coefficient equal to 



wh6‘.re n is the transmission efficiency. This loss can be included in the 

2 ^ 

engine damping coefficient r_ by increasing the factor k by Ak = 1 - r, , 
In the equations for interconnect shaft torsion and engine shaft torsion, 
structural damping can be included as well. 

4.3.2 Equations of motion.- Figure 25 is a schematic of the transmission/ 
engine model considered for asymmetric drive train configurations, such as for 
a single main and tail rotor helicopter. The two rotors are connected by a 
shaft, and the engine is geared to one rotor (rotor #1 in fig. 25). The tor- 
sional flexibility of the drive train is represented by the rotor shaft springs 
and interconnect shaft spring K^, and the engine shaft spring 

The transmission gear ratios are r^ (the ratio of the engine speed to rotor 
speed), and r^^^ and (the ratio of the interconnect shaft speed to the 
rotor speeds) . Thus ^xi^^I2 ~ ^2^^1 ratio of the trim rotational 

speeds of rotor //2 and rotor #1. 

The degrees of freedom are the rotational speed perturbations of the two 

# • • 

rotors (VgX ^s2^ * engine speed perturbation The engine shaft 

azimuth perturbation is defined relative to rotor #1 rotation, so the 

total engine speed perturbation with respect to space is ^£(^gX ’ With 

the rotation of the two rotors coupled by the drive system, it is more appro- 
priate to use the following degrees of freedom: 

M’s 


Here is the differential azimuth perturbation between the two rotors. The 

degrees of freedom and therefore involve elastic torsion in the drive 
train. The degrees of freedom is the rotational speed perturbation of 

the drive system as a whole: both rotors, the engine, and the transmission. 

The differential equation of motion for the rotor speed d>Tnanics are 

obtained fror. equilibrium of the torques on the two rotors and the engine. 

The resultins equations for v* , v* , and v are as follows: 

s i e 






Figure 25* Schematic of rotorcraft transmission and 

engine dynamics model (asymmetric conf iguration) 
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where 0 , = ’.^ and >j; „ = (r^,/r_.) il/ + i(j-.; and with the engine by rotor #1 
Si s s2 II 12 s 1 

as in figure 25, the constants are: 
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The spring constants are normalized by dividing by (NI ) ; and I = 

Q I Cj 

When evaluating the derivatives of ^se in the analysis 

of rotor fr2, it is necessary to account for the difference in time scales: 









In these equations the gear ratio rotors is a posi- 

tive number, regardless of the rotor direction of rotation. Therefore here 
the sign of is not changed for a clockwise rotating rotor. 

With the engine by rotor //2, the constants in the equations of motion are 


-Nv*X. ^ — 


I ~ 


K/saE. ^ — 


c <!.« ^ k/v\^ C Kaai *4- ^X| 

Cr^,/C:^C^ ^6 1^£ I^Maa 

-4- ^ C ^X| ^>\aAa 


Figure 26 sketches a symmetric drive train configuration, as might be 
used for a side-by-side main rotor helicopter. The two rotors are connected 
by a cross-shaft, and there are two engines, one geared to each rotor. The 
degrees of freedom of this system are symmetric and antisymmetric rotor speed 
perturbations, and symmetric and antisymmetric engine speed perturbations. 
Dropping the antis>mmetric engine speed perturbation and antisymmetric throttle 
input, the system has three degrees of freedom as for the asymmetric configura- 
tions above. Let 
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Figure 26 . Schematic of rotorcraft transmission and 

engine dynamics model (symmetric configuration). 
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^Zt. 


again, where here r^^/rj.2 = !• The constants in the equations of motion are: 

1 

^ -K- 2 V^ai 

IC/vv ~H 2-*"^ 

-<- <£^ Vie 

Kj>^+ fe'V-E 

Here the engine inertia is for both engines, as are the damping and 

h 

throttle coefficients (Q^ and Q^) . The resulting symmetric and antisymmetric 
drive train motions can then be obtained from 

-f- 2 

= — t '♦'x 

as required. 
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With just one rotor (as in a wind tunnel) , the equations of motion for 
the rotor and engine speed perturbations reduce to: 

+ r|x^ + rt Oa ( Vi * 








f-t ( Vi H- >vt '+e, 

= Ce -V 


where here 




V-c 

IC/vk •4- 


Hence the equations for the asymmetric drive train configuration can be used, 
dropping the degree of freedom and the rotor #2 torque. 

The case of a rotorcraft in autorotation can be treated with this model 
by dropping the engine speed degree of freedom (4^^) , dropping the engine terms 
from the \|j and equations (helicopters usually have an over-running clutch 
to disconnect the rotors from the engine at zero torque) , and dropping the 
throttle governor control input (0^). The engine out case (engine and rotors 
still connected) requires dropping the engine damping term (or reducing it to 
just the transmission losses contribution) and dropping the throttle governor 
control input. The case of constant rotor speed is modelled by dropping the 
rotor and engine speed degrees of freedom and equations from the system. 
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4.3.3 Rotor speed governor.- When the rotor rotational speed perturbation 
is included in the dynamics analysis, it is usually necessary to also include 
the rotor speed governor for a consistent calculation of the rotor and aircraft 
behavior. The governor model considered is integral and proportional feedback 
of the rotor speed to the throttle, and to the collective pitch of rotors //I 
and y/2. The governor dynamics are represented by a second order lag. The 
control equations are thus: 








Note that c = is the rotor speed error, and ip is the integral of the 
s ° 

error. The integral gains are dimensionless (with 0 and both in radians 
or both in degrees) , and the proportional gains have units of seconds (Kp/K^ 
is the lead in the integral control) . When the throttle control is only used 
for this governor model, it is only necessary that the product of 3P^/3e^ 
and the governor gains be correct: 
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, or in terras of the 
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dimensionless parameters 
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The time constants in the governor equations can be alternatively described 
in terms of frequency and damping ratio: 



5. SOLUTION FOR THE ROTORCRAFT MOTION 

The solution of the equations of motion will be divided into two parts , 
based on the frequency content of the motion. The first part is the solution 
for the rotor motion and the airframe vibration. This motion is periodic 
with fundamental frequency n for the rotor, and Nfi for the airframe, hence 
it is high frequency motion. The second part is the solution for the steady 
state or slowly varying airframe motion, consisting of the aircraft rigid body 
and rotor speed perturbations. The assumption that this airframe motion 
occurs slowly relative to the rotor speed allows the solution of the equations 
of motion to be separated into these two parts. 

5.1 Rotor Motion and Airframe Vibration 

The equations of motion for the rotor and aircraft will be solved for the 
periodic motion by a harmonic analysis method, which obtains directly the 
harmonics of a Fourier series representation of the motion. After a con- 
verged solution for the blade motion and airframe vibration is obtained, the 
rotor performance is evaluated, including the mean aerodynamic hub reactions 
(in particular the rotor thrust and power) . The hub motion includes the 
static or quasistatic contributions from the aircraft rigid body motion. 

The helicopter state is determined by the control positions; the flight 
path angles and trim Euler angles (or test module pitch and yaw for the wind 
tunnel case) ; the quasisteady linear and angular velocity perturbations of 
the airframe; the quasisteady rotor speed perturbation; and the aerodjmamic 
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gust velocities. The rotor motion and airframe vibration are calculated for 
this state. Then the generalized forces acting on the rotor and airframe can 
be evaluated, as well as the various performance parameters of the aircraft. 

5.1.1 Fourier series representation.-- For the case of steady state 
flight, all the rotor blades execute the same periodic motion. It follows 
that the blade motion in the rotating frame can be written as a Fourier series: 





e 



z, 

- I 























bo 




where + mAi/j is the azimuth angle of the m-th blade (Aip = 2 tt/N, 

m 

m = 1 to N) , and is the dimensionless time variable. The complex and 

real Fourier coefficients are related by 6 = (3 - i3 )/2 and 

n nc ns 

e = - i6 )/2 for n ^ 1. The complex representation is most convenient 

n nc ns 

for solving the equations of motion, while the real representation is best for 

>(k) 


interpreting the motion. The notation 3 


n 


is used for the harmonics of the 


k-th bending mode. With the modes ordered according to natural frequency, 

is usually the fundamental lag mode and the fundamental flap mode. 

Similarly 0^ are the harmonics of the k~th torsion mode, with G rigid 

n n 

pitch and the remaining modes elastic motion of the blade. The Fourier 
representation of the gimbal or teeter motion is discussed in the next section. 

The degrees of freedom in the nonrotating frame are the aircraft rigid 
body and elastic motion, and the rotor speed perturbations. These degrees 
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of freedom are excited by the net rotor hub reactions, obtained by summing the 
root forces and moments from all N blades. Ideally, the rotor hub acts as a 
filter, transmitting to the nonrotating frame only those harmonics at multiple., 
of N/rev, The vibratory motion in the nonrotating frame is then also periodic, 
with fundamental frequency and can be written as a Fourier series: 





(K') 




for the body motion, and 




e. 


for the rotor azimuth perturbation (similarly for and ) . The static or 

1 e 

mean terms are obtained from the low frequency solution of the airframe 
equation. 

5.1.2 Gimbal and teeter motion.- The rotor gimbal motion (if present) 
is in the nonrotating frame, but it is most convenient to solve an equation in 
the rotating frame for the gimbal motion, along with the other rotor blade 
equations. From section 2.2.18, the gimbal equations of motion are given by 
equilibrium of the net longitudinal and lateral moments on the rotor hub: 

= o 

where is the flap moment at the blade root. All harmonics of the longi- 

tudinal and lateral hub moments cancel within the hub, except those at multi- 
ples of N/rev. The pN harmonics of the gimbal equations of motion are: 


-i£, 
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where (C ) is the n-th complex harmonic of C 
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Note that with the restriction that only the pN ± 1 harmonics of exist, 

this is equivalent to the relation 
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Substituting for the harmonics of and in terms of the harmonics of 

B„, the gimbal equations of motion give: 

G 

L *3 ” ^ 

or 

— _j ■+' i 

■ 4 - "X^ ( ^ +-Z (^?6^pf0+l 
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which we note are just the pN ± 1 harmonics of the equation 




(except, for the effects of unsymmetric gimbal springs or dampers* and the 
fact that the damping is in the nonrotating frame) . 

The equation for the teeter motion of a two-bladed rotor is 

~ 2 ^ I ^ ^ f-r ^ O 

All the even harmonics of the root flap moments cancel within the rotor hub. 
The summation operator only transmits the odd harmonics to produce teeter 
motion. Hence the solution for the teeter motion can be obtained by solving 
the equation 

(ix =0 

for the odd harmonics (i.e., for the pN ± 1 harmonics, just as for the 
gimbal motion) . 

Thus the gimbal or teeter motion can be obtained by solving the rigid 
flapping equation in the rotating frame: 

for Just the pN ± 1 harmonics of g . Then the harmonics of the gimbal 
motion are 
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from which the gimbal motion can be evaluated 



and so 


The harmonics of the teetering motion are 

for the odd harmonics, from which 

f»6i = (^T = e. 


5.1.3 pitc^ motion.- In the limit of infinite control stiffness, 

the equation of motion for the blade rigid pitch degree of freedom reduces to 
MpAg “ 0 (see section 2.2.9), with the solution 
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^ Pr 5= C'V -V (V*+- 

— ^ \^\ %\, — VC^^ ^<i 

-\- z^©c^vc- “4“ 

Write the total root pitch motion p as the sum of p and the motion due 

*^o r 

to elastic distortion of the control system: Pq ” ^r ^d* Substituting for 

p^, the rotor equations of motion will be solved for the harmonics of p^. 

The case of infinite control system stiffness then requires only that the 
equations for p^ be dropped from the solution procedure. Writing 
p^ = p^ + Pj introduces terms due to p^, p^, and p^ in the equations of 
motion. 

Allowance will be made for different stiffnesses in the collective and 
cyclic control systems by using different natural frequencies for the collec- 
tive motion cyclic motion (0^^^), and the reactionless motion 

n « 2) . Given the collective, cyclic, and reactionless control stiffnesses, 
the dimensionless natural frequencies are obtained from 



or the natural frequencies can be specified directly. 

The control system damping will be specified for the collective motion, 

the cyclic motion, and for the rotating frame motion. The total damping is 

then Ca + Co ^ for the collective mode, Co + Co ^ for the cyclic 
^coll ^rot ^cyc '^rot 

modes, and Cn ^ for the reactionless modes (n i 2). 
t>rot 

5,1.4 Harmonio analysis solution,- A harmonic analysis method will be 
used to integrate the differential equations of motion, solving directly for 
the harmonics of the motion. Consider equations of the form 
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where 6 Is the degree of freedom, K and M are the appropriate stiffness 
and mass, and g is the forcing function (usually nonlinear). To avoid the 
singularity of the resonant response at harmonics near the natural frequency, 
it is necessary to include the damping terms on the left-hand-side of this 
equation. Thus the term C3 is added to both sides, giving 

-t-Kp = V 


where C is the doping coefficient. For good convergence the damping co- 
efficient used should be close to the actual damping of the particular degree 
of freedom, including structural, mechanical, and aerodynamic damping sources. 
The damping estimate does not have to be exact however, since it is added to 
both sides of the equation. In fact the actual damping in the forcing function 
g will often be time varying and even nonlinear, so the viscous damping 
coefficient has to be an approximation. The sole function of this damping 
term is to avoid divergence of the solution near resonance, and the value of 
C has no influence on the final converged solution. 

Now the function F is evaluated at J points around the rotor azimuth; 


where (j = 1 to J and = 2ir/ J) . Then the harmonics of a complex 

Fourier series representation of F are 
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With = 1 these harmonics would give a Fourier interpolation representation 

of While it matches the function exactly at the known points F(ij^,) 

(or with least squared-error if the number of harmonics used is less than 

(J-l)/2), the Fourier interpolation gives a poor representation elsewhere, 

with large excursions due to the higher harmonics. In particular, poor 

estimates of the derivatives of the function F are obtained. With the above 

values for (which reduce the magnitude of the higher hatinonics) , and an 

infinite number of harmonics (n * - «> to <») , a linear interpolation between 

the known points F(i|jj) is obtained. By truncating the Fourier series 

(n = -L to L) the representation of F Is smoothed, the comers of the 

linear interpolation being rounded off. Usually L = J/3 is satisfactory, 

that is the number of azimuth stations should be about 3 times the maximum 

harmonic of interest. The azimuth step thus should be Ail/ * 2tt/J ^ 120/n 

max 

degrees. Then the solution of the equation of motion for the harmonics of B 
is obtained from the harmonics of F by 



where 


H ® K - Mn + Cin. 
n 


The iterative solution, required because the nonlinear forcing function 
F depends on 6 and 6, proceeds as follows. At a given azimuth , the 
blade motion is calculated using the current estimates of the harmonics : 


p = Ac 

M .f-c ‘ 


The forcing function F^ is evaluated next. The estimates of the flapping 

harmonics are then updated to account for the difference between the current 

value of F. and that found in the last revolution; 

J 
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After adding AB to the flap harmonics 3 , the azimuth angle is incremented 
n n 

to y The calculation proceeds around the azimuth in this fashion until 

the solution converges. The test for convergence is performed once each 

revolution. Requiring that the root-mean-squared change in the blade motion 

from one revolution to the next be below a specified tolerance, the criterion 

is (A6) < e for all degrees of freedom, where 

rms 

and 

This test is applied to all degrees of freedom, for both the rotor and the 
airframe. 

In the present problem, the system of equations and degrees of freedom 
can be separated into two sets: the rotor motion, consisting of flap/lag 

bending, rigid pitch/elastic torsion, and gimbal or teeter flapping; and the 
aircraft rigid body and elastic motion, with the engine/transmission as an 
independent subset. The coupling between these sets is accounted for in the 
nonlinear forcing functions. As long as the coupling is weak, it is possible 
to solve the two sets of equations separately, in parallel. Within each of 
these subsystems, it is necessary to solve all the equations simultaneously. 
Including in particular the inertial coupling on the left-hand side. Thus a 
vector equation must be solved for each harmonic of the motion. The solution 
proceeds as follows. At a given azimuth station, the blade motion and hub 
motion are evaluated using the current estimates of the harmonics. Then the 
generalized forces are evaluated, and the rotor equations are solved to update 
the harmonics of the blade motion. Next the rotor hub reactions are evaluated 
(for which the updated blade motion harmonics can be used), and the aircraft 
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equations solved to update the harmonics of the body motion. The azimuth 
angle is then incremented, and the calculation repeated until a converged 
solution is obtained. 

5.1.5 Motion evaluation.-' To begin the solution at a new azimuth station, 
the deflection, velocity, and acceleration of each degree of freedom must be 
evaluated from the harmonics. For the rotor blade bending: 






t•^ 





- * 

= r ?I!L «•»'*> 

v.s> 


Similarly for the blade pitch/ torsion, p, , p, , and p are obtained from the 
(k) K K K 

harmonics • . Note that for rotor #2 these time derivatives are based on 

n 

For the gimbal/teeter motion, recall from section above 

St - SI 

\ V\ =. ^K> ^ 1 

and similarly for and 

Recalling that only the pN harmonics are excited in the nonrotating 
frame, the rigid body and elastic airframe motion of the aircraft is 


— 

s: 4>!T 

vt=^»o ^ 


Vt = 

- . V 

vv. e-. 

v» = 



VI 



- 233 - 



where are the six rigid body degrees of freedom, and q for 

sk 

k _ 7 are the airframe elastic modes. The steady state or slowly varying 

rigid body motion contributes the static velocity terms (q ) (k < 6* 

’sk static - ’ 

this motion is static compared to the high frequency rotor motion and airframe 

vibration). The static elastic airframe deflection gives (q ,) (k > 7). 

sk static 

The rotor hub motion is then 



where c is given in section A. 2. 2. Recall that in the evaluation of a 

x’ 

a^, and (for the aerodynamic analysis) the contributions of the rigid body 
Euler angles (}ip, 6^ and < 1 )^ (q^^ to q^^) are not included; also the linear 
hub displacements (x^, y^, and z^^) are not used in the rotor analysis. 

Hence a is evaluated due to the elastic airframe modes only (q , , k > 7). 

S K 

The velocity and acceleration of the hub are 

For rotor #2, a and a are multiplied by and respectively. 

Also for rotor #2, the aircraft motion harmonics are at n = pN 
(relative to the time scale of the nonrotating frame, fi^) . We can write the 
hub velocity as follows therefore: 
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+1"’ - ] 


— ^ 

z - 

= ^ veN §* +^ e. '' J 

f <* ' IVfNN 
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"] 


Hence by evaluating the hub motion as a sum of harmonics at n * pN, with the 
azimuth angle of rotor #2, the time scale will be automatically accounted for. 
Similarly for the acceleration 




ro-»oc 2. 


=• C 


i - f C(^y <C } 

-4- ^ C ^ ^ ipfo I 


A factor of Q^/Q^ is still required in the second term of a, to account 
for the scaling of the aircraft velocity in c. 

The acceleration due to gravity, considered as an equivalent linear 
acceleration, is 




2 2 

For rotor //2, the factor (fi R)^/(fi R)^ is also required. 

The rotor shaft angular motion perturbation is evaluated from 


=. 
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Mis p|0 
^ m 
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6-v**^ ^ 




Vk a 

From section 4.3.2, the angular motion perturbations of the two rotors 
then obtained from 




= 

Vs 


As for the airframe motion, the time scales of the velocity and acceleration 
perturbations for rotor #2 are accounted for without additional factors of 

The solution for the helicopter trim or transient motion (section 5.3) 
supplies the static motion of the airframe rigid body degrees of freedom 
(^sk k = 1 to 6) and of the rotational speed degree of freedom (iji ) . 

The solution for the static elastic deflection of the airframe and drive train 
(Qsk k - 7, and is given in section 5.1.10. 

The dimensionless time variable is i|», the azimuth angle of rotor #1. 

The azimuth angle of rotor #2 is 



M> 

fz. 




21 


where is the angle when tp = 0 at rotor #1. Hence the analysis of 

rotor #2 must account for this phase difference of the two rotors, by evaluat- 
ing the blade motion, the airloads, and the hub reactions of rotor #2 at 
<P + ^^21’ phase of the time variable for the airframe motion is the same 

as that for rotor #1. Note that if 02 /^! ^ 1 the rotors do not maintain a 
fixed azimuthal phase difference. For that case the rotors will effectively 
be analyzed separately however, so ^'^21 ~ ^ used. 
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Hence the harmonics of the rotor hub motion are obtained from the har- 
monics of the aircraft degrees of freedom by the following expressions: 





— VI 




a nonzero multiple of 

N; and the "static" 

components are 









• • 

^ 1 


r- 


+• 

L 0 


for the displacement (a^ and summation is over the elastic air- 

frame modes only; also, only the angular displacement components are required 

2 2 

(a^, Uy, rotor #2, a factor of R) R )2 is required in the 

gravity term, static multiplied by matrix c is 

multiplied by Also, for rotor #2 the harmonics are multiplied by 

e. 


since for the azimuth angle of rotor #2 equal to the hub motion at 
^ ^^21 required (only if * ^ 2 ) • harmonics of the rotor azimuth 

perturbation are obtained from the harmonics of the drive train degrees of 
freedom by the following expressions: 
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for n a nonzero multiple 

of 

N; and the ’’static” components are 


*s=r O 


— 







CVs 

For rotor #2. (4'g)gtatic 

has 

been multiplied by also for rotor #2 


the harmonics are multiplied by 


Lv% 


The blade rigid pitch motion requires the pitch increment due to the 

governor and due to the rotor mast bending. The harmonics of the governor 
pitch increment are obtained by summing the contributions from the two rotors; 
and for rotor //2 multiplying by 
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The pitch increment due to mast bending is 

* 4 - 
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So the harmonics are obtained from the harmonics of the airframe elastic 
motion by the following expression: 


^^WMcrtT^ c=. 
yvi 


.(k). 



(O 
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4» 

^n— 1 

fcs-T 
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no 

when 

n is a 


kMC|^ — I ^Mlvt 


KmC|l 


n-fi 


(k) 


vention ^ (q ) is used here. For rotor #2, the harmonics 4> 

o sx statxc n 

are multiplied by 

— Uw 

e. 


in this calculation. 

5.1.6 Rotor equations^-- The differential equations of motion for the 
rotor degrees of freedom are given in section 2.2.18. For the n-th harmonic, 
these equations take the form 


- 239 - 




where is the n-th hamonlc of p,, and B_ is only present for the 

n ' a t>n 

n “ pN ± 1 harmonics for gimballed or teetering rotors. From the equations 
of section 2.2.18, the transfer functions matrix and the forcing function 

F are as follows. 
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where 










C2. w 

~Y 

Cv,-0 


C VA -4- 


v\. 

yjULXsUT' 


Note that estimates of the aerodynamic spring and damping forces have been 
added to both sides of the equations. In the matrix these terms must be 

multiplied by ( (J/n7r)sin(niT/J) ) to be consistent with the Fourier analysis 
of the forcing function F (see section 5.1.4). 
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The solution for the blade motion requires the Inverse of H for each 

n 

harmonic. One approach Is to Invert every azimuth step. A more 

efficient approach Is to Invert H once and store the result. It will still 

-1 “ 

be necessary to update occasionally however, because it depends on the 

bending solution (q . terms). The blade motion harmonics should be completely 

^ -1 

recalculated whenever is updated. 

Hence at each azimuth step <1;^ the forcing function Fj is evaluated. 
Then the blade motion harmonics are updated by adding the following increment: 




(P; - ^ 



5.1.7 Hub reactions.- The generalized aircraft forces due 

T 

hub reactions are Q « c F, where c is given in section above 


to the rotor 
and 
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From the results of section above, the required hub reactions are as follows. 
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where 
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These hub reactions are harmonically analyzed in the evaluation of the 
vibratory airframe motion. The effect of the summation over all N blades 
then is to suppress the harmonics not at multiples of N/rev. An equivalent 
approach therefore is to omit the summation operator and only evaluate the 
pN/rev harmonics. Hence the aerodynamic forces are only evaluated for one 
blade at azimuth angle v/;. 

The mean hub reactions, required in the solution for the steady state 
or slowly varying aircraft motion (section 5.3), are obtained by averaging 
the above results over one period. Note that only the aerodynamic terms 
contribute to the mean values of the hub forces and torque. 

The hub forces due to the linear acceleration of the rotor mass 



have not been included here, since the airframe inertias include the rotor 
mass. For the static elastic airframe deflection, evaluated from the 
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mean hub reactions, the ^’static” hub acceleration must be included: 




which consists of the gravitational force, and the centrifugal force due to 
the angular velocity of the body axes. 

To solve the equations of motion for the engine and drive train, the 
rotor torque is required in the following form: ' 



tJ 




5., 

- ciiT-'-ftC 5., 

— ^ C *^4 ©S 


(see section 5.1.9). 

5.1.8 Aircraft equations.- The differential equations of motion for the 
aircraft degrees of freedom are given in section 4.2.4. For the n-th harmonics 
of the rigid body motion, these equations take the form 
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The forcing function is due to the rotor hub reactions 


f 



rwVfcr- 





From the equations of section 4,2,4, the transfer matrix is as follows. 





O 


— y\^ 


For the n-th harmonics of the aircraft elastic motion, the equation of motion 
is 


-2US- 



where Qj^ is the generalized force due to the rotor hub reactions (k > 7). 

The only aircraft aerodynamic forces included in the high frequency response 
are the damping coefficients of the elastic modes (see section 4.2.7). 

Only the harmonics of N/rev are excited in the nonrotating frame. The 
aircraft response to each of the two rotors is evaluated separately. For the 
response to the rotor #1 hub reactions, the harmonics at n = 2N^, 

etc., are required. The time scale of the aircraft equations is the rotational 
speed of rotor #1, so the harmonics of motion due to rotor #2 are at 
n = For the response to rotor #2 hub reactions the harmonics at 

n = ( 02 / 02 ^^ 2 ’ required. The equations of 

the airframe motion are not solved here for the static response (n * 0) . The 
hub reactions of rotor #2 are evaluated as a function of its azimuth angle, 
i|» 2 ; to obtain the response of the airframe these hub reactions must be used 
at »(» = i |)2 ~ ^’^21' 

The hub reactions are evaluated at azimuth stations \I> as the rotor 

j 

equations are being solved. Then the airframe vibration motion is obtained 
from 


for n 


a nonzero multiple of 

“3 



z: 


N; where 

e. 

» i 



are the harmonics of the rotor hub reaction. The motion is evaluated 

n 

for excitation from rotor #1 and for excitation from rotor #2. For rotor //2 
the harmonics must be multiplied by 
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are 


since for the airframe motion at the hub reactions at + ^^^21 

required (only if • 

5.1.9 Transmission and engine equations.- The differential equations of 
motion for the rotor and engine speed perturbations are given in section 4.3.2. 
The equations for the n-th harmonics are 



with the forcing function and transfer function matrix as follows. 
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where 


Sr. = 
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The torque is defined in section 5.1.7. In summing over all N blades, 

all the harmonics of the torque cancel except those at N/rev. Hence the 
drive train equations are only solved for the n = N, 2N, 3N, etc., harmonics. 

The rotor torque (y ?p/aa) is evaluated at the azimuth stations 
as the rotor equations are being solved. Then the transmission vibratory 
motion is obtained from 



are the harmonics of the torque. For rotor #2 the harmonics are multiplied by 

e. 

The transmission motion is evaluated for excitation from rotor Itl and for 
excitation from rotor //2* 
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5«1«10 Statio eZastia defZeatioK,- The equations of motion for the static 
elastic deflection of the airframe and drive train are 

m: + II X [V‘ ( 


- VC-MX. ( ^ -+ 


(from sections 4.2.4 and 4.3.2). Here Qj^ is the mean generalized force 
due to the two rotors, and is the mean rotor torque. 

Hence the solution for the static elastic airframe motion (k > 7) is 




V^-Vic 


JSSS- . t r ^ ^ #1 
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T T 

where c, is the k-th row of c ; and 


f, - 3 S, ^ 




is the mean hub reaction including the rotor mass inertial reaction. The 
static elastic drive train motion is 
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and 


yCg/oa 


is the mean torque. 


5.1.11 Tb)o~-rotoT aircraft, In the present model, the two rotors of a 
helicopter can influence each other through excitation of vibratory airframe 
motion. The analysis proceeds as follows (see fig. 27). The rotor analysis 
calculates the hub reactions of each of the two rotors. From these hub re* 
actions, the aircraft equations of motion are solved for the harmonics of the 
airframe rigid body and elastic motion. Then the hub motion can be evaluated 
at each rotor, due to the forces of each rotor. The motion at each hub due to 
the two rotors is summed. Then the rotor equations are solved for the rotor 
motion and for the hub reactions again^ 


It is useful to be able to suppress the feedback of the nonrotating frame 
vibration to either or both rotors. The coupling can be suppressed by 
omitting the summation of the two hub motion components at one or both hubs 
(the dotted line in fig. 27). The entire vibratory hub motion can be suppressed 
by setting it to zero at one or both hubs (the static or low frequency hub 
motion and the acceleration due to gravity should be retained however) . 
Suppressing the entire vibratory hub motion is equivalent to dropping the 
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Figure 2?. Outline of dynamic interaction of the two rotors 
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where 
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and 


YCp/aa 


is the mean torque. 


5.1.11 Tb)o-rotor aircraft,-- In the present model, the two rotors of a 
helicopter can influence each other through excitation of vibratory airframe 
motion. The analysis proceeds as follows (see fig, 27). The rotor analysis 
calculates the hub reactions of each of the two rotors. From these hub re- 
actions, the aircraft equations of motion are solved for the harmonics of the 
airframe rigid body and elastic motion. Then the hub motion can be evaluated 
at each rotor, due to the forces of each rotor. The motion at each hub due to 
the two rotors is summed. Then the rotor equations are solved for the rotor 
motion and for the hub reactions again. 


It is useful to be able to suppress the feedback of the nonrotating frame 
vibration to either or both rotors. The coupling can be suppressed by 
omitting the summation of the two hub motion components at one or both hubs 
(the dotted line in fig. 27). The entire vibratory hub motion can be suppressed 
by setting it to zero at one or both hubs (the static or low frequency hub 
motion and the acceleration due to gravity should be retained however) . 
Suppressing the entire vibratory hub motion is equivalent to dropping the 
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aircraft degrees of freedom as far as the rotor analysis is concerned, but it 
still may be useful to evaluate the aircraft vibration due to the hub reactions. 

It should also be possible in a similar fashion to suppress the hub motion 
due to the static elastic deflection of the airframe and drive train. 


The procedure described above is based on the assumption that the entire 
system is periodic, which in fact is true only if both rotors have the same 
rotational speed “ 1) • When the two rotors do not turn at the same 

speed, the motion in the nonrotating frame is not periodic even in steady 
flight. The most important example is the single main'-rotor and the tail-rotor 
configuration. In order to analyze a periodic system still, it is necessary 
to neglect the mutual interference of the two rotors. The analysis proceeds 
as described above, except that the hub motion of one rotor due to the vibratory 
airframe motion produced by the other rotor is always suppressed (the dotted 
line in fig. 27) . Effectively the helicopter is then being analyzed as two 
single rotor systems, except for the coupling through the aircraft ’’static” 
motion (steady state of slowly varying, including the airframe static elastic 
deflection) . 


5.1.12 Circulation convergence.^ The blade motion will be calculated 
for a given induced velocity distribution over the rotor disk, uniform or 
nonuniform. When the converged solution for the blade motion is obtained, 
the rotor loading (C^ or bound circulation) is re-evaluated. Then the induced 
velocity estimate can be updated, and the blade motion solution repeated. 

This procedure continues until the root-mean-squared change in the bound 
circulation from one iteration to the next is less than a specified tolerance 
level: 



S. ^ ^ -s: ) 

U 


where r is the maximum bound circulation and the summation is over the 

j 

azimuth. When uniform inflow is used, the criterion is 
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To Improve the convergence of the iterative calculation of the rotor 
loading and wake induced velocity, a lag is introduced in the thrust coefficient 
and used to calculate 


where C^Q2.d thrust used to calculate X^ in the last iteration, and 

^Tnew thrust calculated using that value of X^. The factor f should 

have a value equal to the thrust lift deficiency function 

1 

^ I _• .21a 

* A a*T- 


(see section 6.1.5). Similarly, for a nonuniform inflow calculation a lag is 
introduced in the blade bound circulation used to evaluate the induced velocity. 

5.1.13 Calculation procedure.- In suninary, the calculation of the rotor 
motion and airframe vibration proceeds as follows. The input quantities are 
the linear and angular velocity perturbations of the rigid body motion; the 
rotational speed perturbation; the collective and cyclic pitch control angles 
of the two rotors; the aircraft aerodynamic control positions; and the gust 
velocity components at the two rotors. The output quantities are the general- 
ized forces due to the mean hub reactions of the two rotors; and the converged 
solution for the blade motion and aircraft vibration. 
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The outermost loop is an iteration on the rotor induced velocity and 
bound circulation evaluation. The next loop is an iteration on the rotor and 
aircraft motion calculation. 

One cycle of the blade and aircraft motion calculation consists of the 
following steps. First the transfer matrices are evaluated. Then there 

are a number of cycles of successive evaluation of the rotor and airframe 
motion by the following procedure. First the hub motion harmonics are evalu- 
ated. Next there is an azimuth loop for the rotor. At each azimuth step the 
rotor blade motion harmonics and the aerodynamic hub reactions are updated. 
Next the total hub reactions are evaluated; the aircraft vibration and drive 
train vibration harmonics are updated from the hub reactions; and the static 
elastic deflection is evaluated. 

Within the azimuth loop of the rotor motion calculation there are the 
following steps. First the hub motion and blade motion are evaluated from 
the harmonics. Then there is an integration over the radial station. Within 
the radial station loop the blade section pitch, velocity, angle-of-attack, 
and Mach number are evaluated; the lift, drag, and moment coefficients are 
evaluated; and the section aerodynamic forces are evaluated. The generalized 
aerodynamic forces of the blade modes are evaluated by integrating the section 
forces over the blade radius, and then the blade motion harmonics are updated. 
The aerodynamic hub forces and moments are also updated. 

After each cycle of the blade motion calculation, the convergence is 
tested by comparing the blade and airframe motion harmonics with the values 
at the beginning of the cycle. 

Finally, after the converged blade motion is obtained, the induced 
velocity and circulation convergence is tested by comparing the rotor bound 
circulation with the values at the beginning of the iteration. 

5.2 Rotor Performance, Loads, and Noise 

Once the solution for the periodic motion of the helicopter has been 
obtained, the performance, loads, and noise of the rotor can be evaluated. 

The rotor loads of interest include the tension and shear forces, bending 
moments, and torsion moment at various blade radial stations; the control 
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loads; the blade root forces and moments; and the net rotor hub reactions. 

The rotor- induced vibration can be evaluated at various points in the aircraft. 
From the rotor aerodynamic loading the rotational noise can be calculated. The 
rotor loads at radial station r will be calculated by integrating the aero- 
dynamic and inertial forces acting on the blade outboard of r. 


^ j-‘ 


5.2.1 RotoT perfomanae.~ To evaluate the rotor performance, the mean 
rotor hub reactions are required: 

-4- ca. 4H»: -Sr \ 
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where the simnnatlon operator averages the forces and moments over the azimuth. 
These quantities are directly available from the rotor analysis (see section 
5.1.7), where they are used also to calculate the generalized forces acting 
on the aircraft. 

The rotor performance is determined in particular by the thrust and torque. 
The power delivered to the rotor through the shaft is P * fiQ. The propulsive 
force is the component of the net rotor force in the direction of the aircraft 
velocity: 
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and then the rotor lift, normal to the aircraft velocity, is 


L 




—1 


V"- 
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7 . 


So L and X are the wind axis components of the net aerodynamic force of the 
rotor. 

The hub reactions relative to the tip-path plane are obtained by multiply- 
ing the force vector by the rotation matrix 





Pi 


where 6 and 3 are the tip-path plane tilt angles. Also of interest are 
c s 

the magnitude of the net force of the rotor, and its tilt angles relative to 
the reference plane: = tan ^ (H/T) and <\> = tan*^(Y/T). 


It is useful to split the rotor power according to the type of energy 
loss. The induced, interference, profile, parasite, and climb power losses 
are obtained as follows. 


*= -VX/^AC-s^y C.M/Sc^'^ 
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= <=-^ — C<=ff Cf^> 

<=-<>. = C C4-. -4- o. -► Cf^"> — 


The Induced and interference power losses are obtained by integrating the in- 
flow velocities over the rotor disk. 

The ideal power loss, consisting of the parasite, climb, and minimum 
induced losses, is defined as follows: 


where 


^ C&m 




Z J A* 


Then the nonideal power loss, consisting of the profile and excess induced 
losses, is defined as 

«= CCf. -KZf, C.T- 

A measure of the rotor efficiency is the figure of merit, defined for hover as 








{ 


This will be generalized for axial and forward flight to 
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AA = 


\ — 






The ratios ^Pint'^^T considered equivalent induced velocities, 

A measure of the excess induced losses is the ratio 


K = 






and the profile losses can be expressed in terms of an equivalent section drag 
coefficient: 


* * 




= 's 


In forward flight, the rotor drag is defined by 


br = ^ H-X ~ 




V 


or in terms of an equivalent drag area f = • A measure of the rotor 

efficiency in forward flight 
the total drag is defined as 


efficiency in forward flight is the rotor lift to drag ratio, L/D^. Similarly 


1 _ 

V 


and the total lift to drag ratio Is ^/^total* 
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The section power loading can be split into components in a similar 


fashion: 




<4 /•- 


V- ^ 






-H -p 


— 7 U 


5.2.2 Section force.- The total tension and shear forces on the blade 
section at radial station r are F(r) = F* - where 

A X 

== Sr 


SV -•- (-fit ') 


(see section 2.2.10; here a Includes the gravitational acceleration). Hence 


=. ■» 
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Substituting for 

O VSA. 

r 




a from section 2.2.4, the inertial term becomes 

.1 


— Lft 


“-Ha 



- 4 - 2 - 4 - 

-f- <c zt4- 




- s!-'~^ iu 

(C.5,+24^'»wH' 


i 



or 




-»- “♦“Ys — 

■+ 4 ^*^* 


■*"^a 


— C ^ cj^v -I- C^ 

H- ‘2- C -+- Vs 

+ ^ 2 .-^ . 


H- s 



^ %% 




- sC 


The components of the tension and shear forces in the local blade axes are 
obtained by multiplying ^gjjg^j. rotation matrix 
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where 



i% (s^ ■+* WA| — Wa-j^ ^ 

-t- ( 


The force is multiplied by 
principal axes, where 


to obtain the components in the section 




& 




O ^<^0 


O I 

O 



-I 


The Integrals of the aerodynamic forces will be evaluated as follows 




r*4- 




i»X4-l 


T;-** 
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i*X 

where I is the maximum index (over i = 1 to M) such that 

rj. - ijAr^ < r (I «= 1 if r^ - I 5 A r^^ i r) ; and hj^ = Ar^ for i = I + 1 to M, 

except 

X.X * <■*.■** t < r 

— Zsk^X. VV-S^ ^ <"x. — % 3 ^ ^ 


-O -I- 




The inertia constants are evaluated by trapezoidal integration over the blade 
properties defined at r = (j - 1) Ar (j = 1 to M + 1, Ar = 1/M): 


Sr^''0’^ -= S<-*' 

■+- z 

X. 

cr 



X- r> +2^*- - 

* ■ 2 A.r 

+ V*. L - 

2. 

M 

HT ^ Ar- 

• A a ”* • 

1 — a:H-l 

a. 


t 

= s: !Uv 

L = :t~i 


where I is the index (i = 2 to M + 1) such that _ 1 - ^ 
h^ “ Ar, except 
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2_Ar 


C<~x — C ^ 2-Ar 

^ * 

2. 


— 0 AC 




Similarly 


x-2. 




Ar 


= 7* Z 


X-Z 

1- c 


c ^ 

^X.-I 2 


r-*!- 


Ar 


1= ' 


H- Ox-. + «iif!^;:^2^*-'’J-i^ 
^(nViivv 


with I as defined above; and h^ = Ar, except 


JL ji “ 


(t- — Cj^ -4-Ar'^*’ 

21«- 

^r* 0*x- 

^ + — — 


-Ar^ 


Jt. JU, — ±1 

-^1 2. 


5.2.3 Section bending moment.- The bending moment on the blade section 
at radial station r is 
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where M = M 1 + M k = (ii „ + kk )M. As in section 2.2.7, the aero- 

* A 2^0 wCS 

dynamic and inertial moments are as follows. 


g" . t c a - - fe ^ 





V 


- S'f ^ ’f 


Substituting for the acceleration from section 2.2.4, the inertial moment 
becomes: 
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+ ir^E^" - <-s* ■> <"4 cci, 
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The components of the bending moment in the principal axes of the section are 
then obtained by 


M,\ 

►'i Is 


C-OS0 


J \ 4 


The integrals are evaluated as described in section 5.2.2. 

5.2.4 Section torsion moment.- The torsion moment on the blade section 
at radial station r is 


M 


-Wri 


Mr 





AVa. 
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where . if - + (,xt + zjk)') * "m . As in section 2.2.8, the aero- 

dynamic and inertial moments are as follows. 


Mr, 






=. JL 


— >0?y \rj 

^ ^ Vr ^ ^ 

+ Sr 
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Then substituting for the acceleration from section 2. 2. A gives for the 
inertial term: 




E f; Sc 

^ Sc ^ ^ 

*' I 

4- f ^i. S*c 

IM C»fA ^CA V«Aj 

< S€A»-AfA»^ 3 

+ '^s. Y|/ (ri •T? ^ — v|; v»» 

X r* ^ [^: --^l: (cO ^ 

r\^-€ -Xoi V Irl •‘4 IPi • ^ 
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where 


X = 


— (^- '■')(-?,-e- X y 


or 




=> ? ^N>?» f*' 

'' %0 

+ =4^* • <" Tf» ( ?o -I- p&'t — li Vj 

+ f ^ y. 

4- * ( TCf -+- ^ c-*o4>^'^ 

-4- 

4- 


The Integrals are evaluated as described in section 5.2.2. 
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5.2.5 Control load,- The control load at the blade root is given by the 
moment about the feathering axis: 






ir 





where Mp^ = e^^ • M. From section 2.2.9, the equation of motion for rigid 
pitch is 

Hence the control load can be evaluated from 

to -+- C N ^ p: 

L 

-+ T‘^.c*+ 

-h (st, -+- ^ 

— — Jr 


The inertial constants are defined in section 2.2.19. 
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5.2.6 Root forces and moments.- Following section 2.2.18, the forces 
and moments at the blade root are 


Coetl 


'v«ot 



= V ( Srtk -4- 

0 

If 

k Ti ^ -iv. 1- 

(1 



^■4 - »■ A-r — ^ *• 


Substituting for the acceleration from section 2.2.4 gives the components of 
the blade root force and moment in the rotating frame: 






^ ('6 


- st Vi 

•4- c. -"^4 C \x, — 





— aC ( -4- ^ 

+ 2-St ^'^2 

-4- S 2 -t, ii. 
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The inertial constants are defined in section 2.2.19. Here the gravitational 
forces is included in the linear acceleration terms, and the root reaction 
forces due to the blade mass and weight are retained. 

To be consistent with the Fourier analysis used in the equations of 
motion, the aerodynamic forces must be operated on as follows. Let F be 
the aerodynamic force evaluated at = jA(|; (j = 1 to J, Ai(i = 2ti/J). The 
function F is harmonically analyzed as in section 2.1. A, and then the 
function is re-evaluated (at = kAi)/) using those harmonics: 


^ VI 

t. 
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f, 
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where 





This operation is applied to the aerodynamic load used to evaluate the section 
force, the section bending moment, the section torsion moment, the control 
load, and the root forces and moments. 


5.2.7 Hub veactions • - Following section 2.2.18, the total force and 
moment acting on the rotor hub, resolved in the nonrotating frame (the S system) 
are 



The hub reactions are obtained from the root forces and moments by resolving 
the rotating components into the nonrotating frame and summing over all the 
blades: . - . . rr . 
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Applying the summation operator gives directly the hub reactions in the time 
domain, over a period 2 tt/N. Alternatively, the rotating frame forces and 
moments can be harmonically analyzed, and then the harmonics of the rotor hub 
reactions obtained from the harmonics of the root forces and moments, by 


1 



fc. ' 



Or 

V»mp|sl 




k. ^ ^ I 

r~^i 







Only harmonics at multiples of N/rev are transmitted by the hub to the non- 
rotating frame. It is useful however to evaluate all harmonics of the root 
forces and moments in the nonrotating frame, according to the above relations, 
since a real rotor will not accomplish this filtering exactly. 

5.2.8 Vibration.- Following the evaluation of the hub motion in the 
section 4.2.2, the linear acceleration in the aircraft at an arbitrary point 
r is 

- 4 - {^\ 

Yr ' 

(inertial acceleration in the F coordinate frame) . The harmonics of the 
vibration are thus 
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»y 

Dividing by the dimensionless acceleration due to gravity (g/fi k) gives the 
vibration in g^s. 

5.2.9 Fatigue damage assessment.- The fatigue damage due to the blade 
bending loads ^ torsion loads, and control loads is determined principally by 
the mean and 1/2 peak-to-peak values. An improved estimate of the fatigue 
damage can be obtained by using Miner's rule, with the method of reference 29 
for counting the loading cycles. 

The counting method applied to periodic loading consists of the following 
steps. First all the relative minima and maxima in one revolution are 
identified. The absolute maximum and absolute minimum give one loading cycle; 


So 







The absolute maximum and absolute minimum are discarded then, leaving a set 
of L peaks. 

Consider the first group of K peaks (K ® 3 or 4 usually, unless there 
is a lot of high frequency variation in the loading) . The 1/2 peak-to-peak 
value within this subset of K peaks gives the amplitude of one loading cycle 
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Consider the next group of K peaks, using the last peak of the previous 
group as the first peak of this group, to identify the loading cycle amplitude 
82 * The counting procedure continues in this fashion, taking K peaks at a 
time to identify loading cycles. 

Each group uses K-1 new peaks to obtain one loading cycle, so in K-1 
revolutions of L peaks each, L cycles will be identified. This means that 
each loading cycle identified occurs 1/(K-1) times per revolution. For 
periodic data, it is equivalent to consider the L groups of K peaks, 
starting at each of the peaks in the set over one revolution. 

Then for one revolution of the rotor, there has been identified the load- 
ing cycle 

Z 


occurring one time per rev; and the loading cycles 


each occurring 1/(K-1) times per rev. 

Miner's rule for the damage fraction is 



where n is the nxunber of applied cycles at level S, and N is the allowable 
number of cycles at level S. The S-N curve will be approximated by 


N 
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where is the endurance limit; and C and M are constants depending on 

the material. Hence 




c S 





Using the loading cycles amplitude and frequency determined above, gives the 
damage fraction 










for one revolution. 


The damage fraction is required for an analysis of an actual rotor design, 
but for more general investigations a parameter emphasizing the applied load- 
ing would be useful. Often in fact the loading will be below the endurance 
limit of the specific design considered, but still an assessment of the in- 
fluence of the loading waveform and amplitude on the relative fatigue damage 
is required. Such an assessment can be obtained by considering the damage 
fraction obtained from Miner’s rule for small endurance limit. In that case, 
the damage is proportional to 













M 


A more useful representative of the fatigue characteristic of the rotor load- 
ing is 


= C 


which is the equivalent 1/2 peak-to-peak loading amplitude which would alone 
give the same damage number as the actual loading. 
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5* 2*10 Rotatvonat noise,- The rotational noise due to the blade aero- 
dynamic loads and thickness will be calculated using the following equations 
for the far field harmonics of the sound pressure disturbance. Consider a 
rotor moving through the air with velocity components 


Using a tip-path plane coordinate frame, these velocity components are 

N/ ^ k 


An observer position relative to the hub, moving with the rotor, is defined 
by the components x^, y^, and (positive aft, to the right, and upward); 
or in terms of a range S^, an elevation angle 0^ (positive above the tip- 
path plane) , and an azimuth angle (defined as for the rotor azimuth) ; 

— So 
a So 


Write the blade section 

forces 

as 





= 
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where F^TPP’ ^xTPP* ^rTPP section aerodynamic forces; and 

^x’ chordwide loading distribution functions. The section forces 

relative to the tip-path plane are obtained from the forces relative to the 
shaft axes as follows: 







^ coi*V •+ oos'i*— 

ooiH> -4-^S 


This loading can be written as Fourier series: 
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e. 


The blade thickness distribution is written as 

±(^') =1 


where is the section area, and a(x) is a chordwise distribution func- 
tion. The velocity normal to the section, *= J^r + also be 

expanded as a Fourier series 
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= s^r 
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*N=-» 






Sl*r 


where . i, - (V^ - iV^)/(2fir), and V_^ - (V^ + iV^)/(2«r). 

The sound pressure at the observer has a periodic component, which is 
the rotor rotational noise: 







for t “ 0 to 2Tr/Nfi. Then the harmonics of the far field rotational noise due 
to the rotor blade lift, drag, and thickness are: 
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where the argument of the Bessel function 
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the sound pressure harmonics can be presented in dB: 


SPL^ = \o JUi^ 




for a single sided spectrum “ *00002 N/m ). The overall sound pressure 


level is 


The time history is also of interest for rotational noise. 


OSPL = 


The chordwise distribution functions for lift and thickness give 
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A simple and generally conservative approximation is impulsive lift and thick- 
ness, for which ~ ^n ~ airfoil type loading 





gives 


J? 


vy 


— ?=■ 

e 






where x is the midchord distance from the spanwise axis of the blade. For 

2 

NACA 4- and 5-digit airfoil series, A = .685tc where t is the thickness 

xs 

to chord ratio. For these airfoils 


O Y\ 




CO^ ^ 




where 


O.C = [-zyiAjf - -«i.o^ _ .1CH. 


A rectangular loading distribution gives 

=- 


2r . Sn 

^ 2-<” 


which might be appropriate for the drag. 


5.3 Steady State or Slowly Varying Aircraft Motion 


The linearized equations of motion for the rotorcraft rigid body degrees 
of freedom are 









' 5 . 










"dc - 4 - 



H- M 




from section 4.2.4. Here 

-kT 

g = first 3 elements of 6th row of c 

f = first 3 elements of 6th column of 

T 

with the matrices c and c defined in section 4.2. The terms involving 

•k k 

and account for the rotational moments of inertia of the two rotors. 

The rotor mass is included in the aircraft gross weight and moments of inertia. 
Note that M*g = y2Cy/aa. Hence if these equations are divided by 2y/a (rotor 
#1 parameters) , they will be in the form of rotor coefficient to solidity 
ratio, with the components in the body eixes (F system). When the rotor speed 
perturbation is included, the equations of motion become 
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where Includes the rotational moments of inertia of the two rotors: 





The solution for the generalized forces due to the airframe aerodynamics is 
described in section 4.2.6. The solution for the generalized forces due to 
the two rotors is described in section 4.2.5. 

In the trim analysis, these equations are solved for the case of steady 
state flight. The controls are adjusted until the desired operating condition 
is achieved. 

In the transient analysis, these equations are numerically integrated in 
time. A non-equilibrium flight path is produced by a prescribed control or 
gust input. 

In the flight dynamics analysis, the stability derivatives in a linear 
expansion of the rotor and airframe aerodynamic generalized forces are obtained 
by prescribed perturbations of the body motion and controls. 

5.3.1 Trim anatysis.- The helicopter trim calculation determines the 
control positions and aircraft orientation for the specified flight condition. 
For steady state flight the perturbation rigid body motion is zero, so the 
net force and moment on the aircraft must be zero. Thus the rigid body 
equations of motion give six equations to be solved for the six trim variables, 
consisting of the four pilot’s controls (6^, 6^, 6^, and 6^) and the two trim 
Euler angles (6pj. and • The controls are adjusted until equilibrium 
flight is achieved for the specified flight condition. For level flight 
(9pp = 0) or a specified climb velocity, it is assumed that the engine can 
supply whatever power is required to maintain the rotor rotational sneed. 
Alternatively, the aircraft power available can be specified (such as for 
power-off descent). Then there is an additional trim variable (the flight 
path angle ®pp) and an additional equation to be solved (the power required 
equals the power available) . 
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The helicopter can also be trimmed in a steady turn by prescribing the 

turn rate A coordinated turn is obtained if zero sideslip C'/'™ = 0) is 

“ F F 

specified. In forward flight the resultant bank angle should be 


where n is the normal load factor. 


Hence the trim analysis solves the equations of motion for a specified 
steady flight speed and rotor speed (and possible a specified turn rate) . 
Setting the perturbation rigid body motion and rotational speed to zero gives 
the following equations: 




( 






CVWr Z. 






^ ^ rv H>p 





c> 






f 5 sN 





The contributions to the force and moment are from the hub reactions of the 
two rotors, the airframe aerodynamics, the acceleration due to the turn rate, 
and the aircraft weight. 

The pilot’s controls are collective stick 6^, lateral cyclic stick 6^, 

longitudinal cyclic stick 6 , pedal 6 , and throttle 6 . The controls 

s p r 

of the two rotors and aircraft are related to the pilot’s controls by 
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V = 


c-pe. 


where is the control input with all sticks centered (see section 4.1.6). 

The throttle control variables (6^. and 6^) are not used for the trim analysis. 

For some rotor craft configurations the pilot's collective stick (6 ) does 

o 

not control the rotor collective pitch, rotor trim being handled by a rotor 

speed governor using collective pitch feedback. In such a case the static 

component of the blade pitch governor, (A0 ) ^ for either or both 

govr static 

rotors, can be used as the trim variable in place of 6 . Hence A6 is 

o . govr 

add^d to the rotor collective pitch 9-^ obtained from v and v . 

/5 p o 

The table below summarizes the options considered for the trim analysis. 

For each case there are a number of trim variables, which are adjusted to 
achieve the target values of an equal number of trimmed quantities. In the 
free flight cases, the helicopter is trimmed to force and moment equilibrium. 

In the wind tunnel cases rotor //I is trimmed to a prescribed operating condi- 
tion. The trim option and the degrees of freedom representing the aircraft 
can be specified independently; hence it is possible to use a free flight trim 
option with an analysis of a helicopter in a wind tunnel. The options called 
wind tunnel cases are however more typical of wind tunnel test configurations, 
particularly with only one rotor rather thm the complete aircraft. The trim 
variables consist of the four pilot’s controls, aircraft orientation parameters, 
and wind tunnel orientation parameters. The aircraft orientation parameters 
consist of the trim Euler angles, flight path angles, and turn rate; they are 
used only for the free flight cases. The wind tunnel orientation parameters 
consist of the test module yaw and pitch angles. 

The free flight cases include the following options. In the level flight 
case the pilots controls and the aircraft Euler angles are used to trim the 
six components of the net force and moment on the aircraft to zero. In the 
climb or descent case, the flight path angle is used in addition to trim the 
power required to a specified value. (In vertical flight e^p = ±90® however, 
so the parameter which would have to be varied to achieve the specified power 
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Helicopter Trim Options 


Free Flight Cases 
No trim 

F, M 

F, M 

F, M, CpA- 
F, M, Cp/tr 



Wind Tunnel Cases 
No trim 

Cp/®^ 

?<=• K 

Pc’ Ps 

Cx/v-, Cy/o- 
^c’ 

V— . CxV 

Cj^/*“ . Cx/«»“ 

ClA- , yo- . 


t 4. 


trim variables 

^Ft'*FP ^p"^F ''^T 



T “ trim variable 
Z = zero 

aircraft velocity and rotor 
rotational speed fixed 


- 297 - 


requirements is the helicopter vertical speed.) Optionally the sideslip angle 
ij^pp can replace the roll angle as a trim variable. A useful alternative 

is to trim the longitudinal variables only. The net vertical and horizontal 
forces and pitch moment are trimmed to zero using collective and longitudinal 
cyclic stick controls and aircraft pitch attitude. This longitudinal trim is 
exact for the case of a laterally symmetric aircraft in a symmetric flight 
condljtlon. It is also a useful approximation which may converge better than 
a full six degree of freedom trim analysis; and when neglecting the tail rotor 
in the analysis. 

The wind tunnel cases include the following options. The rotor thrust or 

power can be trimmed with collective pitch. The rotor tip-path plane tilt can 

be trimmed with lateral and longitudinal cyclic pitch. The tip-path plane is 

defined as the first harmonics of the tip deflection z , ; 

tip 

€ • 


The rotor lift and drag (in wind axes) and the side force can be trimmed using 
collective and cyclic control. Either the drag coefficient to solidity ratio 
C /a can be specified, or the equivalent drag area X/q (so Cjc = (X/q)i5 
(V/nR) /Aj^). As an alternative, the shaft angle of attack can be used in 
place of longitudinal cyclic or collective pitch as a control variable. It is 
also possible to trim only the longitudinal variables. 


The trim iteration can also be omitted. In this case the helicopter or 
rotor performance is evaluated for a specified control setting. 


In the free flight cases, the criterion for convergence of the trim 
iteration is that the net force and moment be less than a certain tolerance 
level as specified by the parameter e : 


\^C H- (Cr^lr-y- ] * 

CZ^I<r 


<€ 
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C -+■ -f- ^ 

.o^ Cv*,^ 


and 

1 <^/o- •— C ^ ^’^-hxx-yrt 1 

C ^ • e*c» I 


< ^ 


when the power is triinmed as well. For the wind tunnel cases, the following 
criterion is used: 


\s - ^ 








where f = C^/a, Cp/a, C^/cf, Cj^/a, or (X^/a as appropriate (with “ -Ol* 

.001, .01, .001, or .001 respectively). The criterion for the flapping is 


\ (ic - 1 < C 

- CPs')^,^ \ < € 


(with g and 3 in radians), 
c s 

The trim problem is to find the values of the control variables such 
that the target values of certain trim quantities (lb are achieved. The 
complex, nonlinear equations involved require an iterative solution procedure. 
A first order expansion of gives 




NV 






= M„ 


bv 


or 
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V„+, =- v„ + fc. 


The partial derivative matrix required is 

. i/A r ^/A ^ 

h = 57 = [ ■•• dT; ■•• J - 


[■•■ 




.V£, 


1 


The factor F i 1 is Included to avoid overshoot oscillations in the trim 
iteration by reducing the step size. 

—1 

The matrix D is constructed as follows. First M , , is calculated 

^ old 

using the initial values of v. Then each control v^ is decreased by the 

increment Av_, , and M is calculated; then the i-th column of D is 
i new 

given by 









and the control is restored to its initial value. Finally the matrix D is 
inverted, and all elements multiplied by the factor F. The partial derivative 
matrix can be recalculated occasionally as the iteration proceeds, to improve 
the convergence. Generally a step size of about A = 1 degree is satisfactory 
for all control variables. 

5.3.2 Transient analysis*- The helicopter transient analysis involves 
numerically integrating the equations of motion for the rigid body and rotor 
speed degrees of freedom. A non-equilibrium flight path is produced by a 
prescribed control or gust input as a function of time. The assumptions of 
this analysis are that the aircraft motion is slow compared to the rotor 
speed, and that the perturbed rigid body motion is small. The assumption of 
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quaslstatlc body motion allows the periodic rotor motion solution to be used 
with the transient analysis. The small motion assumption arises because it 
was assumed that the perturbation rotor hub motion is small; it is consistent 
therefore to integrate the linearized equations for the rigid body motion. 

The degrees of freedom considered are the six linear and angular rigid 

body motions (x„, y„, z„, 6„, and the rotor speed perturbation (i ) . 

r r r r r r ® 

The input parameters are the aircraft controls (6 , 6 , 6 , 6 , 6 ) and aero- 
dynamic gusts. Optionally any of these degrees of freedom can be held constant. 
The transient analysis with all seven degrees of freedom fixed produces the 
rotor response to control and gust inputs. The equations to be integrated are: 



where AQ is the rotor and aerodynamic generalized force, less the trim 
value; and AC^ is the rotor torque, less the trim value. The initial 
conditions are zero (except for when the helicopter is trimmed in a 

steady turn) . 
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The transient rotor speed perturbations will produce throttle and rotor 
collective pitch increments due to the governor: 

and 

for rotor #1 and rotor //2. The rotor azimuth perturbation will also 

produce cyclic pitch Increments due to the trim swash-plate tilt: 



For rotor #2 these cyclic pitch Increments must be multiplied by An 

autopilot is also included, since the transient rigid body motions can be 
divergent in some flight conditions: 

At, — It, ( ’c, B» ep') 

Hence the pilot *s control positions consist of the trim setting, the transient 
term, and the autopilot term; and the Individual control positions are obtained 
from 


with the governor contributions added to the elements of v as required. 
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The transient calculation begins with the trim solution, at time t = 0. 
The pilot’s controls or gust input are specified as a function of t. The 
equations of motion above are in the form 

A fourth order Runge-Kutta method will be used to numerically integrate these 

equations from time t to time t . , » t + h: 

n n+ 1 n 

It, = > V > , 

'3 

V.I 4 = i > >Jo •*''^3') 

~ '^■o V' "to ( R| “1" 

^„+. = V ^ ( It, + zk» + zKj + fc,,') 

Note that it is necessary to solve for the periodic rotor motion four times 
per integration step, 

5.3.3 Flight dynamics analysis.- The flight dynamics analysis here con- 
sists of a calculation of the helicopter stability derivatives and an analysis 
of the resulting linear differential equations. As for the transient analysis 
it is assumed that the body motion will be slow (compared to the rotor rota- 
tional speed), so the quasistatic rotor solution can be used. The assumption 
that the perturbation body motion has small magnitude is here consistent with 
the stability derivative representation of the rotor. 

The equations of motion are the same as considered for the transient 
analysis (section 5.3.2). Here the rotor hub forces and the aircraft aero- 
dynamic forces are expanded in terms of the stability derivatives. By making 
successive perturbations to the Inputs for the rotor and aircraft analysis, 
the generalized forces can be expanded as follows; 
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The coefficients of the matrices are the aircraft stability derivatives. 

There are contributions from each rotor, and from the wing/body, horizontal 
tall, and vertical tail. The airframe terms Include the effects of the rotor- 
induced interference velocities. The result for is similar (rotor 

contributions only) . The gust velocity here is uniform throughout space . 

The rotor speed governor, defined by the following control laws 

Vs — V4 

= Kf. Vs 

will be directly included in the stability derivatives. It is Also necessary 
to account for the cyclic pitch perturbations due to the rotor azimuth 
perturbation: 
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For rotor //2 these cyclic pitch increments must be multiplied by 

The result is a set of linear differential equations of the form 






describing the flight dynamics of the aircraft. The state vector x consists 
of the rigid body and rotor speed degrees of freedom; 

X [ 4* % V y -e* Vj 


The control and gust vectors are defined as follows 


— ^ ^-t 3 




— [ 


W# 


3 
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(the gust components are in wind axes). Optionally any of the degrees of 
freedom can be omitted from the analysis. Using these equations the helicopter 
flying qualities can be examined, in terms of the eigenvalues and eigenvectors. 
The transfer function (pole-zero set or frequency response) or the transient 
response to a prescribed control or gust input can also be obtained. The 
transient response can be calculated by numerically integrating the equations 

^ ~ ■— AjA,h -t- 


in a manner similar to the transient analysis (see section 5.3.2). In this 
case only a gust that is uniform throughout space can be considered. 

The stability derivatives are obtained in body axes (the F system) rela- 
tive to the aircraft center of gravity. There are contributions from each of 
the rotors, and the aircraft aerodynamic components (wing/body, horizontal 

tall, and vertical tail) . The following notation is used for the stability 
derivatives . 


Equation 

Notation 

Variable 

Subscript 

roll moment 

* 

^x^ 


P 

pitch moment 

A 

I M 

y 

•f 

q 

yaw moment 

A 

I N 
z 


r 

longitudinal force 

M*X 


u 

lateral force 

M*Y 

^F 

v 

vertical force 

M*Z 

^F 

w 

torque 

A 

ij, 

n 

aircraft inertias are 

introduced so 

that the coefficient 

of the highest 

1 derivative in an equation is unity. 

The derivatives are 

defined with 


positive signs on the right-hand side of the equation of motion. 
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Application of this procedure to the wind tunnel case will give the wind 
or shaft axis derivatives of a single rotor alone. 

5.3.4 Transient gust and control.- Transient gust and control inputs are 
required when the equations of motion are Integrated for the transient or 
flight dynamics analysis. The pilot's control increment is required at 

time t. A simple form is 


where v is a constant vector and C is a scaler function of time. More 
po 

generally, the control input can be a function of the aircraft motion as well. 

T 

The gust velocity g = (u v_ w„) (in wind axes) is required at the 

O 0 G 

time t, at the location of both rotors and the airframe aerodynamic components. 

Consider a convected gust field, defined by a function G(x , y , z ). The 

S S S 

gust coordinates have origin at the center of gravity when t = 0. The air- 
craft velocity is V , in the x (wind axis) direction. The gust field has 

a 

velocity V in the negative x direction; the gust is coming from azimuth 
S S 

angle ip relative to the aircraft (see fig. 28). Hence given the position 
r (in wind axes) , the location in the gust field is 





'Hi 

~ 


1 




The position vector for the wind /body is 
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with In the F frame, relative to the aircraft center of gravity; the 
position vectors for the horizontal tail, vertical tail, and rotor hub are 
obtained in a similar manner. The positions on the rotor disk are 


(neglecting the tilt of the tip-path plane relative to the hub; the sign of 

the ^ component is changed .for a clockwise rotating rotor). A one- 

s 

dimensional convected gust field is defined by 




V 


Note that the gust is convected at the rate V relative to the aircraft if 

g 

V is not used, and at a rate V relative to the fixed frame if V is 
a g a 

used. With V ■ 0 the gust field is stationary (relative to the aircraft 

o 

or the earth if V is not or is used) . 
a 

Alternatively, a uniform gust field can be used, which is a function of 
time only. A simple form is 


1 -= 


where 


is a constant vector and G is a scaler function of time. 
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5.3.5 Calculation -procedure.- The solution for the rotorcraft aero- 
dynamics and dynamics proceeds as follows. The job begins with data input and 
initial calculations. Next the trim solution is obtained. Then the aero- 
elastic stability analysis, flight dynamics analysis, or transient analysis is 
performed as required. An old job can be restarted in any of these four tasks 
(trim, flutter, flight dynamics, or transient). 

In the trim analysis the controls are Iterated until the required operat- 
ing state is achieved. Since the nonuniform Inflow Influence coefficients 
depend on the rotor thrust (through the wake geometry) it is necessary to 
Iterate between the influence coefficient calculation and the trimmed motion 
and forces calculation (unless the rotor thrust is specified as part of the 
definition of the required operating state). The trim analysis is performed 
first for uniform inflow, then for nonuniform inflow with a prescribed wake, 
and finally for nonuniform Inflow with a free wake gecnnetry. After obtaining 
the trim solution, the aircraft performance and loads can be calculated (this 
is the trim restart entry point for an old job)-. 

In the transient analysis, the rigid body equations of motion are numeri- 
cally integrated. At each time step there is an Iteration between the non- 
uniform inflow influence coefficient calculation, and the calculation of the 
rotor and airframe motion and forces. 

In the flight dynamics analysis, the stability derivatives are calculated 
and the matrices are constructed that describe the linear differential equa- 
tions of motion. At each motion or control increment in the stability deriva- 
tive calculation there is an iteration between the influence coefficient 
calculation and the calculation of the motion and forces. Finally the system 
of linear differential equations is analyzed (optionally including a numerical 
integration as for the transient analysis) . 

In the flutter analysis the matrices are constructed that describe the 
linear differential equations of motion, and the constant coefficient or 
periodic coefficient equations are analyzed. Optionally the equations are 
reduced to just the aircraft rigid body degrees of freedom by assuming quasi- 
static response of the other degrees of freedcmi, and the equations are 
analyzed as for the flight dynamics task. 
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6. AEROELASTIC STABILITY 


The objective of the aeroelastlc analysis Is to derive a set of linear 
differential equations describing the perturbed motion of the helicopter from 
the trim flight condition. The stability of the system is defined by the 
eigenvalues of these equations. 


6.1 Rotor Model 

The differential equations of motion for the rotor blade have been 
deirved in section 2.2.18. Here it is necessary to linearize the inertial 
and aerodynamic forces in these equations. 

6.1.1 Rotor degrees of freedom.- The rotor blade motion is described by 
coupled flap/lag bending, rigid pitch and elastic torsion, and optionally the 
glnibal pitch and roll motion (or teeter motion for the two-bladed rotor case) . 
The blade degrees of freedom are written as the sum of trim terms and pertur- 
bation terms. The trim solution is described in section 5.1; the perturbation 
motions are the degrees of freedom for the aeroelastlc analysis. In particu- 
lar, the generalized coordinate of the 1-th blade bending mode is written 

v -4- 


or for the bending deflection 

After substituting for q^, the delta notation indicating the perturbed 
motion can be omitted. 

The rotor equations of motion have been obtained in the rotating frame, 
with degrees of freedom describing the motion of each blade separately. In 
fact, however, the rotor responds as a whole to excitation from the nonrotat- 
ing frame — shaft motion, aerodynamic gusts, or control Inputs. It is 
desirable to work with degrees of freedom that reflect this behavior. Such a 


C V 


.Po I- V 
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representation of the rotor motion simplifies both the analysis and the under- 
standing of the behavior. 

The appropriate transformation to obtain the degrees of freedom and 
equations of motion in the nonrotating frame is of the Fourier type. There 
are many similarities between this coordinate change and Fourier series, 
discrete Fourier transforms, and Fourier interpolation; the ccmnnon factor is 
the periodic nature of the system. A Fourier series representation of the 
blade motion is appropriate for dealing with the steady-state solution 
(section 5.1.1). Here we are considering the general dynamic behavior, 
including transient motions; hence the Fourier coordinate transformation is 
required. This coordinate transformation has been widely used in the clas- 
sical literature, although often with only a heuristic basis. For example, 
it has been used in ground resonance analyses to represent the rotor lag 
motion (ref. 30) and in helicopter stability and control analyses for the 
rotor flap motion (ref. 31). More recently, there have been applications of 
the Fourier coordinate transformation with a sounder mathematical basis 
(e.g., ref, 32). 

Consider a rotor with N blades equally spaced around the azimuth, at 
“ 'I* + (where = 2Tr/N and the blade index m ranges from 1 to N) . 
Hence tfi ■ At is the dimensionless time variable. Let q^™^ be the degree 
of freedom in the rotating frame for the m-th blade, m « 1 to N. The Fourier 
coordinate transformation is a linear transform of the degrees of freedom 
from the rotating to the nonrotating frame. Thus the following new degrees 
of freedom are introduced: 




- ^ 

v*. ~ 1 





v-^«. I i 
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Is 


Here is a collective mode, and Bj^g are cyclic modes, and ^N/2 

the reactionless mode. For example, for the rotor flap motion, B^ is the 
coning degree of freedom, while B^^^, and Bj^g are the tip-path plane tilt 
degrees of freedom. The inverse transformation is 

* ' Vl ' 


which gives the motion Of the individual blades again. The summation over 
n goes from 1 to (N-l)/2 for N odd and from 1 to (N-2)/2 for N even. 

The Bjij /2 of freedom appears in the transformation only if N is 

even. The corresponding transformation for the velocity and acceleration are 





-4- C ■2-V»^V\4.— 


Note that transformation to the nonrotating frame introduces Coriolis and 
centrifugal terms. 


The variables 


«o- ®nc’ ®ns’ ^ ®N/2 


are degrees of freedom^ that is, 
(m) 


functions of time, just as the variables q'“"' are. These degrees of free- 
dom describe the rotor motion as a whole, in the nonrotating frame, while 

describes the motion of an individual blade in the rotating frame. Thus 
we have a linear, reversible transformation between the N degrees of free- 
dom q^™^ in the rotating frame (m = 1, ...» N) and the N degrees of 


freedom (& , B , B„_, 
o nc ns 


^N/2^ in the nonrotating frame. Compare this 


coordinate transformation with a Fourier series representation of the steady- 

state solution. In that case, Is a periodic function of . so the 

m 

motions of all the blades are Identlclal. It follows that the motion in the 
rotating frame may be represented by a Fourier series, the coefficients of 
which are steady in time but infinite In number. Thus there are similarities 
between the Fourier coordinate transformation and the Fourier series, but 
they are by no means Identical. 

This coordinate transform must be accompanied by a conversion of the 
equations of motion for q^ ^ from the rotating to the nonrotating frame. 
This conversion is accomplished by operating on the equations of motion with 
the following summation operations: 

The result is equations for the , P , and 3„,_ degrees of freedom, 

respectively. Note that these are the same operations as are involved In 
transforming the degrees of freedom from the rotating to the nonrotating 
frame. Since the operators are linear, constants may be factored out. Thus 
with constant coefficients in the equations of motion, the operators act only 
on the degrees of freedom. By making use of the definition of the degrees 
of freedom In the nonrotating frame, and the corresponding results for the 
time derivatives, the conversion of the equations of motion Is then straight- 
forward, Complexities arise when it is necessary to consider periodic coef- 
ficients, such as due to the aerodynamics of the rotor in nonaxlal flow. 

The total force and moment on the hub have been obtained by summing the 

contributions from the individual blades. The result Is operators exactly of 

the form above, for obtaining the total hub reaction in the nonrotating frame 

from the root reaction of the Individual blades in the rotating frame. The 

origin of the summation operation is clear, and the sini|j or costi; factors 

m m 

arise when the rotating forces are resolved Into the nonrotating frame. In 
fact, the equation conversion operators in general may be viewed as simply 
resolving the moments on the individual blades into the nonrotating frame. 
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The Fourier coordinate transformation is often associated in rotor 
dynamics vith the generalized Floquet analysis. The latter is a stability 
analysis for linear differential equations with periodic coefficients. 

Indeed, there is a fundamental link between these topics because both are 
associated with the rotation of the system. They are, however, truly separate 
subjects — either can be required in the rotor analysis without the other. 

For example, a rotor in axial flow on a flexible support (or with some other 
relation to the nonrotating frame) requires the Fourier coordinate transforma- 
tion to represent the blade motion, but is then a constant coefficient 
system. Alternatively, for the shaft-fixed dynamics of a rotor in forward 
flight, a single-blade representation In the rotating frame Is appropriate, 
but there are periodic coefficients due to the forward flight aerodynamics 
which require the Floquet analysis to determine the system stability. 

For the present Investigation, the degrees of freedom to be transformed 
to the nonrotating frame are blade bending, blade pitch, and glnfbal motion. 

The nomenclature for the corresponding degrees of freedom In the rotating 
and nonrotating frames are as follows: 


bendng 

rotating 

nonrotating 
jO CS' C:> 0> 

torsion 

Pi”’ 

lO CO CO c.:^ 

©w/a 

gimbal 




The notation Is used for the 1-th bending mode In the nonrotating 

frame. With the modes ordered according to frequency, is thus usually 

( 2 ) 

the fundamental lag mode, and 3 the fundamental flap mode. Similarly, 
9^^^ is the i-th torsion mode, with 0^^^ rigid pitch and the remaining 
modes elastic torsion. The collective and cyclic modes (0, 1C, IS) are 
particularly Important because of their fundamental role in the coupled 
motion of the rotor and the nonrotating system. When the transformation of 
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the equations and degrees of freedom Is accomplished, for axial flow there 
is a complete decoupling of the variables into the following sets: 

(a) the collective and cyclic (0, 1C, IS) rotor degrees of freedom 
together with the gimbal tilt and rotor speed degrees of freedom 
and the rotor shaft motion 

(b) the 2C, 2S, nc, ns, and N/2 rotor degrees of freedom (as 

present) 

Thus the rotor motion in the first set is coupled with the fixed system, 
while the second set consists of purely internal rotor motion. Nonaxial 
flow couples, to some extent, all the rotor degrees of freedom and the fixed 
system variables, primarily due to the aerodynamic terms; still the above 
separation of the degrees of freedom remains a dominant feature of the rotor 
dynamic behavior. 

For a two-bladed rotor, the blade bending degrees of freedom are coning 
and teetering type modes: 


II 

JL 

10 

vir=i \ = 

t( 



= 

J- 

lO 

i (- = 

^ « 1 

i ( 




The pitch/torsion degrees of freedom 6^ and are similarly defined. 

The teetering degree of freedom is also Included for the two-bladed 

rotor (in place of the gimbal degrees of freedom) . The teetering motion is 
defined in the rotating frame, hence 

Be, = O 

The special characteristics of the two-bladed rotor dynamics are reflected 
in the appearance of the teetering-type degrees of freedom 6^^, and B,^) » 

rather than the cyclic motions (1C and IS) as for N > 3. The coning and 
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teetering equations of motion are obtained by applying to the rotating frame 
equations the following operators: 








6.1.2 Rotor equations and hub reactions,- The equations of motion for 
the coupled flap/lag bending and for elastic torslon/rlgld pitch motion of 
the blade in the rotating frame (section 2.2.18) are linearized for the 
aeroelastlc analysis. The result is: 
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M>«. 


'*s 








The Inertia constants are defined In section 6.1.3. Only linear lag damping 
has been considered here, 2 ind for convenience the lag damper term Is Included 

if 

In the coefficient Ig^q^. Introducing the Fourier coordinate transformation 
for the blade degrees of freedom, the rotor hub forces and moments derived In 
section 2.2.18 become: 
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The equations of motion for glmbal tilt and roll, or for the teeter motion 
of a two-bladed rotor, are obtained from 








- 


or 


— Tt 




Css Ocs -o 


<zT i 




vhere 




Z C-.')" V 5=’ 

2. ^ A I ^Ck 


«> 


— 31^*^ ( f ^ -1- 



.* 

_ 

— ^ 

Cp. 



-- 

( 




^ f p 


r.) 


(section 2.2.18). Also, for a two-bladed rotor, the hub reactions take a 
somewhat different form: 
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The aerodynamic forces required are 
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6.1.3 Inertial aonstants.- The Inertial constants for these linearized 
equations of motion are obtained from the constants defined In section 
2.2.19 as follows: 
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6.1.4 

Aerodynamic forces.-- The blade section forces and pitch moment as 

derived in 

section 

2.4.1 are: 



OkO 
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U (•'f -t- Ut g ^ 
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Each component of the velocity seen by the blade has a trim term and a small 
perturbation term, so we write 



© H— 

IAt- ^ t 

Up -4- 


It follows that the perturbation of angle of attack, resultant velocity, and 
Mach number are: 

=- /CA 

StA •=• ^‘^-r Sut -f- 
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and the perturbations of the aerodynamic coefficients are 






i-h 






with similar results for c and c,* The perturbations of the section aero- 

m d 

dynamic forces may then be obtained by carrying out the differential operation 

on the epxressions for F , F , F , and M , using the above results to express 

z X r a 

the perturbations in terms of 66, 6u^, 
the perturbation quantities are evaluated at the trim state. The results for 
the perturbation forces are: 
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The blade trim velocity components are defined in section 2.4.2. The 
perturbation velocity components are due to the blade degrees of freedom, the 
shaft motion, and the aerodynamic gust velocity; 
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The gust velocity components are here assumed to be uniform throughout space. 
Perturbation Inflow components X^, X^, and X^ have been Included in 6up 
(see section 6.1.5). Recall that the body Euler angle contributions are not 
Included in the evaluation of a^, a^, and here. The perturbation quanti- 
ties required for the unsteady pitch moment are 
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(see section 2.4.8). The derivatives of the blade section aerodynamic coeffi- 
cients with respect to angle of attack and Mach number are obtained from 
steady, two-dimensional airfoil characteristics with corrections for tip flow, 
yawed flow, and dynamic stall effects, as described in sections 2.4.4 and 

2.4.7. 
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Combining the expansion for the section forces and moment in terms of 
the velocity perturbations, and the velocity in terms of the motion of the 
rotor and shaft, the perturbation aerodynamic blade forces expanded linearly 
in the degrees of freedom are obtained. Giving names to the aerod3mamic 
coefficients at this point in the analysis, the results for the required 
aerodynamic forces on the rotating blade are as follows. The aerodynamic 
force for flap/lag bending is: 

H- UoT) ^ 
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The radial force is 
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The aerodynamic force for blade torsion and pitch Is: 
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Finally, the aerod3mamlc hub forces and moments are similar to the result 
for the blade bending, but with the following changes In the Integrands and 
notation: 


Flap moment 
Torque 

Blade drag force 
Thrust 


Integrand 


rF 


rF 

X 

F 

X 

F 

z 


Coefficient Notation 
M 

Q 

H 

T 


- 330 - 



Combining the results for the expansion of the aerodynamic forces and the 

expensions of the velocities, the aerodynamic coefficients can be evaluated. 

The aerodynamic coefficients are constant for axial flow, but for nonaxial 

flow they are periodic functions of ^ . The coefficients for flap/lag bend- 

m 


ing are: 
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The aerodynamic coefficients for the flap moment are 
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The aerodynamic coefficients for 
pattern of the flap moment, with 
integrands : 

the other hub forces and moments follow the 
the following changes in the notation and 


Integrand 

Coefficient 

Flap moment 

rF 

z 

M 

Torque 

rF 

X 

Q 

Blade drag force 

F 

X 

H 

Thrust 

F 

T 


z 
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The radial force coefficients are 

^ = So ^ 

. I ^ 

= S. 

K = ^ 

^ = si Pc, 

^ — <*e ~ 

■= *i>o (Pct''^-'^; -t- 4^ 

"* S, (Ppr*^ ^ ' -t- Pr, "4 •'!''> 4 j- 

+ S'„Pc, ■] ^ 
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where 

— ^xt t ^^3 "^' ”+ -^“^y 3 

and Fj.p and Fj.g are similarly defined. Finally, the aerodynamic coeffi- 
cients for the blade pitch and torsion are 

Vp^txTz ^ 4.T- 

-I ^a. 

3i 0-+-4 ^!^^ua 

it ii ^ )2.Ufc c-,a »igf> V ^ ^ 

~ ^vh Wl (l-t-'+^Y -i^ 

■'■ ^ (^l + 4^'>Ufc<a>sO &y.V S-ir- 





C,* U» - f ■*- V^'> 1 

-- i-ft C UK, - rr,^-4 + f„t»'>- 

***• ^‘eik t U*^» — "ii 

(§t • '/"K'V" 

Sr^ U "Iz. W I C * ■*■ '^ "^ • f“/ ■£. 

“•■ SrB^ U fz.("\ i*# '*-'Y‘*3*’J’'''^ — 

= Sr„,U,K, - CK,T*-+-P-*e‘^'>-%^l^i.lr 

^< 5 * U TC ‘^A \ (^l+ 4 ^')(' + 2 .^') ^ If— 

/^f,K « __ S^Ull -ft W1 C'+^%>6+ 2.%^ I;;. 4^ 


where 
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6.1.5 Inf loti} dynamioB.- The aerodynamic forces on the rotor result in 
wake-induced Inflow velocities at the disk, for both the trim and transient 
loadings. The wake-induced velocity perturbations can be a significant 
factor in the rotor aeroelastic behavior; an extreme case is the influence 
of the shed wake on rotor blade flutter. The rotor inflow dynamics should 
therefore be included in the aeroelastic analysis. However, the relationship 
between the Inflow perturbations and the transient loading is likely more 
complex even than for the steady problem (nonuniform wake- induced inflow 
calculation) , and models for the perturbation Inflow dynamics are still under 
development. In the present analysis, an elementary representation of the 
inflow dynamics is used. The basic assumption is that the rotor total 
forces vary slowly enough (compared to the wake response) that the classical 
actuator disk results are applicable to the perturbation as well as the trim 
velocities . 

A contribution to the velocity normal to the rotor disk of the following 
form has been Included in ^u^: 

=. -4- 


where X is the inflow perturbation component iinlform over the disk, while 
the and X components vary linearly over the disk. The inflow dynamics 
model must relate these inflow components to the transient aerodynamic 
forces on the rotor, specifically to the thrust, pitch moment, and roll 
moment; and to the transient rotor velocity perturbations 6p , 6y and 6y . 
Following reference 33 we use: 
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where 


<^-r ~ Cr-r- 






and 

<^\ 

<^“v 


dM 



Here we have included linear velocity perturbations due to the thrust, con- 
sistent with the trim Inflow model of section 2.4.3, lAlch gives expressions 

for the constants k and k . These relations for the Inflow perturbations 

X y 

imply the following lift deficiency functions: 



for moments 


for thrust 



forward flight 


moments In hover 


thrust in hover 
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(see ref. 33). A time lag In the inflow response to loading changes will 
also be Included: 






f&r 1 

^ 1 


V. 





i j 





1 ‘^1 


i / 


using 6X/6T and 6X/6M, with the constants ■ .85 and 

K = .11 (refs. 34 to 36). These relations give a dimensional time lag of 
M 

T S! .22/Xn in hover, and in forward flight x^ ■ .42/y{2 and x^ ■ ,22/yfi. 

The effect of the ground on the inflow dynamics is to add a perturbation 
due to changes in the rotor height above the ground: 

SX - a s. 

As for the trim inflow analysis, the result of reference 4 for the ratio of 
the Induced velocity in and out of ground effect is used: 

^ ZL_ ^ I 

^ I / 

which gives , 

di. y e* 

Expressions for cose and z are derived in section 2.4.3. Since 6X/6z > 0, 
ground effect introduces a positive spring to the rotorcraft flight dynamics. 

A decrease in the rotor height above the ground produces a decrease in the 
induced velocity, hence a rotor thrust increase that acts as a spring against 
the vertical height change. 

For the side-by-side helicopter configuration, the antl 83 nnmetric dynamics 
exhibit an unstable roll oscillation due to interaction of the rotor wake 
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and ground. Such behavior can be Included in the ground effect model derived 
here by using a negative value for 6 X/ 6 z (a negative roll spring), which must 
be obtained from experimental data. In this case the Inflow perturbations 
of the two rotors are related to the symmetric and antisymmetric height 
perturbations : 

+ tell, 

5.\ » ^ 

* » bii ^ 5 55 * 

X / \ _4_ -L f ^ ~ Ss, 


where and 6 Z 2 are the height perturbations at the two rotor hubs. A 

form applicable In general Is 






^■2 i-aL 



S »\2 







including the factor Rj^/R 2 si^ce the hub motion is normalized using Rj^. 

Finally, the rotor /rotor interference is Included in the inflow dynamics 
model, using the same interference factors Kj ^2 *^21 ^ trim 

induced velocity model (see section 2.4.3). 

In summary, the differential equations for the inflow perturbations 
X and X are: 

X V 
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for rotor #1 and for rotor #2. The coupled equations for the uniform inflow 
perturbations of the two rotors are: 



The velocity perturbations 



and 5y 

z 


are required for 





- V 2Ct 
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In shaft axes, the shaft motion and gusts give 

We, 

The rotor height perturbation 6 is obtained from the vertical component 

z 

of the displacement at the rotor hub: 

S-t = — Vt • ( -4- 

= -"^6 - ( f 

expressing the hub motion in terms of airframe .degrees of freedom. For the 
rigid body degrees of freedom the mode shapes are: 

it - -Ik 1 = [ 3 

(section 4,2.1). Hence 
“ 

S^&r-r 

— ^ op 

p ('cos©^^ 

-h ([c-cxad^pX" 
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Note that a spring is Introduced Into the z_ equation (if cosB coscfi ^ 0) , 

F FT FT 

and also possibly into the y^, and Xp equations. 


The time lag is often not an important factor, so a quasistatic model 
for the Inflow dynamics is generally sufficient. Dropping the time lag terms, 
the equations for X^, and X^ reduce to linear algebraic equations. 

Hence in the quaslstatlc case the Inflow perturbations do not Increase the 
order of the system. The wake Influence reduces to an algebraic substitution 
relation, which if Incorporated analytically would lead to lift deficiency 
functions; with large-order systems. It is more practical to accomplish the 
substitution numerically . 


6.1.6 Rotor equations of motion.- The linear differential equations of 
motion for the rotor model can be constructed now. The equations of motion are 
In the nonrotating frame, that is the Fourier coordinate transformation has 
been applied to the bending and torsion degrees of freedom of the blade. For 
now only a three-bladed rotor Is considered; the equations are extended to an 
arbitrary number of blades below. The equations of motion for the rotor, and 
the hub reactions, take the following form: 


+ AiiXr + AoXit -i-ktP 

CT c — ^ •• ^ 

= '-'Z Kr -f- -I- Vvuu» 

The coefficient matrices are constructed from the results of section 6.1.2. 

Here the matrices only Include the structural and inertial terms; M and 

aero 

F are the aerodynamic forces. The vectors of the rotor degrees of 

freedom (x^) , shaft motion (a) , rotor blade pitch input (v^) , aerodynamic 
gust (g » in shaft euces) , and the hub forces and moments (F) are defined as: 


[ 

■ ^ 

r>i ©o ©,4 ^sc 

An 


r 



Vj = 

L©. 



= 





- 3 ^ 2 - 





C. '^JU Ofyi 







Note that in the rotor degrees of freedom the notation and 

is Intended to cover as many bending and torsion modes as the analysis 
requires. Also^ the degrees of freedom used for the inflow dynamics model 
^y by A *= X (in order that the highest order 

derivatives wil be A, in the acceleration matrix) . The inertial matrices 
are defined in section 6.4.1. 


The aerodynamic terms M and F are required to complete the 

aero aero ^ 

differential equations of the rotor model. They are obtained by summing 
over all N blades the aerodynamic forces in the rotating frame (section 6.1.4) 
and introducing the Fourier coordinate transformation for the blade bending 
and torsion degrees of freedom as required. The result for the required 
aerodynamic forces is 

H- 4-Aoot — 



For the case of a rotor operating in axial flow (y * 0) the aerodynamic 

coefficients for the blade forces in the rotating frame are constants, 

independent of the blade azimuth angle • The coefficients are also then 

m 

entirely Independent of the blade Index (m) ; hence the summation over the 
N blades operates only on the system degrees of freedom, not on the aero- 
dynamic coefficients themselves .(which factor out of the summation). The 
resulting coefficient matrices, which are constant for axial flow, are 
defined in section 6.4.2. 

For the case of a rotor operating in nonaxlal flow (y > 0) the aerodynamic 
coefficients of the rotating blade forces are periodic functions of 

m 

because of the periodically varying aerodynamics of the edgewise moving rotor. 
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It follows that the rotor In nonaxlal flight Is described by a system of 
differential equations with periodic coefficients. It Is possible to express 
the aerodynamic coefficient of the rotating blade forces as Fourier series, 
and then to obtain the coefficients of the nonrotating equations In terms 
of these harmonics. However, the simplest approach for numerical work with 
large-order systems Is to leave the coefficients of the nonrotating equations 
In terms of the summation over the N blades of the rotor. The summation 
Is easily performed numerically, and It Is found that this form Is also 
appropriate for a constant coefficient approximation to the system. For 
nonaxlal flow, the coefficient matrices are periodic functions of the blade 
azimuth angle *" 'I' + = 2 tt/N. The period is ■=> (2/3)u 120® 

for the N ■= 3 case considered here. The coefficient matrices for nonaxlal 
flow are defined In section 6.4.3. 

The rotor equations as constructed here are not entirely complete. 

First, the rotor aerodynamic thrust and hub moments: 

where . i >. 

^ (iSZCrl^\^ 

have been put In place for the X , X , and X equations. Because of the 

u X y 

rotor /rotor interference and ground effect. It Is appropriate to finish the 
construction of these equations at a later stage (section 6.3.1). 
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Secondly, we have substituted in the equation of motion for rigid pitch: 




However, the rotor pitch control here still Includes the governor and 

mast bending terms, as well as external control Inputs. 

Thirdly, the rotor torque yC^/aa has been put in place as the equa- 
tions of motion for the rotational speed perturbation The drive train 

couples the rotational speed perturbations of the two rotors, so it is 
necessary to construct these equations at a later stage (section 6.3.1). 

Consider now the case of a rotor with four or more blades. Each rotat- 
ing degree of freedom of the blade (bending or torsion motion) must result in 
N degrees of freedom for the rotor as a whole. Thus increasing the nunter 
of blades adds degrees of freedom and equations of motion to the rotor de- 
scription. In axial flow these additional degrees of freedom do not couple 
with the collective and cyclic degrees of freedom of the rotor. Hence the 
equations given above remain valid for rotors with N > 3 also, and we need 
be concerned here only with the equations of motion for the additional degrees 
of freedom. These additional degrees of freedom are not coupled inertially 
with the shaft or glmbal motion. The additional equations of motion for 
bending and torsion of a rotor blade with four or more blades are then: 

A2 -4- ”4- ~ Kcuuxi 

with the vectors of the degrees of freedom and blade pitch control here 
defined as follows: 

fW> W W (tv _-r 

— L ©M, ©W/2 J 

The inertial coefficient matrices are given in section 6.1.4. 
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The aerodynamic forces are required to complete the equations of 
motion. In axial flow the aerodynamic forces still do not couple the addi- 
tional degrees of freedom for N > 4 with the shaft or gimbal motion. Hence 
the aerodynamic forces for axial flow take the form 


— R 


our^ 


— - 4 - Kgi 


with the coefficient matrices defined in section 6.1.5. 

The aerodynamic forces in nonaxlal flow (y > 0) couple all degrees of 
freedom of the rotor with each other and with the shaft and gimbal motion. 
Then not only are additional degrees of freedom and equations of motion 
Involved if N > 3, but the number of blades also Influences the equations 
and the hub reactions given above. Rather than directly presenting the aero- 
dynamic matrices for the general case of three or more blades In nonaxlal 
flow, the analysis is extended by means of an observed pattern In the coeffi- 
cients. In the nonaxlal flow equations (section 6.4.3), note the repeated 


occurrence of the 

following submatrices: 









\ 


s, 




2C, 

2Lcf 

2lc,S, 















1 







o 

-s, 

<=1 


bb= 


o 

-ZC.S, 





D 


2C.S, 




- 




(using 

the notation S * 
n 

sin n4^ and C “ 
m n 



s.') 



-S. 


cos ni|< ) . 
m 


These matrices are a 


direct result of the introduction of the Fourier coordinate transformation 
(columns) and the application of the summation operators to obtain the non- 
rotating equations (rows) . The matrix DP arises from application of the 
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Fourier transformation to the time derivatives (q. or p ). In the B and X 

lx G 

matrices, only some columns of P and DP appear, while in the C matrices 
only some rows appear. The extension to an arbitrary nuinber of blades 
(N > 3) is then simply 



(-\T I 



f ' I 

f aQ 

I 

I 





I 




V^Cvi 



Rotors with three or more blades may be analyzed within the same 
general framework, but the two-bladed rotor is a special case. The rotor 
with N > 3 has axi-symmetric inertial and structural properties and hence 
the nonrotating frame equations have constant coefficients in axial flow. 
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In contrast, the lack of axi-symmetry with two blades leads to periodic 
coefficient differential equations, even in the Inertial terms and in axial 
flow* Only in special cases are the dynamics of a two**bladed rotor described 
by constant coefficient equations. The equations of motion again take the 
form 


with now for N => 2 the rotor degrees of freedom and pitch input defined as 
follows : 

= t P- h d, ©. f , Hi jU JU JL., ] 


The Inertial coefficient matrices are defined in section 6.4.6. Note that 
there are periodic coefficients in the matrices coupling the rotor and shaft 
motion (A, C , C) . 


The required aerodyanmic forces for the two-bladed rotor case again take 
the form 







C«-» ^ 


The aerodynamic coefficient matrices are defined in section 6.4.7. 

An Independent blade analysis is useful for problems not involving the 
shaft motion or other excitation from the nonrotating frame. The only rotor 
blade degrees of freedom involved are the bending and torsion nratlon. The 
shaft motion, gimbal motion, and the rotor speed perturbation are dropped 
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from the system • Only a single blade need be analyzed, in the rotating frame. 
The equations of motion for the bending and torsion modes are then: 


- f <pv ^ r- 

~ f V - ~ ^ 

( h + M f h -<- ?• 

- F <T T' - F - ><r : '=i'. 

- F - f ^ 

or 

it 

T 

with * (q, p- ) and - 6 . These equations can also be obtained by 

R k k R con 

dropping all degrees of freedom except the collective modes from the analysis 
above; a separate construction for the independent blade case is more effi- 
cient however. 

6.1.7 Constant coefficient approonmation*- The rotor dynamics in non- 
axial flow are described by a set of linear differential equations with 
periodic coefficients. A constant coefficient approximation for nonaxial flow 
is desirable (if it is demonstrated to be accurate enough) because the calcu- 
lation required to analyze the dynamic behavior is reduced considerably 
compared to that for the periodic coefficient equations, and because the 
powerful techniques for analyzing time-invariant linear differential equations 
are then applicable. However, such a model is only an approximation to the 





correct aeroelastlc behavior. The accuracy of the approximation must be 
determined by comparison with the correct periodic coefficient solutions. 
The constant coefficient approximation derived here uses the mean values of 
the periodic coefficients of the differential equations in the nonrotating 
frame. 

To find the mean value of the coefficients, the operator 

^ C 


is applied to the periodic aerodynamic coefficients (given in section 6.4.3), 
resulting in terms of the following fonn for the N * 3 case: 


c- 


Zrr 





__ 1 , £ 

— to 4 






\ 


tc 



^ < 1 


1 / 


the rotating blade aerodynamic coefficient M: 

.= -4- M 

VI d 


In the present case, these harmonics must be evaluated numerically. The 
aerodynamic coefficient M is calculated at J points, equally spaced around 
the azimuth. Then the harmonics are calculated using the Fourier interpola- 
tion formulas : 
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J- -5- 

— ^ -f- 

Mll'P 


II 

Mh 

M. 




= jA'l' = 

j(2-ir/j) (j = 1, .. 

. , J) . The number of harmonics required 


is n = N-1 for N odd and n = N-2 for N even (N is the number of blades) . 
Good accuracy from the Fourier interpolation requires at least that J = 6n. 
Using these Fourier interpolation expressions, the required harmonics are 





2^ 


I 

/ 



oosH>j 

2 


It follows then that the constant coefficient approximation is obtained 
from the periodic coefficient expressions by the simple transformation; 


KJ 






-L 






The summation over N blades (m « 1, N; = 2*rT/N) for a periodic 

coefficient is replaced by a summation over the rotor azimuth (j®l> ...jJ; 

= 2tt/J) for the constant coefficient approximation* This is quite con- 
venient since the same procedure may be used to evaluate the coefficients for 
the two cases » with simply a change in the azimuth increment. The periodic 
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coefficients must be evaluated throughout the period of = 0 to 2 it/N; the 
constant coefficient approximation (mean values only) is evaluated only once. 

With the substitution (1/N) I (1/J) T, the results given in section 
6.4.3 for the periodic coefficient matrices are directly applicable to the 
constant coefficient approximation as well. 


For the case of a rotor with four or more blades, the constant coeffi 
dent approximation involves the transformation of higher order harmonics: 
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So the periodic coefficient results are still applicable to the constant 
coefficient approximation if the summation over the N blades is replaced 
by a summation around the rotor azimuth. 

This transformation is also applicable to the case of a two-bladed 
rotor, but the constant coefficient approximation is not as useful or as 
accurate for N = 2 as for N > 3. With three or more blades, the source 
of the periodic coefficients is nonaxial flow, hence the periodicity is order 
\i or smaller. At low advance ratio then, the constant coefficient approxi- 
mation may be expected to be a good representation of the correct dynamics. 
The two-bladed rotor has also periodic coefficients due to the Inherent lack 
of axi-symmetry of the rotor. This periodicity is large even for axial flow, 
and neglecting it in the constant coefficient approximation may be a poor 
representation of the dynamics. In particular, it is not possible to use the 
constant coefficient approximation as formulated here for the flight dynamics 
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analysis of a two-bladed rotor helicopter, since this averaging eliminates 
the coupling between the rotor and the shaft motion, 

6.2 Aircraft Model 

6.2.1 Aircraft degrees of freedom.- The aircraft motion is described 
by the rigid body and elastic airframe degrees of freedom: 

~ ~ ^ V 'y ^ 

as defined in section 4.2.1. The aircraft controls consist of flaperon, 
elevator, aileron, and rudder deflections: 

V* = L ^ 5. ^ 

The rotor hub motion is obtained from 

«=r C 

where c is defined in section 4.2.2, including the sign changes for a 
clockwise rotating rotor and scaling for rotor #2. Recall that the Euler 
angles do not contribute to a , a , and a however. In addition there is a 
linear acceleration due to the rotation of the velocity vector in body axes 
by the Euler angular velocities, written 

(see section 4.2.2). 

The feedback of the airframe elastic motion to the rotor cyclic pitch is 

~ ~ Si 3 

as defined in section 4.2.3. 
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6.2.2 Aircraft equations of motion.- The linearized equations of motion 
for the aircraft rigid body and elastic motion are: 



(see section 4.2.4). The generalized forces due to the rotor hub reactions 
are 



where c^ is defined in section 4.2.5, including the sign changes for a 
clockwise rotating rotor and the scaling for rotor #2. 


The generalized forces due to the aircraft aerodynamics can be linearized 
by succesively perturbing the inputs to the analysis described in section 
4.2.6. Hence for the rigid body degrees of freedom we obtain: 
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(The coefficients of the matrices are the aircraft stability derivatives, due 
to the wing/body, horizontal tail, and vertical tail.) A gust velocity 
uniform throughout space is considered; hence the gust velocity components 
are the same at the wing/body and tail. The mean Inflow perturbations 
influence the airframe through the rotor-induced aerodynamic interference, 
which is modelled as described in section 4.2.6. So the interference 
velocities at the wing and tail are obtained from 








From section 4.2.7, the generalized forces for the airframe elastic degrees 
of freedom (k > 7) are: 



M j.r_ 

<r« A L 


V 



Hence the aeroelastic motion of the helicopter airframe is described by a set 
of linear, constant coefficient differential equations of the form 


•• 

'll 


Xc -+■ X, 


Wvj -4- 

“I- eT R -4- 
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where x Is the vector of aircraft degrees of freedom, v is the vector of 
s s 

aircraft control variables , g is tha vector of gust velocity components 
(in wind axes), and and F 2 are the hub reactions of the two rotors. 


6.2.3 Drive train equations of motion.- A model for the transmission 

and engine dynamics was derived in section 4.3.2. The degrees of freedom 

involved are the rotational speed perturbations of the two rotors (ij and ii ) 

si s2 

and the engine speed perturbation (degree of freedom i(> , defined relative to 

e 

the rotation of rotor #1) . With the coupling of the rotors by the drive 
train, it is more appropriate to use the degrees of freedom defined by 




where ratio of the trim rotational speeds of the two 

rotors. So ip is the rotational speed perturbation of the rotors and drive 
train as a whole, while represents differential rotation of the two 

rotors. The degrees of freedom and i|/^ therefore involve elastic deflec- 
tion of the drive train. The engine model introduces the throttle control 
variable 6^. The linearized equations of motion for the drive train are 
then : 



/^Cq\ 

V ' 2 . 
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rjda C*V'i + V*,') 

*+*«. ~ — 

The constants appearing in these equations are defined in section 4.3.2 for 
several drive train configurations. 

For the autorotation case the degree of freedom is dropped , and the 

Q., Q , and I„ engine terms in the i|^ and equations are omitted. For 
at t E Si 

the engine out case, with the engine and rotors still connected, the engine 
terms and Q are omitted. The case of constant rotor speed is treated 

db tZ 

by dropping the and ip degrees of freedom from the system. 

S i 

The rotor speed governor, consisting of i> and ^ feedback to the 

s s 

engine throttle and to the collective pitch of each rotor, is described by 
the following equations: 

" - . 

Tri + n:„ 

(see section 4.3.3). For the tilting proprotor configuration, the variable 
ip is replaced by the symmetric variable ~ 

6.3 Coupled Rotor and Aircraft 

6.3.1 Coupled equations of motion.- The equations of motion have been 
derived for the two rotors and the aircraft body. The rotor equations of 
motion take the form 

A 

4- *4- A|(W 4> Ao^ 
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^ -I- CoXjt 4- Cfc&t -t- ^5 

(see section 6.1.6). The rotor degrees of freedom vector consists of 

flap/lag bending, rigid pitch and elastic torsion, gimbal or teeter motion, 
rotational speed, and Inflow perturbations. The rotor control vector v 

R 

consists of the blade pitch control. The gust vector in shaft axes is 
related to the gust vector in velocity axes by 




(section 4.1.4). The hub motion is related to the aircraft degrees of 
freedom by 


C.X 


S -I- o 


(section 4.2.2; only the 3x3 submatrix in the upper left comer of c is 
nonzero) . For rotor #2 it is necessary to change the time scale to the 
rotational speed of rotor #1. 

a _i 


So the matrices A^, C^, and are multiplied by (fij^/fi 2 ) ; and the 

matrices A 2 , A 2 , C 2 , C 2 are multiplied by 


The aircraft equations of motion take the form 


-»-<XiXs =a Wvc + W 






<ii fi, F* 
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(see section 6.2.2). The vector x consists of the aircraft rigid body 

8 

and elastic airframe degrees of freedom. The vector consists of the 

aircraft controls. 

The equations for the rotor and aircraft can now be combined to con- 
struct the set of linear differential equations that describes the dynamics 
of the complete system. These equations take the following form: 




ky i 









The state vector x, control vectors v and Vp, and the gust vector g are 
defined as: 

^ ?€ic. 

Pw: ^ 

Yr 

V — [ ^ ^ 

=- ^ SkX "J 

^ =: [ U<w Vg Wg 


The vector of the degrees of freedom for the entire system (x) consists of the 

degrees of freedom of the two rotors and the aircraft. The rotational speed 

degrees of freedom of the two rotors are replaced by the coupled degrees of 

freedom and and the engine speed degree of freedom is added. The 

si ® 

governor dynamics introduce the degrees of freedom ^®govri* ^^govr2* 

The vector of control variables for the entire system (v) consists of the 

blade pitch of the two rotors, the aircraft controls, and the engine throttle. 

The vector of the pilot ^s controls (v ) is related to the individual control 

p 

inputs by the linear transformation 
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cjPe 


dsfincd in SBCtlon 4«X»6« ThS' S6irodynssilc gust vsctor (g) 
is in velocity axes. 

The coupled equations of motion are obtained by substituting the hub 
motion into the rotor equations and hub reactions, and then the hub reactions 
into the body equations of motion. The resulting coefficient matrices for the 
coupled system are: 





















w 
















Recall that the Euler angle contributions are not to be included in 

Oy, and a^. Hence in constructing it is necessary to skip the angular 

hub motion (a , a , a ) columns of A and C for the Euler angle 0 _, 

X y Z O O r r 

ipp) columns of c (and A^) when evaluating A^c and C^c. 

The rotor mass will be included in the helicopter gross weight. Hence in 
constructing A2 it is necessary to skip the linear hub motion (x^, y^, z^) 
columns of C2 for the rigid body 6 ^, 4 ;^, y^, z^,) columns of c 

(and A2) when evaluating C2C. Also the rotor mass is included in the 
generalized mass of the airframe free vibration modes; so the linear hub 
motion columns of C2 will be skipped for all the body degrees of freedom 
columns of c when evaluating C«c« This approach is also required for the 

T-w _ 

c C2C term in A^|^. Since the term involving c is the linear hub accelera- 
tion due to Euler angle velocities, which has already been included in a^ 
if the gross weight includes the rotor mass, it follows that the entire 
c C2C term is to be dropped. 

The construction of the coupled equations of motion is still incomplete. 
The matrices defined above basically account for the coupling of the rotors 
and aircraft through the rotor hubs. It remains to account for the coupling 
which occurs through other paths. 
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Frequently the rotor Is modelled as having a rigid control system. This 
option requires some restructuring of the equations of motion, since the 
rotor equations have been derived assuming that the blade rigid pitch 
degrees of freedom are present in the model and that the blade pitch control 
inputs enter through these degrees of freedom. In the limit of infinite 
control system stiffness, the solution of the rigid pitch equation of motion 
reduces to 



for N > 3; and 


(S:\- (t-L, - f k.; (Sa -<^( 1 ) 
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for N = 2. The blade rigid pitch motion In this limit consists of the 
commanded control Input, feedback of the bending and glmbal motion due to the 
kinematic coupling, and a pitch change due to the azimuth perturbation with 
a fixed swashplate. Substituting for p^, the pitch/bending, pitch/gimbal , 
and pitch/azimuth coupling requires operations on the colinons of A^, as 
defined by the above equations. Next the rotor terms in the control matrix 
B are reconstructed from the rigid pitch columns of since the blade 

pitch motion becomes a control variable rather than a degree of freedom. Then 
the equations of motion for the rigid pitch degrees of freedom are dropped 
from the system. 


So far the aerodynamic hub forces have been put in place for the X^, 

X , and X equations of the two rotors; a = cx has been substituted for 
X y s 

the hub motion; and the 'time scale for rotor #2 has been changed to 
Completion of the equations of motion for the Inflow dynamics requires the 
following steps . 


a. Multiply the thrust by (aa/2Y)3X/3T and the moments by (aa/2Y)3X/8M. 

b. Add K times the X equation to the X equation, and < times 

the X^ equation to the X^ equation. 

c. Construct the ground effect terms (body motion contributions to 
~(3X/3z)^z) in the X^^^ and X ^2 equations. 

d. Account for the coupling of the X^j^ and X ^2 equations due to rotor/ 
rotor interference. 


e. Construct the diagonal terms; 1 in A^^ and x “ 3X/3M or K 9X/3T 


in A 2 (times (Rj^/n 2 ) 


for rotor #2) . 


f. 


Construct the aerodynamic Interference terms 

jC. 




In the aircraft equations. 


So far the torque has been put In place for the rotational speed equations 

of the two rotors; a = cx has been substituted for the hub motion; and the 

s 
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time scale for rotor #2 has been changed to Completion of the equations 

of motion for the drive train dynamics requires the following steps. 

a. Transform the , and i|) - columns of the matrices according to 

si s2 

^S, - 

to convert from and degrees of freedom to and 

degrees of freedom. 

b. Combine C.^ and C„ as required for the ij) and equations. 

Ql Q2 s I 

c. Construct the equation of motion for and complete construction 

of the 4) and equation by adding inertial, damping, and spring 
8 1 

terms. Construct throttle control terms in ip and \p equations. 

s e 

d. Construct the governor equations, and the governor degree of freedom 
terms in the other equations of motion. 

It is still necessary to account for the mast bending and governor feed- 
back terms in the pitch control 6 . The pitch/mast bending coupling is 

con 

(for each rotor) 



for N > 3, and 

60 \ 

A©1 ==^ — ^ C W'H*) 

for N = 2. 

For a two-bladed rotor the pitch control degrees of freedom are 0^ and 
It is also useful to consider conventional cyclic control, which gives 

s=z ©le c-osy -I- Bk 
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The © and O columns of the control matrix can be constructed from the 
Ic Is 

column according to this equation* 

Next, the pilot's control matrix is constructed from the individual 
control matrix by the linear transformation Bp = 

Finally, the unused equations of motion and degrees of freedom are 
eliminated from the model by deleting the appropriate rows and columns from 
the coefficient matrices. Note that a number of the degrees of freedom are 
first order (no spring terms) : all the Inflow perturbation variables , the 

rigid body degrees of freedom ij'p* Xp, yp, and Zp (except possibly in ground 
effect) , and perhaps the rotational speed degree of freedom (in axial flow 
with no integral governor) . 

6.3.2 Quasistatia approximation,- It is frequently possible to reduce 
the order of the system of equations describing the rotorcraft dynamics by 
considering a quasistatic approximation for certain of the degrees of freedom. 
Assume that the equations of motion have been reordered so that the quasistatic 
variables (Xq) appear last in the state vector: 



The quasistatic approximation consists of neglecting the acceleration and 
velocity terms of the quasistatic variables. Thus the equations of motion 

t-V» ^ * 



The quasistatic variables now are not described by differential equations but 
rather by linear algebraic equations. The solution for x^ then is simply 

X. = (^: y L B' V - _ x-; i. -A*j X. 3 
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Substituting for Xq in the x^ equations of motion gives then the reduced- 
order equations for the quasistatic approximation; 

CA‘k-K:co''A» [X:-A':(Ary'A? ix. 

+ [a“-a‘: (ATT* a;’ = [B' - a'.*(AVT‘ B*3v 


In the present analysis, the quaslstatlc approximation can be applied to 
the inflow dynamics of either or both rotors, to the rigid pitch/elastic 
torsion degrees of freedom of either or both rotors, to all the degrees of 
freedom for either or both rotors, or even to all the degrees of freedom 
except the rigid-body motions of the aircraft. 

The quasistatic approximation retains the low-frequency dynamics of the 
eliminated degrees of freedom. Whether that is a satisfactory representation 
of the elastic behavior must be established by comparison with the results of 
the higher order model. 

The quasistatic approximation as Implemented here does not give the low 

frequency response when applied to a two-bladed rotor. The source of this 

(k) (k) 

difficulty is the fact that the teetering equations of motion ’ 

B ) are really still in the rotating frame, so the response of the teetering 
modes to low frequency inputs from the nonrotating frame is not at low 
frequency also, but rather at frequencies around 1/rev. 

6.3.3 SymmetriQ aircraft.- With lateral symmetry of both the aircraft 
and flight state, the symmetric and antisymmetric motions are completely 
decoupled. Hence it is possible to analyze the dynamics by considering two 
sets of equations, each of half the order as the whole system. This case is 
applicable to side-by-side or tilting proprotor aircraft configuration in a 
trim flight condition with no sideward velocity or turn rate. 

The equations of motion must be restructured in terms of symmetric and 
antisymmetric degrees of freedom. The motions of the right and left sides of 
the aircraft are given by respectively the stim and difference of the symmetric 
and antisymmetric degrees of freedom; 
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2: ~ — S[*uJL. 

These conversions apply to the rotor bending, torsion, glmbal, and inflow 
degrees of freedom, and to the rotor pitch control variables. The columns 
of the coefficient matrices are reconstructed according to these definitions 
of the S 3 rmmetric and antisymmetric degrees of freedom. The symmetric and 


antisymmetric equations of motion are obtained by operating on the rows of 
the matrices as follows: 

*= 



=. -^( tvH>r _ X (rurt^ltz 

Finally, the symmetric 
from 

and antisymmetric drive train motions are obtained 

•^5 - 







(from section 4.3.3); the columns of the matrices are reconstructed according 
to these equations . Now the symmetric and antisymmetric degrees of freedom 
are completely decoupled, and may be analyzed separately. The state vector, 
control vectors, and gust vector for the symmetric system are: 

X = [ f'' Vj JU X, 

dp xp^ -£p 

V 
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V|» — ^ Ss 

^ — L 3"^ 

The state vector, control vectors, and gust vector for the antisymmetric 
system are: 

^ C ^ ^ '^x Hv ^^3 

V= L i. y 
\ — £ s«. 3 

3= CvG.^-^ 


Note that it is necessary to designate symmetric and antisymmetric elastic 
airframe vibration modes. 

6.4 Matrices of Rotor Equations of Motion 

6.4.1 Inertial matrices for rotor equations.- The inertial matrices 
for the rotor equations of motion in the nonrotating frame are given below. 
For clarity, the superscript * denoting the normalization of the inertial 
coefficients has been omitted. 


- 368 - 
















sz- 




















- 371 - 





s -to 





s -k 













-S -k 

Pk ® 








^k»-®B 








i 




«l 

o 




1 

I 

O 








I 

o 





























' 3 - 









375 



















377 





































MA- w^- 





































































































385 


-2yQ. -2yQ„ -2yQ 
















































6.4.3 Aevodynamio matvioes for rotor equations in nonaxidl floi},.- The 

aerodynamic matrices for the rotor equations of motion In nonaxlal flow are 

given below. Note that each matrix Is a summation over all the blades , that 

Is, m*°l, ...,N. The notation C ■■ costp and S « slnip Is used In these 

in m 

matrices « 
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6.4.6 Inertial matrices for two-bladed rotor.- The Inertial matrices 

for the two-bladed rotor equations of motion are given below. The notation 

C = cosij^ and S = sinii) Is used, where ■= ip + mir. 
mm m 
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6. A. 7 Aerodynamia matriaee for tU)o-blcded rotor. ~ The aerodynamic 

matrices for the two-bladed rotor equations of motion are given below. The 

notation C = cosil^ and S = sintj/ is used, where ip = ip + mir. 

mm m 
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7. LINEAR SYSTEM ANALYSIS 


The flight dynamics and aeroelastic stability analyses (sections 5.3.3 
amd 6) result in a set of linear differential equations describing the air- 
craft motion, of the form 

A 2.^ A| X + A«x ^0 V 

where x is the vector of degrees of freedom and v is the vector of control 
and gust Inputs. The coefficient matrices (A 2 , A^^, A^, and B^) are either 
constant or periodic in time. The analysis of these linear equations 
proceeds as follows (see reference 37) . 

7.1 State Variable Form 

It is convenient to transform the equations to a standard first order 
form. Let MX be the number of degrees of freedom and MV the number of 
controls (dimensions of x and v) . Assume that MXl of the degrees of 
freedom are first order, that is have a zero column in the spring matrix Aq. 
Reorder the degrees of freedom so these first states occur last in the 
vector X : 



where X 2 are the MX-MXl second order degrees of freedom and x^^ are the 
MXl first order degrees of freedom. Then 

Ao =» La. o] 

where Aq has dimensions MX * (MX-MXl) . The equation of motion can now 
be written in the form 

X = Ax Bv 
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where 


k = 


E> = 



with the state vector 




' • \ 


There are MX2 = 2*MX-MX1 states; the top MX states are the velocities of 
the degrees of freedom, and the bottom MX-MXl states are the displacements 
of the second order degrees of freedom. 

7.2 Constant Coefficient System 
7.2.1 Eigen-ccnalysie.- The transient solution of x * Ax is 

X. ^ I ^ M » 

where and are the eigenvalues and eigenvectors of the matrix A, and 
q^(0) are constants determined by the initial conditions. A is the diagonal 
matrix of eigenvalues, and M is the model matrix with the eigenvectors as 
columns. The system is unstable if Re > 0 for any mode. 

The frequency of a mode is u “ 1 Im X | , and the natural frequency is 
(I) * [xj. The frequency in Hz is obtained by multiplying by n/2ir. The 

period is then T « 1/u) sec, with a> the frequency in Hz. The damping ratio 
is C “ -ReJyfxj (fraction of critical damping). The time constant is 
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T = -l/(fiReX) sec. The time to half-amplitude is then '^■^12 “ *693t (or time 
to double-amplitude for unstable modes) . 

7.2.2 Statio response.- The static response, obtained by setting 

• • 

X = 0 in X = Ax + Bv, is 

X = -A-'Bv 


7.2,3 Frequency response . The frequency response, obtained by setting 

iwt - lo)t . • * , j 

X » XQe and v » v^e in x « Ax + Bv, is 


Vo \\ 

where 

\A = - ( A-iuT‘ B ^ B 


The frequency response can also be obtained from the poles and zeros by 

''i * ’ 


where p are the poles (eigenvalues of A) and z are the zeros. 

7.2.4 Zeros.- In Laplace form the equation x « Ax + Bv becomes 


X — — (Av — 


By Cramer's rule then 

XL 



( A - s') 


where A* “ (A-s) with the 1-th column replaced by the j-th column of B. 
The poles are defined by det(A-s) “0, hence are the eigenvalues of the 
matrix A, as above. The zeros are defined by det(A*) “0. In A* the 
diagonal elements are all of the form (a ^^ -s) , except for the 1-th column. 
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By expanding the determinant of A* It Is possible to reduce It to a form 
with (a^j^-s) for all the diagonal elements: 


^ A'*' = k, 'kS* CA| —s') 


Then the zeros are the eigenvalues of the matrix Aj^. We have then 


V 


_ _ _ Atf A* ^ si y, 'h' (^' 

i AiJt:(A-s') ~ ‘ Ji^cA-s'i ' irCf-J) 


and the static response Is k 2 “= IIz/IIp. 

7.2.5 Transient response.- The general solution of x = Ax + Bv Is 

y = Me Mx. +•) A't 


Consider the case with zero Initial conditions at t » 0 and control with 

o 

time variation of the form v = v f(t) (f “ 0 for t < 0) . Then 

o 


y — 


M 


-t V; 


Uv M~‘ 6v. 


M { fiCti] M~' Sv, 


fdiere F {F^} Is a diagonal matrix depending on the eigenvalues and on the 
Input function f. The cases considered are 

a. step response, f = 1 

b. Impulse response, f ■■ i(t) 

c. cosine Impulse, f = l/2(l-cos^^) , 0 < t < T 

2Trt 

d. sine doublet, f “ si n ^ , 0 < t < T 

e. square Impulse, f=l, 0<t<T 
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f. square doublet, f = lforO<t<y, f = -l for ^ t ^ T 
The function F^(t) is readily evaluated for these and other inputs. 

7.2.6 RMS gust response.- Consider the RMS gust response of the 
system x = Ax + Bg. Here the only input considered is the vector of gust 
components g. The gust is model as a Markov process: 







where w is stationary white gausslan noise, with zero mean and correlation 

2 

For an RMS gust velocity of we have “ 2a^ /t^. In forward flight 

the correlation time is obtained from r^ = L/2V, vrtiere L is the gust 
correlation length and V is the aircraft velocity. The RMS acceleration 
can be obtained by including accelerometers in the model: 

where t is the accelerometer lag. For the acceleration of a particular 
a 

state we have 

i; -I- OL = 


so the row of has a 1 in the column, with the rest of the elements 

and the entire row of C equal zero. For body axis acceleration at the 

-a ° 

point r we have 






where the summation is over the rigid body and elastic modes (see section 
4.2.2). The matrix is here zero except for the rigid body and elastic 

airframe acceleration elements , and the matrix is zero except for the 

Euler angle velocity elements. 
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The aircraft, gust model, and accelerometers are thus described by the 
set of equations 


1 ") 

• 

A O 



m 

O 

cx 



1 

cx 


o 









O o ~ T<w 



OL 


W 


or 

^ — V =2. -4- & VM 


The correlation matrix 

^ 

is then the solution of the equations 

F 2. — o 

The solution is 

Z - •s- 

where M is the model matrix of F, and I is the diagonal matrix with 1 

Cj 

for the gxjst columns and zero elsewhere. Note that only the diagonal elements 
of the correlation matrix Z are of interest normally. Let 



Th^n with the summation k extending only over the gust columns of 



is the RMS gust response of the I-th state. 

7.3 Periodic Coefficient System 

Consider the system of equations x “ Ax + Bv, where A is periodic 
with period T, A(t + T) = A(t) . The transient solution can be written 

X, = ^ U; ft) 

L 

as for the constant coefficient system, but here the modes u^ are periodic 
functions of time. The eigenvalues are obtained as follows. 

The state transition matrix (Ji is calculated by integrating ({> = A4> 
over one period, t = 0 to T, with initial conditions «j)(0) = I. Let 
C = <|>(T); md let X^ and v be the eigenvalues and eigenvectors of the 
matrix C. Then the system roots are 

=. -=^ 5 =: 

(the principal part - a multiple of 2 tt/T can be added to the frequency) 
and the mode shapes can be obtained from 



The system is unstable if ReX > 0 (or jx^j > 1) for any mode. 

It is necessary to integrate the equations 4> “ f(t, <t>) ~ A(t)<j) 
numerically from t = 0 to t = T. Two methods are considered. The first 
method is a modified trapezoidal rule: 
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+«♦, = 4*«') 

* '4’v* ^ ^t>*t 4>«+| ■*" A« +V>'^ 

-- <'1— %*>v.y' ci-i- 

= Cl -I- ^ A,*+|') ( 1+ <|>v, 

= 4>v, -•- %('->■ ^ A„^.,')C^»,-*-A„^.,')<^>^ 

where <(>^ ■ <Kt^) , 4*^+]^ “ ^ ” ^n+1 ”^n* second method Is 

the fourth order Runge-Kutta method: 

+v» C ^ +^*f^ ■»“ ^ 

~ ^v» 4»v\ 

< 4>v* Ka") 

= Av.^1 ^+v, -«- -t k,') 
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